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Viscosity and the Extraordinary Heat Effects 
in Glass 


By Arthur Q. Tool 


In accordance with experience concerning the behavior of glass at temperatures within 
its annealing range, an equation is proposed that relates the various extraordinary heat effects 
to the inelastic deformability and to the degree of superheating or undercooling. By using 
this equation in connection with the thermal-expansion curves of a glass within its annealing 
range, certain constants that are related to the coefficient of viscosity and its changes with 
temperature and the degree of superheating or undercooling have been determined with 
reasonable results. Such results make it possible to estimate the inelastic deformability of 
a glass in its various conditions at all annealing temperatures and are, therefore, valuable in 
connection with problems that are encountered in the process of annealing glass. The 
apparent success achieved in applying the proposed equation to experimental data suggests 
that the concepts underlying this equation are fundamental and must be considered in any 
theory concerning the constitution of glass or that of any other extremely viscous liquid. 
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I. Introduction 


1. Superheating and Undercooling with 
Accompanying Extraordinary Effects 
in Viscous Liquid 


As already pointed out [1, 2, 3],' there is a 
possibility that the extraordinary heat effects 
(including the rapid expansion) of glass in its 
annealing range result because the high viscosity, 
which increases rapidly as the temperature de- 
creases, causes glass to become considerably 
superheated or undercooled even when the rates 
of heating or cooling are very low. Although the 
phenomena are related in some respects, the 
superheating and undercooling of glass are not 
those experienced when melting and crystalliza- 
tion, respectively, are delayed. Rather, they are 
effects that may be envisioned as possible even 
in the most fluid liquid at temperatures far above 
its freezing point if unattainably high heating 
and cooling rates could be employed. When a 
liquid is in equilibrium at any given temperature, 
there is an average statistical relation between 
its molecules that is unique, since any change in 
temperature, even if infinitesimal, will change it 
accordingly. Such a change in the average 
statistical relation, which corresponds to an 
equilibrium condition, presumably requires many 
processes or readjustments involving the mole- 
cules of the liquid. Although some of these ad- 
justments are completed almost instantaneously, 
time is undoubtedly an important factor in the 
development of many others, especially if their 
progress is affected by viscosity. If the time 
required for the full development of the sluggish 
readjustments, that are demanded by a given 
change in temperature, exceeds the time required 
to accomplish the temperature change, the liquid 
can be brought to a temperature in a state of 
nonequilibrium. That is, in the case of heating, 
the liquid can be superheated and, in the case of 
cooling, it can be undercooled. Moreover, if the 
cooling increases the viscosity so greatly that the 
development of the lagging processes ceases al- 
together, the liquid (unless reheated) will remain 
indefinitely in the undercooled state. It may even 
attain the amorphous solid condition if cooled 
sufficiently. 


1 Figures in brackets indicat at the end of this 


paper. 
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When a liquid glass reaches any temperaty 
a state of nonequilibrium, it might, on § 
thought, seem improbable that this under 
or superheated state corresponds closely j 
equilibrium condition at some intermediate 
perature between the temperature reached 
that at which the liquid was in equilibrium 
the cooling or heating was initiated. Doub 
there are many cases in which such a 
pondence is wanting, but it appears to ex 
most glasses if the tests are not too seared 
Considerable discussion of this point is incy 
in @ previous paper [1]. 


2. Equilibrium Curve of a Viscou 
Liquid 

When almost any glass is cooled and he 
through its annealing range at a low enough 
it is always practically in equilibrium on reac 
any temperature in that range. By deterni 
the changes in various properties as these tn 
ments proceed, equilibrium curves can be es 
lished. A rapid heating or cooling from 
temperature at which a glass is in equilib: 
causes departures from such curves. If 
departures are caused by sufficiently high rate 
heating or cooling and if the glass is returnel 
the initial temperature with equal rapidity 
no other lapse of time, it will be found that 
glass is still in practical equilibrium at 4 
temperature. However, if the initial heatix 
cooling is comparatively slow, the departures 
be less marked and a very rapid return of 
glass to its initial temperature will show th 
is no longer in equilibrium there. After 1} 
trials in which the same departure is always} 
duced, it will be found, nevertheless, that 
glass can be brought rapidly to some other « 
perature at which it is in practical equilibn 
This temperature will be higher or lower tha 
initial equilibrium temperature, according & 
heating or a cooling caused the departure. 
appears, therefore, that in certain respects ! 
has been no marked departure of the glass 5 
the equilibrium sequence that it follows on 
heating and cooling. 


3. Equilibrium Temperature 


In view of experiences such as those jus! 
cussed, it has been conjectured that glass, 
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or heated very rapidly in the annealing 
», behaves in many respects as a solid material, 
sas it behaves as a liquid when the rate of 
Further- 
_ it is inferred that the physicochemical 
ition or state of a glass is reasonably well 
only when both the actual temperature 
that other temperature at which the glass 
be in equilibrium, if heated or cooled very 
lly to it, are known. This latter temperature 
been termed the equilibrium, or fictive, 
perature of the glass, and a glass is under- 
or superheated according as the fictive 
perature is reached by the actual temperature 
gh heating or cooling, respectively. Ob- 
sly, the difference between these temperatures 
measure of the departure of the glass from 
ibrium. In the following discussion the 
ibrium temperature will be designated by r+ 
er to differentiate it from the actual tempera- 
T. 
s stated in previous papers [1, 2, 3, 4], the rate 
hich equilibrium is approached, whenever a 
rheated or undercooled glass is held at a 
ant temperature in the annealing range, is 
ly proportional to the departure of the glass 
equilibrium. Moreover, if the departures 
bf the same magnitude, an1 the temperatures 
lifferent, it seems that the rates of approach 
in addition, inversely proportional to the vis- 
y which increases exponentially with decreas- 
emperature. Thus, although equilibrium may 
pproximated in a few minutes at the upper 
of the annealing range, reaching a similar 
ee of approximation, even when the initial 
ure is no greater, may require many months 
e treating temperatures are near the lower 
of the range. 


hanges in Equilibrium Tempera- 
ep and Their Relation to Inelastic 

‘formability Under Elastic Molecu- 
Strains 


pecially when subjected to considerable loads, 
is, to an appreciable extent, still inelastically 
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deformable at temperatures far below its annealing 
range. At such temperatures and also at tempera- 
tures within the lower part of the annealing range 
there are indications that glass has ceased to be a 
truly viscous liquid and possibly has assumed a 
plastic nature [4, 5]. That is, whether the loads 
are of a purely mechanical or a physicochemical 
character, it appears that the number of flow 
planes decreases rapidly with the temperature 
unless the loads exceed “minimum”? limits, which 
presumably also vary with the temperature and 
condition of the glass. When flow in a plastico- 
viscous material is being maintained by stresses of 
either a mechanical or a molecular nature, the 
viscous flow that has been proceeding unhampered 
along certain localized flow planes may become 
completely blocked at any time because of changed 
molecular relations, or it may stagnate in a limited 
locality because conditions in the immediate or sur- 
rounding neighborhood have reduced or eliminated 
the load, at least, temporarily. Likewise, the 
resistance in a flow plane that has been blocked 
may break down because the flow in neighboring 
planes has unkeyed the block or developed stresses 
that exceed the breakdown limit. As a conse- 
quence, some readjustments and molecular proc- 
esses in a glass approaching equilibrium may not 
only develop much more rapidly than others, but 
some may also become completely blocked at low 
annealing temperatures. The discussion [1] of 
some peculiarities observed in the exothermic 
effects yielded by samples of a chilled glass, that 
were only partially annealed at low annealing 
temperatures, was based on this concept of the 
behavior of glass at low annealing temperatures. 
Furthermore, according to this concept, it is 
obvious that any theory based on the assumption 
that glass is a purely viscous liquid will not be 
applicable at temperatures below the annealing 
range. The effects caused by the development of 
a plastic condition are demonstrated in the follow- 
ing pages. However, such effects are simply modi- 
fications of the heat effects, density changes, etc., 
that would have occurred if the elastic molecular 
strains had relaxed because of viscous flow alone. 








II. The Approach to Equilibrium 


1. Differential Equations for the Rate 
of Approaching Equilibrium 


According to the introductory statements rela- 
tive to the rate of approaching equilibrium, 


dr/dt=Ky(T—1). (1) 


In this relation Ky is a factor that is inversely 
proportional to the viscosity as long as there is 
little change in the coefficient of elasticity, and 
dr/dt and T—r are, respectively, the time rate of 
change in the equilibrium temperature, 7, and the 
departure of this fictive temperature from the 
actual temperature, 7. Assuming the empirical 
relation suggested by Twyman [6] for the depend- 
ence of viscosity on temperature in the case of a 
glass in its annealing range, it appears that 


Kr=Ke™"*. (2) 


In this expression, k is Twyman’s constant, and 
K is obviously the value of K; when T=0. Thus, 
K is an extrapolated coefficient that relates to the 
inelastic deformability at 0° C. Experience 
indicates that K is approximately the reciprocal 
of the Maxwellian relaxation time for mechanical 
stresses, and it is, consequently, extremely small 
for glass at ordinary temperatures. This appar- 
ent similarity between the approach to equilib- 
rium and the relaxation of ordinary elastic strains 
first suggested that the relaxation of certain elastic 
molecular strains effects the approach to equilib- 
rium. 
On combining eq 1 and 2, the expression 
dr/dt=K(T—r)e™* (3) 


results. Even in the annealing range, eq 2 and 3 
yield only a very rough fit when they are applied to 
data procured onan annealing glass. For example, 
dr/dt will decrease somewhat more slowly with time 
in an undercooled glass and somewhat more rapidly 
in a superheated one than eq 3 allows. One pur- 
pose of this paper is to demonstrate this inade- 
quacy of eq 3 when it is applied to the extraor- 
dinary contraction and expansion and to the 
exothermic and endothermic effects of a glass that 
is being heated through its annealing range. 

That eq 2 and, consequently, eq 3 are not ade- 
quate is not surprising as the effect of a changing r 
on the viscosity is ignored. From experience it 
seems that the manner in which a property is 
affected by a change in r is very similar to that in 
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which it is affected by a change in 7. Py 
reason it has been assumed in previous cases! 
that eq 2 would be more nearly adequat 
had the form 
Kr=Ke?!*e™ 
and g and Ah were analogous to Twyman’ 
stant k. Equation 3 then becomes 
dr/dt=K(T—r)e™!*er™. 
By this means the above-mentioned inadeq 
in expressing the change in 7 at all ordinary a 
ing temperatures seem to be removed. Ho 
more precise data on the changes that are « 
in various properties and effects must be pme 
before a final decision on this point is rea 
In any case, it is obvious that eq 4 and 5 an 
inadequate as soon as glass ceases to be a} 
viscous liquid and assumes a definite pla 
viscous nature. That is, eq 4 does not take 
account the possibility that the coefficie 
inelastic deformability may depend on the 
Moreover, this equation may not be adeq 
express the changes in viscosity over wide ta 
ature ranges. 


2. Equations Modified for Linear 
ing or Cooling 

In many heating tests, 7 is increased ats 
stant rate that considerably exceeds that at 
r decreases at temperatures below the an 
range even if the initial r-value was high. 1 
in the plastico-viscous range of a glass, the 
tion between 7—r and T diverges compari 
little from a linear one during any given # 
which the heating rate is relatively high 
constant. The absolute value of the factor 
as it appears in eq 3 and 5 may be thought 
possessing some of the qualities of a mech 
load, especially in its relation to the ind 
deformability of a plastico-viscous glass. 
is, as this factor increases, the deformabilit 
creases. If the change in the plastic deform 
is an exponential function of the load, and! 
special heating condition mentioned above ¥ 
in any test, then eq 5 should be modified fe 
plastico-viscous range by introducing the ! 
e‘?’-Tv8’ In this exponential factor, 7” * 
sumably near the upper limit of the range in 
the glass is appreciably plastic. 
When a glass is being heated or cooled 
constant rate (dT /dt=+R, respectively) 
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being held at a constant annealing temperature, 
3 and 5 become 


T. Be dr/dT= + (K/R)(T—r)e™* (6) 
S Cases! 


equate MM = dr/dT=+ (K/R)(T—s)e™*er™, (7) 


if 7 and dr/dT are given, 7—r is definitely 

ermined according to eq 6, and the equation 

also easily solved in terms of exponential inte- 

ls. That is, if y=(T—r)/k and +2= 
Tk Lee 

naden /R)ke™, dr/dT= +2y and 


nary Bi( +:2)—Bi(+2,) =ye**, 8) 


How co dy-+yde—=(dz)/z. In this equation, 2 is the 
t are QMMMlue of « for some particular temperature, 7», 
be POSE which the glass was in equilibrium. (For the 
1S Te@dition that equilibrium does not exist at 7%, 
id 5 ani right-hand member becomes ye**—yoe**o. 
be aj viously, the negative sign before z is used in the 
te ple of cooling, and y is positive or negative 
ot takeMording as the glass is superheated or under- 
efficiel led. 
mn the ith respect to eq 7, assigning values to 7’ and 
acleq dT leaves T—r indeterminate, since, under 
ride te ain conditions, r may have, from a physical as 
as from a theoretical standpoint, either of 
ear Hao values,? depending on 7, the initial value of r. 
hough no simple solution for this relation is 
od at a@@meilable, the course followed by T—r during any 


rat at When y or (T'—r)/h, has the value +1, the product ¢-7/*dr/dT has only one 

ie anne the maximum possible. If, in addition, dr/dT=+-1, then g is at a 

oh oe and eT/'=Re/Kh. Furthermore, T=k In (Re/Kh) is the minimum 

ign. perature, T'», for which dr/dT' can be +1. That is, ye-*<e- if positive y 

3s, the ber greater or less than unity, and 7>k In (Re/Kh) if dr/dT=1 at the 
time. 


mparae ben ye-* has a positive value less than ¢~, it is obvious that y has a pair 
iven té ossible values, one greater and one less than unity. If 7 is the same in 
high cases, dr/d T' will also be identical and, moreover, {t is obvious that the 

y of », under such conditions, correspond to different heating curves. 
factor is, the larger value corresponds to a curve procured on a comparatively 
ed sample of glass (re<7\.), whereas the lesser corresponds to a 


ought obtained, under the same standard conditions of test, on a less well 
mech ed sample of the same glass (rs>T7..). Furthermore, if dr/dT=1 
he im hese values of », it follows that the common temperature, T, for the two 
e ine is greater than Ta. When the usual heating rate (3° to 6° C/min) is 
lass. oyed, it will be found that the condition y=<1 and dr/dT=1 usually 


bili when the sample has been moderately well annealed. Since the heat 
MADUN' EM (endothermic effect, for instance) is Just getting under way when 
formal 's unity, it often appears that the beginning of the endothermic (or 
and F d expansion) effect in a moderately annealed sample occurs at a somewhat 
’ . temperature than in either less or more thoroughly annealed samples. 
bove 8 mples of series of curves that present this appearance are found in figures 
fied for d8 of a previous publication {1}. 
or the chilled and annealed samples that were used in the expansion tests 


- the f are described in this paper, the temperatures at which dr/dT equaled 

T’ is y were estimated to be 586° and 578° C, respectively. When 7. was 
p= 5 puted by using the values found for k, A, and K (or R/K), the value 
ie 1D C was obtained. 


€q 8 were applicable, there would be no T. above which the maxima of 
endothermic effects must be reached, and consequently they would not 
eooled onfined to a limited range. Experience shows that the temperature 
ly) in a maxima occur is limited, even if the treatment of a glass 
‘4 : Mi) 
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heating or cooling is quite definite for any given set 
of initial conditions. Since (1/g)+(1/h)=1/k, as 
previously shown [4], eq 7 becomes 1/z—(h/k)dy/ 
dz=+ye~" if substitutions,» y=(T—r)/h and 
+2=(K/R) he™*, are made. However, from a 
practical standpoint, neither this nor any of the 
numerous other forms that can be obtained by 
various substitutions and transformations seem 
to possess advantages over eq 7. Although a 
solution of this equation is lacking, it as well as 
eq 8 can be tested by the use of thermal expansion 









































data and also by a study of their significance with 
respect to the endothermic and exothermic effects 
observed in glass. 


3. Relation of Equilibrium Tempera- 


ture to Changes in Length 
When 1, that is 7—ky, has been computed by 


means of eq 8 for various temperatures through 
the use of assumed or determined values for k, 
K, %, and R, the computed results may be com- 
pared with experimental data obtained on any 
property that is subject to the extraordinary 
changes observed when a glass is heated or cooled 
through its annealing range. For example, the 
linear thermal expansion curves may be used, 
although it is difficult to determine precisely the 
expansivities for those values of 7’ and + that are 
in that part of the practical annealing range where 
softening begins. It appears that no great errors 
are introduced by assuming that the linear expan- 
sivities (a for varying 7 and constant 7, and a 
for varying r and constant 7) are approximately 
constant throughout the annealing range and have 
values that can be determined with reasonable 
precision near the lower limit and sometimes even 
throughout the lower half of that range. Thus, 
it appears that the linear expansion per unit length 
through the annealing range is expressed approxi- 
mately by 


(L—Iy)/In=a(T—T) + a(r—1) (9) 


if ZL, is the length of the ylass sample when 
T=T, (a temperature that will usually be chosen 
at some point within or just below the annealing 
range and that may or may not be that for which 
z=, and r=7)). The use of eq 9 in conjunction 
with eq 7 is considered in the following section. 


3In making these substitutions, use is made of the identity e?/ser/4= 


eT/te(-T/) , However, when appearing in an exponent as in this case, 
T—r, although its magnitude is presumably proportional to certain molecular 
strains, cannot be considered as necessarily having an effect that is analogous 
to that of a mechanical load in its relation to inelastic deformations. 
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III. Application of Equations to Thermal Expansion Data 


1. Thermal-Expansion Curves 


In a progress report [2] made to the Glass Sec- 
tion at the Cincinnati meeting of the American 
Ceramic Society in 1942, the inadequacy of eq 
6 and 8 was shown by applying them to experi- 
mental expansion data. It was also pointed out, 
although not shown definitely, that much of the 
discrepancy between experimental and computed 
results could be eliminated by assuming that the 
viscosity of glass is affected by changes in r very 
much as it is affected by changes in 7. In the 
following, a résumé of that report and a fuller 
discussion of the effect of changes in + are 
presented. 

The expansion data used for this purpose were 
obtained on two samples of a borosilicate crown 
that was of the usual kind employed in ther- 
mometers that are graduated to temperatures that 
exceed 500° C. One of the samples was treated at 
a high temperature (approximately 850° C) and 
cooled rapidly in moving air in order to prevent r 
from falling below 650° C. In fact, the cooling 
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TEMPERATURE 

Figure 1.—Comparison of computed and experimentally 

determined thermal-expansion curves for a borosilicate 
thermometer glass. 


1 and 1’, Experimental curves for an annealed sample (circles) and a chilled 
sample (dots), respectively; 4 and B, tangents to curves 1’ and 1, respec- 
tively, at points near 390° C; EE’, the approximate thermal-expansion 
curve of the glass when it does not depart from equilibrinm because the 
heating and cooling rates are very low; C and C’, computed thermal-expan- 
sion curves (broken lines) corresponding to 1 and 1’. These curves were 
computed in accord with eq 7, and the departure of C’ from 1’ below 500° 
C results presumably because the equation is inapplicable while a glass is 
in its viseo-plastic stage. 
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rate of this treatment was such that the red 3 
of the small sample in a comparatively dark y 
disappeared in about 5 seconds. The other say 
was treated for several weeks at 500° C in one 
reduce r to that temperature. Each sample, 
sisted of three small pyramids of a type thy 
suitable for expansion measurements by the iz 
ferometric method and that makes very 
cooling possible after treatments at high » 
peratures. The expansion measurements 
tended from room temperature to the softe 
point near 620° C and were made by Jame 
Saunders according to the weighted-spacer met 
that he has described [7]. 

The portions of the resulting expansion cw 
between 260° and 620° C are shown in figuy 
and 2 (curves 1 and 1’). The relative placey 
of the two expansion curves (curve of the 
sample above that of the annealed) was determi 
by making them practically coincide at 4 
615° C, where both samples were just entering 
range of slight but noticeable softening. In 
range the r-values of the two samples were} 
sumably almost equal, and the values of 7—; 
doubtless not much greater than a degree, ss 
heating rate used in these tests was only 3 deg 
min. That is, within the range of apprec 
softening, it is difficult to maintain T—r at vi 
exceeding 1 degree unless the rates of heati 
cooling are exceptionally high. In other w 
if deformation of the samples had not taken p 
the expansion curves presumably would 
practically coincided with the equilibrim a 
from this point upward. This curve is represet 
by the line ZE’, and it intersects the expan 
curve of the annealed sample at 500° C bee 
up to that point there was no appreciable ch 
in the r-value of that sample. However, if! 
heating rate from this point had been made " 
low, the resulting expansion curve on the anned 
sample would have coincided approximately 4 
EE’ from 500° C until the temperature of ap 
ciable deformation was reached. Doubtles 
line (such as EE’) representing the sequence 
equilibrium conditions in terms of chang 
length should be somewhat curved; but int 
purely demonstrational presentation, relati 
minor details of this nature will be ignored for! 
purpose of simplicity. 
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TEMPERATURE 


pre 2.—Comparison of thermal-expansion curves deter- 
ined experimentally and computed in accord with 


ation 8. 


, B,and EE’, Reproduced from figure 1; 2and 2’ (large circles and large 
, respectively), correspond to curves | and 1’ and were computed on the 


mption that the necessary corstants were k=15 degrees, Kae —*"4, 
T,=560° C, and also that the initial r-values were 500° C and 653° C 


the annealed and chilled samples, respectively; 3 and 3’ (small 
Hes and small dots, respectively), similarly computed except that the 
tants k= 13.57 degrees, K=e~#4, and 7; =544.1° C were used. This 
that eq 8 is 


parison of experi tal and puted curves indicat 
tisfactory, probably because the viscosity varies as the equlibrium 
perature changes. 








[ethod of Estimating Expansivities 


etween 300° and 500° C the expansivity of the 
ealed sample was found to be about 0.066 
)/deg C, but below this range the average 
less than 0.063. Although the expansivity 
n severely chilled glass is usually increased a 
percent at atmospheric temperatures, the 
ansivity of the chilled sample also appeared 
be about 0.066 at 300° C. If a line coinciding 
the expansion curve of the chilled glass and 
ing this slope (0.066 w/deg) is extended to high 
peratures, it intersects EE’ at 653° C, which 
resumably the initial value of r in the chilled 


bple and also its value during the earlier steps 


the expansion test and until the equilibrium 
ition of the sample began to change appre- 
bly—as it appears to have done immediately 
ve 340° C. This temperature is about 150 
below the lower limit of the practical annealing 


ge of this glass, and is, therefore, far below the 
ge in which purely viscous flow is presumed to 


possible. 
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If the relative placement of the two curves is 
approximately correct, the difference (653 to 500 
deg) between the initial r-values of the two sam- 
ples produced a difference in length of 38.4 »/cm. 
This corresponds to a linear expansivity of 0.251 
(u/em)/deg C change in r. In the usual units, 
this value also corresponds to a volume expansivity 
of 0.75X10~* for a 1-deg change in +, and it is 
about 0.7 as large as the value found by testing a 
borosilicate optical glass. The sum, 0.251+ 
0.066=0.317, is presumably 10* times one-third 
of the volume expansivity of the glass when it 
behaves as a liquid within its annealing range— 
tbat is, when this particular glass follows its equilib- 
rium curve (line EE’). 


3. Determination of Equilibrium 
Temperatures 

























































Knowing the expansivities and the initial r- 
values, it is a simple matter to compute the 7 
value at any temperature, 7, from the expansion 
curves. Thus, for the chilled glass, r=653°— 
[0.066(7'— 300) — (Lr— Lgq)/Lr)/0.251 (see eq 9) and, 
for the annealed glass, the expression is the same 
except that 500° replaces 653° C. In figure 3 
(curve 1, circles, annealed; curve 1’, dots, chilled 
glass), the r-values so obtained are plotted against 
T. 
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Ficure 3.—Curves relating to changes in the equilibrium 
temperature as a glass is heated or cooled. 


1 and 1’ (Circles and dots, respectively), equilibrium temperatures com- 
puted in accord with eq 9 from the corresponding experimental curves 1 
and 1’, figure 1; C and C”’ (broken lines), equilibrium temperatures that 
correspond to curves C and C” of figure | and were computed by a step-by- 
step method based on eq 5; B, A, and D, tangents (dr/dT'=0) to equilib- 
rium temperature curves | and C, 1’, and C’, respectively, at points near 
300° C; EE’, trend of the equilibrium temperature (slope dr/d7T=1) when 
glass remains in constant equilibrium because the rate of heating (or 
cooling) is very low; T (parallel to EE’), tangent to curves 1 and C when 
degree of superheating (or T’—,r) reaches a maximum; V, one of a family 
of curves indicating conditions for constant viscosity; F, @, and H, approxi- 
mate courses that would be followed by the equilibrium temperature 
during cooling if the cooling rates (30, 37.5, and 40 deg C/sec, respectively) 
from comparatively high temperatures prevented the equilibrium tempera- 
ture from falling below 653° C when 7’ became zero. 


79 


4. Insufficiency of Equation 8 Which 
Relates Viscosity to Actual Tempera- 
ture Only 


In view of the temperature range involved and 
the form of the r-curve for the annealed giass, it 
is comparatively easy to select values for K and 
k that make it possible to compute, by means of 
eq 8, a curve that follows the experimental curve 
moderately well. By introducing the r-values so 
computed into eq 9 and thus transforming to ex- 
pansions again, the computed results shown in 
figure 2 (curve 2) were procured. To delay the 
beginning of the rapid expansion to the point in- 
dicated by the experimental curve, it was neces- 
sary in computing these r-values and expansions 
to assume that Ky, possessed the value unity 
(reciprocal of 1 hour) at a rather high tempera- 
ture (560° C). Furthermore, it was then neces- 
sary to assign a rather high value (15 deg) to k 
in order to prevent the rapid expansion from de- 
veloping too rapidly when 7—r reached the rather 
large values indicated by the experimental re- 
sults. As a result, when this difference became 
comparatively small at high temperatures, the 
rate of expansion was too low as a comparison of 
curves 2 and 1 (fig. 2) indicates. Moreover, when 
the values assigned to K and k were used to com- 
pute 2’ for the chilled sample, the contraction ef- 
fect, that accompanies a decreasing r, was dis- 
placed to temperatures that were far too high, 
since its range almost coincided with that of the 
rapid expansion of the annealed sample. As a 
matter of fact, such a result is to be expected 
because it was presumed in the development of 
eq 8 that the viscosity at any temperature was 
always the same, regardless of the value of r. 
Consequently, unless there is a large difference in 
the absolute magnitude of the values of 7—r in 
the two cases, the contraction of the chilled sample 
and the rapid expansion of the annealed sample 
will begin at about the same temperature. Also, 
the maximum rates of contraction and expansion 
in these cases will be attained at temperatures 
that are not materially different. 


5. Determination of Constants in Equa- 
tion 7 Which Relates Viscosity to Both 
Actual and Equilibrium Tempera- 
tures 


As differential eq 7 has not been solved explic- 
itly for + (except in the form of series suitable 


only for relatively narrow ranges in Tj) 
simple procedure followed above cannot be , 
in applying this equation to experimental 4 
However, if the slopes, dr/d7, are well en 
known at three well-separated points in the 
nealing range, the constants, g, A, and K/R, 
be determined. Having made such detery 
tions, it is then possible to compare compy 
values for the ratio (dr/dT)/(T—r) at va 
temperatures with the similar ratios that ea 
obtained by using the slopes (fig. 4) and de 
tures from equilibrium (fig. 5) that are both 
tained from curves 1 and 1’ in figure 3. U 
tunately, the relative error in the slope dete 
nations from these curves is large over mud 
the range of the extraordinary heat effects and 
some cases (especially throughout the rang 
the superheating of chilled specimens), the x 
tive errors in determining 7'—r are also n 
large. Nevertheless, it is considered that 
data obtained for dr/dT and T—r (figs. 4 ani 
are adequate for the purpose, because this py 
is purely a demonstrational presentation wit! 
pretense of high precision. 

The three ranges in which these data seeme 
be comparatively good were near the respect 
temperatures (1) at which dr/d7=1, (2) at wi 
this slope attained a maximum in the case of 
superheated annealed sample, and also (3) 1 
the temperature at which the slope reaches 
minimum (maximum negative value) in the « 
of the undercooled chilled sample. Letting p 
d symbolize slope (dr/d7) and departure 1 
equilibrium (7—r) and the subscripts ¢ ani 
indicate the chilled and annealed samples, re 
tively, the data in the first two cases were | 
Pe=1.00, T.=578.0° C, r—a=513.0° C, and é 
65.0° C; and (2), pe=5.37, T.=603.0° ©, 1 
578.0° C, and d,=25.0° C. In the third « 
the data at a point chosen near that of the m 
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Ficure 4.—Slopes of the equilibrium temperature 

and 1’, figure 3, at various actual temperatures. 
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were (3), Pe=—0.72, T,=542.9° C, r= 
° ©, and d-=—46.8. Use of the first and 
of these sets of data with eq 7 gives the result 
(r,—12)/(T.— T.)=2.185.4 On using this 
in connection with the second set of data, 
liows that g=19.78, h=43.22, and k=13.57. 
value of & is rather high, but it is by no means 


using the value found for h/g in connection 
the second set of data, advantage was taken 
e fact that p.=5.37 was presumably a maxi- 
», At such a point, dp/dT=0, and thus, on 
entiating eq 7, it follows that for this condi- 
h=(p+h/g)d/(p—1) and that p=(1+d/g)/ 
h). These equations also hold for a mini- 
nin p. The fourth set of data was chosen at 
int where such a minimum appeared to be 
was as follows: (4) p-=—0.73, T.—540.0°, 
591.6°, and d,.=—51.6°, and it was found 
the expression relating p and d was reason- 
well satisfied by these data without changing 


g. 

ving determined g and h, any of these sets of 
can be used to determine K/R and also 7,/k 
introducing two sets of data (such as the first and third presented 
Ny=(n—n)/(Ti—T)+(/(Ti—T,)) In (Pidalpads). In the above case 


d set of data was taken at a point such that p,/d; = pa/d; approximately. 
nuently the logarithmic term disappeared. 
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E 5.—Degrees of undercooling and superheating 
~r1) indicated by curves 1 and 1’ in figures 1 and 8. 
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if 7, symbolizes the temperature for which the 
glass, when in equilibrium (7'=7), has a relaxation 
time, 1/Kr, equaling 1 hour. That is, K= 
(pR/dje~*"e""=e-1*, Since R=180 deg/br 
(3 deg/min), the two values for K are K,=e~”, 
and K.=e~* if the second and fourth sets of 
data are used. The corresponding values for 7; 
are 545.6° and 543.5° C. Two such sets of data 
should give the same values for K and also for 
T, if eq 7 is applicable, regardless of the treatment 
received by the glass, and the fair agreement 
found here exceeds what should be expected in 
view of the unreliability of the d- and p-deter- 
minations. Besides, there is the probability that 
R varied more than was apparent during the tests. 

After these constants are determined for a glass, 
the relaxation times, 1/K,y, can be obtained for all 
conditions and temperatures in the annealing 
range. Even the viscosity coefficients, 7, can be 
estimated if the elastic modulus, E, is known within 
reasonable limits. That is, 7=E/Kr. 


6. Testing the Applicability of Equa- 
tion 7 


(a) First Method 


From eq 7 it is apparent that a straight line 
should result if In (d/p)—d/h is plotted against 
In (R/K)—T/k. That is, as long as the viscosity 
varies approximately as assumed (and also pro- 
vided that the relation between it and the extra- 
ordinary expansion and contraction is that 
assumed), there should be within the annealing 
range no departures from this straight line that 
are not ascribable to uncertainties in the data. 
For this test, the p and d values from figures 4 and 
5 and the A, k, and K (average) values given above 
were used. (Since the rate of heating was 180 deg/ 
hr, R=e§™ and In (R/K) =45.320.) 

The results of this test are shown by curve E 
in figure 6. Between 600° and 500° C the depar- 
tures from a straight line are reasonably small, 
except in a few cases in the range where the chilled 
glass became superheated and the procurement of 
data without very large relative errors was impos- 
sible. However, beginning near 500° C and 
extending to lower temperatures, there is a steadily 
increasing departure from a straight line by the 
results for the chilled glass. The sense of the 
departure indicates that dr/dt at temperatures 
below the annealing range has absolute values 
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Fiaure 6.— Test of the degree to which changes in equilibrium 
temperature satisfy eq 7. 


E, Tests of results (curves 1 and 1’ in figs. 3, 4, and 5) derived from thermal- 
expansion data (curves 1 and 1’ in fig. 1), ordinates on the right. C, tests 
of data derived from curves Cand C’ in figure 3 (displaced ordinates on the 
left). (Circles and dotes, annealed and chilled glass in both cases.) 





(at least when r— 7 is large) that are greater than 
they would have been if the deformability had con- 
tinued to decrease with decreasing temperature 
according to its manner of decrease within the 
annealing range. That is, the relatively distended 
glass of the chilled sample under heavy loads 
collapsed more readily at these relatively low 
temperatures than it would have collapsed if it 
had retained an almost purely viscous character 
with a viscosity coefficient that continued to 
decrease with 7 and 7 as in the annealing range. 
This result was expected, because experience in 
attempts at establishing equilibrium conditions 
at such relatively low temperatures led to the 
conclusion that the rate of collapse at these temper- 
atures is no longer roughly proportional to the 
load, as in the annealing range, but decreases so 
much more rapidly as the collapse progresses 
that—in the case of a glass approaching equi- 
librium at a constant temperature—the product of 
e'T-")/*/(7 —r) and d (T—r)/dt decreases rapidly 
(instead of remaining constant, as it would if the 
glass were purely viscous) as (r— 7’) becomes small. 
That is, when a glass is at a temperature such 
that it is definitely plastic, this product exceeds 
or falls short of Ke?* in magnitude according as 
the physicochemical “load” (as represented by 
(T—r)) is large or small. 


(b) Second Method 


As figure 6 indicates that the chilled sam 
not behave as required by eq 7 until 500°¢ 
exceeded, it becomes a matter of interest to; 
mine approximately how much lower than ¢ 
the initial r-value would have been if the 
had behaved according to that equation thn 
out the heating. Such a determination « 
made by laborious step-by-step computation 

If the temperature is constant, eq 5 « 
written in the form e~*dz/z=— KeT’*dt, in 
2=(r—T)/h. The solution of this for any iny 
of time, At, for instance, is Hi(—z;)—Ei(-4 
— Ke?’*at, in which z and 2, are the valus 
at: the beginning and end of the time int 
Consequently, the r-value reached by a treatz 
at J) is 1=hz,+7,. Assuming, now thai 
temperature is changed so sharply from J 
T)+AT at the end of this first step, that 
unchanged, it follows that e7,/* is chang: 
e'7%o-47/® and that 2 is changed from (7-] 
to z,—AT'/h for the second step in which the 
is treated for an interval, Aft, at an increased 
perature. If 2, is the z-value reached at the 
of the second time-interval, the r-value reach 
t2=hz,+7,+AT, and after n such step-by 
treatments, r,=hz,+7,+(n—I1)AT. 

By making AT/At=R, the rate of heatig 
degrees per hour, the above exponential int 
equation becomes 


Ei(— 2) = Ei(— 2) — (K/R)e™* AT. 


Using this equation in the step-by-step com 
tion and choosing some point (7), 79) on the en 
mental curve as the initial condition produ 
computed curve that passes through this i 
point and that approximates a curve that sa 
eq 7, and the smaller AT is made, the close 


approximation will be. Intervals ranging 
0.5 to 6 degrees were tried. In some ranges, 
intervals can be used without developing 
progressively increasing error too much, but® 
r is changing rapidly, intervals smaller tha 
least used would have improved the result 
tained. The intervals used in computing @ 
C and C” in figure 3 were 0.5 and 1 degree, 
T, and rt. were 516° and 608.4° C, for the 
sample, and 530° and 500° C for the ans 
sample. Consequently, the computed curve 
the r-curves derived directly from the 
mental data intersect at these points. 
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uted curves, C and C’, indicate in both cases 
the glass at temperatures above 516° C 
ved very nearly as eq 7 requires. 

hen the previously employed test—made by 
ing In ((d/p)e~*) against In ((R/K)e~*"*)— 
pplied to these computed curves, the results 
n by curve C in figure 6 were obtained. In 
of the fact that intervals no smaller than 
half degree were used, these results are as 
actory as could be expected, and they 
ate that the computed curves, C and C’ 
3), throughout most of their courses are not 
om curves that fully satisfy eq 7. 

















7. Evidences of Plasticity 





figure 3 the trend of curve C’ when compared 
at of curve 1’ at temperatures below the an- 
g range suggests that 653° C is about 28 deg 
er than the initial r-value required if eq 7 
expressed reasonably well the behavior of the 
pd sample at temperatures below as well as 
the annealing range. In accord with 
ously expressed surmises, this difference be- 
the actual and required initial r-values is 
mably the result of the failure of eq 7 to 
into account the increased mobility of plastic 
when subjected to very large loads. In 
onnection, curve £ of figure 6 indicates that 
p points that were obtained below 516° C, 
e T—r is almost proportional to 7, would 
been alined approximately with those 
ned above that point if the abscissas for the 
of these ranges had been computed by 
bs of the expression In (R/K)—T/k—(513— 
2.6. [tis to be presumed that the values of the 
ants 7” and g’ of this added term (513° and 
’in this particular case) depend, at least, on 
itial r-value and the rate of heating. In 
of previous mention that an added term 
tining such constants might be indicated 
certain conditions, the manner in which 
esults below 516° C depart from the basic 
pht line suggests that the deformability of 
in its plastico-viscous condition is an ex- 
ntial function of the load as well as of the 
rature. This possibility merits further 
tigation. 
is discussion of the possible characteristics 
lis particular glass at temperatures below 
C must not be taken as a suggestion that it 
tgoes some sudden and definite change at 
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that point. It is more probable that all changes 
that take place in ordinary glasses while heating 
or cooling do so gradually over a considerable 
temperature range. That is, this point for the 
upper limit of the plastic range has about the 
same significance as the softening point or the 
point chosen as the beginning of the rapid ex- 
pansion range. 


8. Use of Equation 7 in Determining 
Cooling Curves 


Since the magnitude of the cooling rate that 
produced an initial r-value of 653° C was a matter 
of some curiosity, approximations to cooling curves 
(as shown in fig. 3) were also computed by the 
step-by-step method ° and on the assumption that 
r had fallen to 653° C by the time that 7 had 
fallen to 550° C, below which there had been no 
further appreciable change in r. It was thought 
that the cooling rate might have been about 40 
deg/sec, since it appeared that the sample in 
cooling from 850° C had lost all redness in about 
5 sec. However, the cooling curve procured for 
that rate suggests a lower rate, because it indicates 
a crossing of the equilibrium curve at a tempera- 
ture slightly above 670° C rather than at 850° C. 
When a rate of 30 deg/sec was used in the com- 
putation, the resulting curve began at about 655° 
C to diverge from the equilibrium curve instead of 
continuing to approach it. That is, this rate 


+ Equation 8 can be used to determine directly a rate of cooling that is 
definitely lower than the required rate to produce a given +r-value 
(r’) when T becomes zero. When T'reaches zero, z= Kk/R and is very small 
unless R has a very low value. Consequently, e¢=1 and Ei(—z)=In 
(Kk/R)+y approximately. (Euler’s constant, y=0.5772+-.) Thus, eq 8, for 
cooling from an equilibrium temperature, 7, to zero degrees becomes — r’/k= 
In (Kk/R) +7 —Ei(—z0) which, since zo=(Kk/R)e™’*, may be changed so as 
to become [(To—r’)/k]--y =In z»— Ei(—zo). From this equation, ze is easily 
determined approximately when the left-hand member is known and ade- 
quate tables of exponential functions and integrals are available. Thus, R 
is known if 7) and the necessary constants of the glass are known. More- 
over, zo is so large, when the starting temperature, 7», is high and R is not too 
great, that Ei(—zo) is negligible. Consequently, forsuch cases R= Kkev++'/*, 
Using the values previously found for K and k, and assuming that equilib- 
rium was established when cooling began at a starting point near 850° C, it is 
apparent that z is large for this case and that R=19.9 deg/sec, approximately. 
Thus, a rate computed in this manner and known to be too low, is about one- 
half of the rate indicated by the step-by-step computation described above. 

Oceasionally it is desirable to estimate a minimum for the rate at which a 

glass in equilibrium at some temperature, 7», can be cooled to room tempera- 
ture without material change in r from ro= 7». In such a case, R must be so 
large that z» becomes very small. It then becomes obvious, on expanding 
Ei(—zo) in the above equation for determining zo, that Ty—r’ =kzo, approxi- 
mately, and that R= Kk*eTW*/(T)—r’). That is, if eq 8 were applicable, R 
is the minimum rate necessary to prevent ro—r’ from exceeding some small 
assigned value. Ifk is replaced by g, and K by Ke*/*, the result is that which 
would have been obtained if eq 7 had been used for the same purpose and 
on the assumption that r varied so little that e*/* could be treated as constant. 
The ratio of these results for R is g*/k? and the smaller result, at least, is much 
lower than required, but such computations aid in estimating the cooling rates 
required to obtain desired r-values. 
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was much too low unless a much higher rate pre- 
ceded it from 850° to about 660° C, for instance. 
Without using very small intervals in the computa- 
tion, it also appears that a rate of 37.5 deg/sec 


IV. Application of Equations to Endothermic and Exothen 
Effects 


1. Relation of Heat Effects to Inelastic 
Deformability 


As previously stated [4, 5), the manner in which 
strain is dissipated in annealing glass and the 
character cf the exothermic and endothermic 
effects, observed when chilled and annealed glasses 
are heated, suggest the gradual development and 
disappearance of a degree of plasticity in glass as 
cooling and heating, respectively, proceed within 
and especially below the practical annealing range. 
As for the expansion effects, so also in that of the 
endothermic effect, the degree of plasticity that 
develops within the useful annealing range can 
apparently be ignored. However, at unusually 
low and, for most purposes, impractical anneal- 
ing temperatures, it seems necessary to assume 
the existence of a considerable degree of plasticity 
in seeking an explanation for the progressive 
changes which appear in the exothermic effect of 
samples of a glass that, subsequent to a severe 
chilling from a high temperature, have been an- 
aealed for different periods of time at a very low 
annealing temperature. Curves showing such 
effects have previously been presented [1], and 
the progressive changes are believed to indicate 
that some molecular readjustments leading to 
equilibrium conditions develop rather readily at 
low treating temperatures, whereas others develop 
very slowly or not at all until the glass is sub- 
jected to higher treating temperatures. Such 
differences in the rates of development are easily 
explained if it is assumed that only a portion of 
the flow-surfaces normally available at tempera- 
tures in the annealing range are still open at 
lower temperatures for viscous flow while other 
portions are either entirely closed or are practically 
so unless exceptional loads are applied. 

Careful consideration will show that the heat 
generated by the flow, as the elastic molecular 
strains (which are excessive because of super- 
heating or undercooling) are relaxed, is probably 
negligible compared to the heat effects observed. 
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is somewhat too low. Thus, it seems thy 
probable average rate of cooling in the q 
treatment was slightly under 40 deg/see, 
curves F, G, and Z in figure 3. 
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These considerations then lead to the cond 
that the observed heat effects are the result 
rapidly changing specific heat as the equilij 
temperature changes. 


2. Experimental Conditions in Deg); 
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For understandable reasons, the actual 
tions that exist in obtaining curves that g 
exothermic and endothermic effects are : 
very different from the idealized condition 
sumed in discussing the methods of test. § 
erally, it is assumed that the test sample i 
small that any temperature gradients with 
while heating or cooling at a considerable 
can be neglected. Actually, this conditia 
seldom approximated and, as a result, the obs 
heat effects are blunted and spread over a tem 
ature range that is somewhat too broad. 

In tests of this kind the well-known differe 
thermocouple method is commonly used. | 
hot junction of this couple is in the test sam 
whereas the other is either in the wall that 
plies the heat to the heating sample and 
supposedly conforms to an isothermal surface 
it is in a so-called neutral body that is ench 
in contact with the sample, by the same is 
mal wall. Also, supposedly, there is no tru 
of heat from one of these junctions to the 
that would appreciably affect the magnitu 
the heat effects that are sought. At best, | 
ideal conditions are only approximated and, 
result, any improvement in the test condi 
usually accentuates the effects as they are % 
by heating curves. 
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3. Relation Between Heat Effects @@,, ;), 
Changes in Equilibrium Temperat@ erin, 
ily | 


If T is the temperature of the test sample 3! 
instant and if the simultaneous temperatut 
the isothermal wall is 7’,, the rate of heat (H) 
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to the sample is dH/dt=A(T,—T). Under 
ized conditions, the factor A in this relation 
e product of the area and the so-called external 
juctivity of the contact between sample and 
If the sample weighs 1 gram, the rate of 
heat absorption equals this rate of supply and 
'dT/dt+-odr/dt for glass. In this expression, 
the heat absorbed per degree increase in 1, 
s is the normal specific heat. Both of these 
ficients, like the expansivities, are functions 
rand JT. Thus, s’’=s(1+67+r)=s'+8r 
res as a first approximation for the change in 
specific heat with these temperatures. In the 
od form of this approximation, s’ is pre- 
nably independent of + and is, consequently, 
same for all samples of a glass whether chilled 
annealed. The normal specific heat of a mate- 
represents (in the case of heating) both the 
vy required to increase the thermal agitation 
l that expended in expanding the material 
inst the so-called intrinsic pressure that results 
m the molecular attractions. Increases in r+ 


e the intrinsic pressure somewhat, since a 
aration of the molecules usually reduces their 
tual attraction. Consequently, o represents 
ly an energy of separation and 8 is normally 
ative. 
By equating the heat supplied to that absorbed 
i introducing the first approximation for s’’, the 
tion J,—T=(R/A)(s’+88r+oedr/dT) is ob- 
ned and, for the neutral body, which pre- 
mably has no peculiar characteristics in the 
perature range of the tests, the corresponding 
ation is 7,—T7,=(R,/A,)s,’ if R and R, 
the heating rates of the sample and neutral 
y, respectively. It is helpful in practice if 
s’, and A,=A approximately. In that case, 
follows that 
—T,=(R/A) (s’(—1+B,/R) —8B(r—12)— 


odr/dT)=G. (11) 


this equation, r,, the r-value for a very well 
healed sample, is introduced because r never 
h be reduced to zero in ordinary glasses. 


in the following discussion, @ represents the 
perimental results found for T—T,. Ordi- 
ily these experimental results are greatly 
ected by gradients, and consequently seldom 
proximate the values that 7—T, should have 
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under ideal conditions. For example, in testing 
a well-annealed sample, 7—T7', should approxi- 
mate zero throughout the range between the tem- 
perature at which a steady heating condition is 
first established and that at which the endothermic 
effect begins. This is obvious because R,=R, 
T=T., and dr/dT=0. However, G may have a 
fairly large value that is either positive or negative, 
depending on the characteristics of the heating 
furnace and the manner of packing the sample 
and neutral body. When proper care is taken, 
these disturbing gradients do not change materially 
between the time a steady condition is first 
developed and that at which granular samples 
show a considerable degree of sintering. As a 
disturbing degree of sintering seldom takes place 
below temperatures that are well above those 
required for the completion of the endothermic 
effect, a second steady condition of heating is 
usually developed. In these ranges of steady 
heating, the G-values change very slowly and only 
because A, s’, and o change slightly as 7 (and 
also r+ when the second steady condition is 
reached) increases. Thus, when the results for G 
are plotted against 7, the curves show two dis- 
tinct levels, one preceding and the other following 
the heat effects. If these levels are represented 
by G, and G,, repeated tests on identical samples 
show that G,—G, does not vary greatly if proper 
precautions are taken. For this reason it will be 
considered that G properly represents T— 7, plus 
an error that is practically constant and can be 
neglected at least as long as only differences in the 
determinations of G for a single test are con- 
sidered. 

When several samples of a glass are tested, they 
will all have attained approximate equilibrium 
when the second level is reached and the several 
curves can properly be made to coincide approxi- 
mately in this region, although the G,-values may 
be quite different because of the errors just men- 
tioned. Because the variation in G,—G, is usually 
small in the case of identical samples, the initial 
levels of the curves will then almost coincide 
unless there is a wide spread in the initial r-values 
of the several samples. That is, when the spread 
is not large, the effect of the differences in the 
initial 7-values is often masked by the errors in 
determining G,— G@,, but this masking is insufficient 
to hide the effect of differences as large as those 
found between the r-values of severely chilled and 
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Ficure 7.—Ezothermic and endothermic effects. 


27 and 38, Observed effects on heating chilled and annealed samples, respec- 
tively, of a borosilicate optical glass; 1 and 1’, indicate approximately the 
heat effects which, according to curves 1 and 1’ in figure 1 and curves 1 
and 1’ in figure 4, should be found for annealed and chilled samples of the 
borosilicate thermometer glass if the decrease in specific heat with increasing 
equilibrium temperature were insignificant. The separation between 
curves 27 and 38 below 400° C is mainly the result of a comparatively low 
specific heat when r is high. 





well-annealed samples (fig. 7, curves 27 and 38). 


(See also ref. 1, fig. 8.) 

The significance of the difference between the 
initial levels of these curves is at once apparent 
from eq 11 when it is considered that dr/dT is 
approximately zero and unity in the ranges of the 
first and second levels, respectively. To make this 
significance more apparent, it will be assumed that 
the change in A with temperature is negligible. 
In that case, G,—G,—(R/A)(sB(12—1)+<¢), since 
R,/R=1 at both levels. Thus, by neglecting the 
relatively small 8-term, the difference between 
the levels is approximately Re/A, but as £6 is 
negative, there is some increase in G,—G,, as the 
initial r-value, 7;, is increased and thereby causes 
the B-term to decrease. The r-value, 7, is that 
of the temperature at which all of the various 
curves, when plotted to show the results of a series 
of tests, are brought into coincidence in the range 
of the second level. Unless R is very low or the 
glass has been chilled very severely, 7, is greater 
than Ti. , 

When dr/dT7 reaches its maximum negative and 
positive values in the temperature range between 
the first and second levels, the heating curves (as 
plotted) attain, respectively, their maxima (exo- 
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thermic effects), which are observed in ¢j 
glasses only, and their minima (endothe 
effects), which are very pronounced in all anne 
glasses. During the endothermic and exothe 
effects, (R,/R)—1 takes on small positive 
negative values, respectively. Consequently, 
both cases, the effect of the s’-term, eq 11, j 
detract somewhat from the effect of the o-te 
However, as when developing the approxin 
relation, G,— G,—Ro/A, the s’- and 8-terms, by 
relatively small, may again be neglected for g 
plicity. Then, if the values obtained for 7- 
at the peaks of the endothermic and exothen 
effects are designated by G, and G,, respectiyé 
it follows that G,;—G,=— (Re/A)(dr/dT).,, » 
G.— G, = — (Re/A)(dr/dT) max roughly approxim 
actual conditions. As estimated from curve 
and 38 in figure 7, it appears that the three dif 
ences, G,—G,, G,—G,, and G,—G, are about | 
1.0, and —8.0 deg C, respectively. Consequent 
(dr/dT ) min = — 0.55, and (dr/d7)max= 4.4. 
compared to the respective values —0.73 
5.37 found from the expansion curves of the th 
mometer glass, these values obtained from 
differential heating curves of the borosili 
optical glass appear rather low. However, 
result is to be expected since experimental cou 
tions prevent the heating curves from showing! 
full magnitude of the heat effects. 

If the full magnitude of the endothermic ef 
and the corresponding change in + could be © 
determined (especially in the case of a thorou 
annealed glass), it seems that it should be possi 
to determine the ratio s’’/c, since the relat 
between the relative cooling of the glass and 
change in + should be expressed approximatd 
by 8’'(G,—G,)=e(r4—17;), in which 7, is the m 
nitude of + when the relative cooling reaches 
peak. A rough estimate, which could be 30 » 
cent or more in error, gave the result s’’/o=14) 
In arriving at the estimate of 140 deg as the vw 
of r,—7;, it was cons.dered that the initial equi 
rium temperature, 7,, of the sample that yiell 
curve 38 was probably soraewhat below 440’ 
and that 7, was not far below 600° C, because! 
difference between the actual and equilibm 
temperatures is not great when the endothem 
effect reaches its maximum. 

As the heating rate in obtaining the heat 
curves under consideration was about 6 deg 
it follows that «/A=(G,—G,)/R=1.8/0.1 a 
from this result and the above ratio, that s”/4 
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Thus, for any probable value of s’’ at tem- 
atures in the annealing range, the order of A 
bears to be in the neighborhood of 10~*. This 
h approximation is not a wholly improbable 
ue, since A (instead of possessing its idealized 
ificance in these tests) was concerned with the 
duction of heat from some undermined surface 
he furnace through a wall that was of uncertain 
a and was composed of a thin layer of burned 
bestos paper and an undetermined thickness of 
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d for ¢ nular glass or alundum. However, this uncer- 
for T. ‘Mapaty is of minor significance, as the chief purpose 
xothe ntroducing the very rough approximations, that 


re obtained as described in the immediately 
eding paragraphs, is to give a clearer picture 
the significance of the curves that result from 
rather complicated type of test which, although 
ellent in a qualitative sense, does not lend itself 
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4. Wa Estimated Heat Effects Based on 
0.73 al Equilibrium Temperatures Obtained 
the it From Thermal Expansion Curves 

bs ~ Besides the reproduced curves 27 and 38 (which 
ever. mr’’ obtained by tests in which R=6 deg/min), 
tal conmere 7 also contains curves obtained for com- 
owing ison by plotting the data, shown in figure 5 for 


dT, according to the manner and scale used in 
ptting the experimentally determined differential 














“vty ting curves. Although the two glasses were 
sail like, the treatments that they had received 
7 pen pre not too dissimilar. Consequently, it seemed 
- relatiqame’’ their observed exothermic and endothermic 
; and e's should have had about the same magnitude 
simedl the heating rates had been the same. To make 
mana —G@, and G,—G, approximate 1 and —8 deg, 
mete pectively, for the synthetic curves (right-hand 
ee ile of fig. 7), it was necessary to multiply the 
ait lues of dr/dT by 1.48. This factor corresponds 
the va Re/A, since dr/d T is zero in the range of the 
endl t and unity in that of the second level. How- 
t viele when the difference in the heating rates is 
- 460" nsidered, it is seen that the factor in this case 
aus a” very little definite significance, since halving 
sila heating rate changes the magnitude of the 
athe at effects appreciably. Nevertheless, a survey 
the heating curves obtained on a number of 
re linary glasses indicates that the average 
degia agnitude of G,— G, is near 1.5 deg. 
Te In preparing curve 1’, the $8(r2—7;)-term was 
ld lected. As a result, G,—G, is the same for 
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both synthetic curves, 1 and 1’ (fig. 7). That is, 
the decrease in s’’ (at low temperatures) as 7, 
is increased was ignored. In all other respects 
the two sets of curves, experimental and synthetic, 
are quite similar in form. 

In figure 8 the synthetic differential heating 
curves 2 and 2’ (for chilled and annealed samples, 
respectively) correspond to the curves with similar 
designations in figure 2, and are replicas of those 
shown before the Glass Section at the previously 
mentioned Cincinnati meeting of the American 
Ceramic Society [2]. As then, these curves are 
now presented to show that the form of the 
exothermic effect, if the viscosity of glass de- 
pended only on the actual temperature, would 
not be that of the experimentally determined 
effect. Also, it is apparent that, contrary to 
experience, the exothermic and endothermic effects 
of a chilled sample would then occur exactly 
within the temperature range of the endothermic 
effect of a well-annealed sample. In plotting 
these curves, dr/d7=-+zy (see eq 8) was also 
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Ficure 8.—Comparison of computed exothermic and 
endothermic effects. 


land 1’, (Large and small circles, respectively) same as 1 and 1’ of figure 7; 
Cand C’, computed effects in accord with computed curves C and C” of 
figures land 3. The failure of eq 7 to account for the plasticity of glass below 
the annealing range presumably causes the divergence between C’ and 1’ 
below 500°C. 2and 2’, (broken lines) computed heat effects in accord with 
curves 2 and 2 in figure 1. The wide divergence between these curves and 
land 1’ presumably is caused by the failure of eq 8 to account for both the 
plasticity of the glass below 500° C and the decrease in the viscosity as 
the equilibrium temperature increases. 
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multiplied by the factor 1.48 to bring them into 
conformity with the other curves of figure 7. 

Curves C and C”’ correspond to the similarly 
designated curves in figures 3 and 1. These 
curves were obtained from the values of dr/dT 
that were procured by the previously described 
step-by-step computations, and they show a 
definite approach to the form and relative place- 
ment of the experimental curves. They also 


closely resemble the synthetic heating cur, 
(fig. 7) that were derived more directly from 
expansion data, except that much of the hy 
developed in the plastic range is not represent 
That is, the dependence of the viscosity on + be 
broadens the exothermic effect and displaces ; 
somewhat toward temperatures below those , 
the endothermic effects of annealed glass 
even if plasticity is not taken into account. 


V. Speculative Discussion 


In addition to the external pressure which is 
comparatively quite small, there are two pressures 
to be considered in relation to liquids, (1) the 
thermal pressure arising from the collisions of the 
thermally agitated molecules, and (2) the so-called 
intrinsic pressure arising from the mutual attrac- 
tion of the molecules. On the average, these 
opposed pressures are equal when a liquid is in an 
equilibrium condition at any temperature. Cool- 
ing or heating a viscous liquid rapidly from a 
temperature at which equilibrium has been reached 
disturbs the balance between these pressures 
because the liquid is then no longer in an equilib- 
rium condition at the temperature reached and 
may be said to be in a condition of “distension” 
or “compaction,” respectively. If the change in 
the thermal pressure during a heating or cooling 
is always greater than that in the intrinsic pressure, 
it follows that (until a balance is reached between 
the pressures) the intrinsic pressure will gradually 
compress the relatively distended liquid after an 
undercooling and that the thermal pressure will 
gradually expand the relatively compacted liquid 
after a superheating. 

The work expended when such a compacted 
liquid with a density, D, expands at a constant 
temperature and against the intrinsic pressure, 
P, equals the mechanical equivalent of the heat 
absorbed because of the expansion. Thus, 
Pdv=JDedr. From the previous discussion of 
expansion tests, the result dv/dr=0.7510-* was 
obtained for the glass tested. For the same glass, 
D=2.5 approximately, and according to the 
previous very rough estimate, =8''/17.5. 
Although s’’ for this glass is unknown, it appears 
from a summary [11] of the results of a number of 
investigations that the specific heat for the 
average glass at temperatures in its annealing 
range is somewhat in excess of 0.2. On the basis 
of these rough values, P=1.610" dynes/cm., 


That is, the intrinsic pressure of glass in its anne 
ing range appears to exceed 1.5 10* atmospher 
and, as previously pointed out [3], this is the ss 
order of magnitude as that of values that he 
been mentioned as possible for the intring 
pressure of water [8]. If it becomes possible 
obtain really dependable values for the intring 
pressure of glasses and other viscous liquids } 
this or other means, the values will be of conside 
able interest, since this theoretical pressure | 
closely related to surface tension and to certai 
other properties. It is also to be hoped th 
speculations of this sort may help in those theori 
that are based on X-ray studies of glass. Thi 
hope rests on the assumption that the ve 
excessive compactions and distentions found } 
extremely viscous liquids when superheated an 
undercooled are in some comparatively simp! 
manner concerned with the character and size 
the larger molecules in these liquids. 

Such molecules not only add to the viscosity anf 
thus hamper the readjustments of smaller molecu! 
to temperature changes but, when present i 
relatively large numbers, probably form com 
paratively extended frameworks [9] (supermol 
cules) that also become too compacted (or dis 
tended) whenever the temperature is rapidl 
raised (or lowered) and that are delayed in co 
recting this incompatibility by the viscosity am 
plasticity to which they presumably contribu 
so greatly. As previously pointed out [3], chav 
or ring-like molecules [10] and some other form 
of extended molecular structures that often hav 
been discussed as possibilities (for references 
see [10, 11]) might easily be subject to suc 
incompatibilities whenever there is a temperat 
change. 

Within the annealing range, these extendet 
lattice-like structures formed from simpler mol 
cules are presumably unstable in the sense tha 
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y (individually but not on the average if 
jibrium exists) are being continually altered 
form and size by the thermal agitation of the 
mponents. However, the rate of alteration is 
sumably so low that the structures not only 
ange comparatively little during a rapid heat- 
» or cooling over even considerable temperature 
ervals but also persist in a compacted or dis- 
nded state for appreciable periods after the 
dof the heating or cooling, respectively. 
4s the thermal agitation decreases in a cooling 
nuid, the persistence of such structures pre- 
mably increases rapidly and the rate of increase 
pubtless parallels that at which the viscosity 
creases. Furthermore, the structures gradually 
ome permanent and probably attain a con- 
erable elasticity of form if the cooling is 
tended to temperatures considerably below 
e practical annealing range. At this stage 
be liquid is overly distended and the intrinsic 
essure exceeds considerably the thermal pres- 
Consequently, much of the contraction on 
bling consists of something akin to an elastic 
pression in which the structural members 
jer various excessive elastic distortions or molec- 
lar strains. 
As the thermal pressure decreases toward zero 
n continued cooling, the elastic resistance of the 
ructures becomes a greater and greater factor 


in balancing the intrinsic pressure. However, it 
seems that their elastic coefficients should also be 
increasing quite rapidly at very low temperatures. 
In that case, the increase in compression per degree 
decrease in temperature decreases, and the in- 
crease in the elastic coefficients may even be 
sufficient at low temperatures to cause the strains, 
that have been building up, to diminish. Further- 
more, the expansivity of the vitreous solid may 
become negative if the expansion caused by this 
diminution in strain more than offsets the con- 
traction which accompanies the decreasing thermal 
pressure. That is, many, if not all, vitreous 
materials may yield expansivity curves that, in 
this respect, resemble the expansivity curve of 
vitreous silica at very low temperatures [12]. 

The compacted condition as caused by rapid 
heating can be produced to an appreciable degree 
only in the annealing range. In this condition 
the excessive thermal pressure may also be re- 
garded as the cause of elastic molecular strains. 
That is, the excess of the thermal over the intrinsic 
pressure is balanced by the elastic resistance of 
molecular structures to distortion. Even after 
equilibrium is established, molecular strains are 
doubtless present and continually varying about 
an average as one after the other of the pressures 
becomes temporarily the greater in limited elements 
of the whole volume of liquid. 


VI. Conclusion 


Practically speaking, glass at temperatures 
ithin and above its annealing range is a purely 
scous liquid that is easily undercooled or super- 
eated because its viscosity is so very high. 
owever, at temperatures far below its annealing 
uge, glass behaves as an elastic solid although 
here is an intermediate range in which it possesses 
¢ characteristics of a plastico-viscous material. 
Although the condition of an undercooled glass 
snot one of equilibrium at the actual temperature 
{that glass, it usually corresponds to one of 
jilibrium at some temperature within the 
mealing range and, if the actual temperature is 
ed quickly enough to that point, the glass 
il be practically in equilibrium there. A cor- 
sponding relation exists between the condition 
fa superheated glass and some equilibrium 
mperature within the annealing range. The 
gree of undercooling (or superheating) is 
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measured by the difference between the actual 
and equilibrium temperatures. 

The properties of a glass are affected not only 
by changes in the actual temperature but also by 
changes in the equilibrium temperature. Thus, a 
glass expands when the equilibrium temperature 
increases as well as when the actual temperature 
rises, and lowering either temperature increases 
the viscosity very rapidly. 

Ordinarily, the equilibrium temperature changes 
at a perceptible rate only when the actual temper- 
ature exceeds the lower limit of the annealing 
range; but, if the degree of undercooling is very 
high as in quenched glasses, it decreases at ap- 
preciable rates even when the actual temperature 
is in the immediately lower range within which 
glasses take on a plastic character. 

The rate at which the equilibrium temperature 
changes is proportional to the degree of under- 
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cooling (or superheating) and also to the inelastic 
deformability. Moreover, for the same departure 
from equilibrium, the rate increases very rapidly 
as the inelastic deformability is increased by 
raising the actual temperature. As this deform- 
ability also increases as the equilibrium temper- 
ature increases, the manner of approaching equili- 
brium from an undercooled condition differs 
somewhat from the manner of approach from a 
superheated condition at the same actual temp- 
erature. 

In approaching equilibrium from an _ under- 
cooled condition, a glass contracts and evolves 
heat, whereas it expands and absorbs heat in 
approaching equilibrium from a superheated con- 
dition. The rates at which these heat effects 
develop depend on the rate of approaching equili- 
brium and are therefore also controlled by the 
inelastic deformability. 

In view of the foregoing conclusions, an equation 
has been proposed that relates these various rates 
to the inelastic deformability. It is found that 
this equation is applicable to experimental data 
on the extraordinary expansion and contraction 
effects and on the exothermic and endothermic 
effects of glass in its annealing range. As the 
equation is based on the concept that glass is a 


purely viscous liquid, it is not applicable to data 
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An Arrangement with Small Solid Angle for 
Measurement of Beta Rays 


By Leon F. Curtiss and Burrell W. Brown 


An arrangement using a Geiger-Miller counter witn small aperture with a radioactive 
source at some distance from it is described for counting beta particles. Sources emitting 
only beta rays with maximum energies above 1 million electron volts of the order of 1 milli- 
curie can be measured in terms of the disintegration rate from the observed counting rate 
and the solid angle as calculated from the dimensions of tne apparatus. An independent 
check of the arrangement shows that this can be done reliably. A suggestion for discarding 
the use of the curie and substituting a unit consisting of 10° disintegrations per second to 
be called the “‘rutherford” is made. The curie is properly applicable only to members of 


the radium family. 


I. Introduction 


The development of improved methods for the 
broduction of radioisotopes, resulting in an increase 
n the use of such radioactive materials, has led to 

demand for better methods of measurement. 
nthe simple case where a radioelement emits only 
beta particles of fairly high energy, as in the case 
i P®, it should be comparatively easy to measure 
he strength of radioactive sources in terms of dis- 
itegrations per second, even for fairly strong 
wurces. This has not yet been accomplished, 
however. Disagreements of the order of 300 per- 
ent have come to the attention of the authors. 
The need for reliable methods of measuring the 
trength of sources of P® is accentuated by its 
hewly developed use in the treatment of leukemia. 
tis of considerable importance that the amounts 
bf P* as measured by various laboratories should 
be in reasonable agreement. Otherwise, clinical 
tsults obtained in various parts of the country 
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from the use of this element cannot be compared. 

The writers have undertaken to develop an 
arrangement, comprising a Geiger-Miller counter 
and source holder, with a sufficiently low solid 
angle that sources of the order of 1 millicurie can 
be counted in it directly. The number of disin- 
tegrations per second will then be obtained from 
the observed counting rate and the solid angle 
determined from the dimensions of the arrange- 
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II. Geometrical Arrangement 


Figure 1 is a diagram of the counter, chamber, 
and source holder built to test the effectiveness 
of this method of measuring beta-ray sources. 
A brass tube approximately 60 cm long and 4-cm 
internal diameter constitutes the chamber. In 
this a number of baffles, B, are mounted with a 
movable shutter, Sh, of 3-mm brass at the end 
where the source, S, is mounted. This shutter 
permits a measurement of background with a 
source in place. The opposite end of the chamber 
is closed by a steel plate, P, 6 mm thick with an 
aperture, A, 4 mm in diameter in the center. A 
counter, C, is mounted concentric with A, which 
is covered by a mica window of approximately 
2 mg/em*®. The chamber can be evacuated during 
observations. The construction of the source 
holder is shown in figure 2. It is made of Plexiglas 
and designed to reduce back-scattering from 
various parts. 

The most important spurious effect to be antici- 
pated in this arrangement is small-angle scattering 
from the walls of the tube. The very small solid 
angle used makes this a serious difficulty. Con- 
sequently, a series of studies was made of the 
effect of baffles of various materials, at different 
locations in the tube, on the counting rate. These 
studies revealed that Plexiglas and soft-wax 
baffles were equally effective in producing a 
definite minimum counting rate when the source 
was in full view of the aperture. 

Typical results are shown in table 1 for equally 
spaced baffles, consisting of disks 4 in. thick 














Ficure 1.—Arrangement of counter, vacuum chamber, and 
source holder. 





Fiaure 2.— Beta-ray-source holder 


closely fitting the internal diameter of the » 
with a central hole % in. in diameter, 
diameter of the source was 1 cm. 

The optimum diameter of the circular ope 
in the baffles was found to be slightly grew 
than the diameter of the source. Small cha 
at this approximate diameter had no effect on 
counting rate. 


TaB_e 1.—Relative scattering for different arrangeme 
of baffles 








Relative 
counting 
rate 

















aperture distance from 60 to 40 em by mounti 
the source on a rod extending 20 cm into the tub 
The counting rate observed under these conditi 
agreed within less than 1-percent with 
calculated from the ratio of the calculated » 
angles and the counting rate observed at 60 a 
From the actual dimensions of the apert 
diameter =4.06 +0.01 mm, and _ source 


aperture distance equals 60.0 +0.2 cm, the sil 


angle divided by 4x is computed as 2.861" 
Multiplying the observed counting rate by 
reciprocal of this number, 3.49105, gives 
number of disintegrations per second. If it 
desired to express this value in millicuries, t 
observed count is multiplied by 9.44 x10 
assuming the curie as 3.710" disintegrati 
per second. This practice of expressing 
activities in terms of curies is to be condemned, 
will be explained later. 

There are at least three disturbing efle 
common to all Geiger-counter measurements 
beta rays, which must be considered in connect 
with the present arrangement. These are bec 
scattering from the source and support, absorpti 
in the mica window and in some cases 5 
absorption in the source. In the present case 


back-scattering from the support for the sou" 


was reduced to a minimum by the selection 
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exiglas with a considerable cavity behind the 
nurce to reduce this effect. The sources were 
posited on aluminum 0.00035 in. thick. How- 
er, the material of the source itself would cause 


considerable back-scattering, as the phosphorus 
used did not have a high specific activity, and 
there also would be some self-absorption. These 
effects tend to counteract each other. 


III. Independent Check on Solid Angle 


A number of P® sources were measured in the 
rrangement as described. The values obtained 
ere in general agreement with those expected, but 
became obvious that some independent check 
as required to establish the reliability of the 
nyeasurements. For this purpose a thin mica- 
indow counter of the bell type was constructed 

id arranged so that a source could be mounted 
ery close to the window in an accurately reprodu- 
ible location. This was achieved by mounting 
he counter with the mica window up and laying 
he sources, properly centered, directly on the 
punter. The curvature of the mica window 
rovided a small but definite separation between 

and the source. This arrangement is repre- 
ented in figure 3, where S is the source and M the 
ica window. 

To determine the solid angle for this position of 
he source a sample of RaD was prepared from 
ad extracted from pitchblende to be used as a 
tandard. The radium content of this ore was 
arefully measured, and the total amount of lead 
as determined quantitatively. From the above 
lata and the known weight of lead in the standard, 
he number of beta rays per second from the RaE 


n equilibrium with the RaD could be calculated. 








Figure 3.—Mica-window counter. 


unter for Beta Rays 


The ratio of the observed number of counts to this 
calculated value gives the effective aperture of the 
counter. With this information available, an ali- 
quot was taken from one of the P® solutions, of 
such amount to give a convenient counting rate 
in the inverted mica-window counter. The count- 
ing rate multiplied by the above ratio gives the 
number of disintegrations per second in the aliquot 
of P®. Multiplying again by the aliquot ratio, 
the number of disintegrations is obtained in the 
original P* sample, which had previously been 
measured in the chamber with small solid angle. 
In making these measurements in comparison with 
the RaD+E standard, care was taken to deposit 
the P® on the same backing material as was used 
for the standard and over the same area. This 
reduces the back-scattering correction to a small 
value, dependent only on the difference in energies 
of the RaE and P® beta rays. The mica window 
had a thickness ~ 6 mg/cm ?, so that the absorp- 
tion corrections would also be nearly the same. 
When comparing the measurements in the 
arrangement with large solid angle with those 
obtained in the apparatus with small solid angle, 
it is to be noted that whereas the effect of back- 
scattering in the former has been reduced, this is 
not true for the small-angle arrangement as there 
were several milligrams of phosphorus in the 
sources of approximately 0.5-millicurie strength. 
Therefore, we would expect the value obtained in 
this case to be somewhat higher than that obtained 
from comparison with the RaD-+-E standard. On 
the other hand, for the large solid angle arrange- 
ment there is almost no self-absorption of beta 
rays in the material of the source due to the 
smallness of the sample, whereas in the small 
solid-angle apparatus a source of approximately 
0.5 me weighs about 30 mg and extends over a 
little less than a square centimeter. The self- 
absorption will be considerable under these 
conditions. A measurement gave 0.489 me for 
30 mg in the high, large solid-angle arrangement 
and 0.488 mc for the low solid-angle arrangement. 
Under the conditions of these measurements the 
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results can only be explained if the back-scattering 
effect was approximately offset by self-absorption 
within the source. The agreement might be 
regarded as fortuitous. Therefore, the compari- 
son was repeated from the beginning, using freshly 
prepared sources from the same original solution 
of P® in both pieces of apparatus. The measure- 
ments gave 0.504 me in the large solid-angle 
counter and 0.502 mec in the low solid-angle 


arrangement. This makes it quite certain 


the low solid-angle chamber gives a reli 
measurement of beta rays from P® or any soy 
emitting only beta rays of approximately , 


sameenergy. Furthermore, this counting array 
ment is independent of any standard and g 
the disintegration rate simply by multiplying 
observed counting rate by the solid angle 
puted from the dimensions of the apparatus, 


IV. Use of the Curie 


It has become the practice to express the 
strength of all radioactive sources in terms of 
curies. Actually, the curie can be applied logi- 
cally only to members of the radium family. This 
arises from the fact that the curie was originally 
defined as “That quantity of radon in equilibrium 
with 1 of radium.”” Therefore, only members of 
the radium family in equilibrium with 1 g of radium 
can have this disintegration rate. A further 
disadvantage in the use of this unit is that it is 
uncertain to at least 4 percent and values are in 
current use which fall well outside these limits. 

The solution of this difficulty is remarkably 
simple. The intensity of radioactive sources is 
measured actually in terms of disintegrations per 
second. Therefore, all that is required is to 
designate some convenient rate and give it a 
name. We suggest that a convenient rate is 10° 
disintegrations per second, and that the unit be 
called a “rutherford,” abbreviated as “‘rd.””. The 


microrutherford then is the amount undergoing 
disintegration per second and a kilo rutherfy 
10° disintegrations per second. In the futur 
is anticipated that even those nucleii havi 
fairly complex disintegration schemes will 
measurable in terms of disintegrations per secon 
and many can already be treated in this wy 
Thus all confusion that has arisen from the {uy 
application of the curie can be eliminated. 

This subject has been discussed by the Natio 
Research Council Committee on Radioactivity 
which, in turn, has suggested that the Nation 
Bureau of Standards recommend the use of th 
“rutherford” as a general unit for the measu 
ment of radioactive sources. This recommend 
tion was published in Science 103, 712 (1946), ani 
will appear shortly in other journals. 


WasHINGTON, May 23, 1946. 
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feats, Equilibrium Constants, and Free Energies of 
Formation of the Alkylbenzenes' 


William J. Taylor,? Donald D. Wagman, Mary G.Williams,? Kenneth S. Pitzer,’ and Frederick D. Rossini 


For benzene, toluene, ethylbenzene, the three xylenes, normal and isopropylbenzene, 
the three methylethylbenzenes, the three trimethylbenzenes, and the higher normal alkyl- 
benzenes, values are presented for the following thermodynamic properties for the gaseous 
state to 1,500° K: The heat-content function, (H°—H})/T; the free-energy function, 
(F°—H)/T; the entropy, S°; the heat content, H°—H the heat capacity, Cp; the heat of 
formation from the elements, AHf°; the free energy of formation from the elements, AFf°; 
and the logarithm of the equilibrium constant of formation from the elements, logis Kf. 

Equilibrium constants and concentrations are given in tabular and graphical form for 
the isomerization of the four CsHy» alkylbenzenes and for the eight CyH alkylbenzenes as 
a function of the temperature to 1,500° K. Equilibrium constants are also given in tabular 
and graphical form for some reactions involving alkylation (addition of olefin to benzene to 
form alkylbenzene), cyclization (conversion of paraffin to alkylbenzene plus hydrogen), and 
trimerization (of acetylene to form benzene and of methylacetylere to form 1,3,5-trimethyl- 


benzene). 
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I. Introduction 


As part of the work of the American Petroleum 
Institute Research Project 44 and the Thermo- 
chemical Laboratory at the National Bureau of 
Standards, values have been calculated for the 
heat-content function, free-energy function, en- 
tropy, heat content, and heat capacity for the 14 
alkylbenzene hydrocarbons through CyHp», and 
for the higher normal monoalkylbenzenes, in the 
gaseous state to 1,500° K. These data have been 
combined with values of heats of formation at 


25° C to calculate values of the heats, free ene 
ties, and equilibrium constants of formation , 
these compounds in the gaseous state to 1,500°} 
Equilibrium constants have been calculated 4 
some reactions involving alkylation (addition ¢ 
olefin to benzene to form alkylbenzene), cyclin, 
tion (conversion of paraffin to alkylbenzene ply 
hydrogen), and trimerization (of acetylene 4 
form benzene and of methylacetylene to fom 
1,3,5-trimethylbenzene). 


II. Units and Constants 


The unit of energy, atomic weights, and values 
of the fundamental constants used in this report 


are the same as those previously described {|}! 
‘ Figures in prackests indicate the literature references at the end of thispere 


III. Heat-Content Function, Free-Energy Function, Entropy, 
Heat Content, and Heat Capacity 


The thermodynamic functions for benzene, 
toluene, o-, m-, and p-xylenes, and mesitylene 
were published by Pitzer and Scott [9]. The pres- 
ent calculations include minor revisions in these 
values to take account of more recent data, but 
they primarily concern the remaining alkylben- 
zenes through C,H» and the higher normal 
monoalkylbenzenes. The detailed methods, in 
which contributions from translation, over-all 
and internal rotation, and vibration are considered, 
were used for benzene and all its methyl deriva- 
tives and for ethylbenzene. The calculations for 
ethylbenzene and all methylbenzenes with ortho 
methyl groups contain certain additional approxi- 
mations. The functions for n-propylbenzene and 
the higher normal alkylbenzenes, isoproplybenzene, 
and the three methylethylbenzenes, were esti- 
mated by a method of increments from the func- 
tions for the normal paraffin hydrocarbons and 
the lower alkylbenzenes. 


l. Translational and Free Rotational 
Contributions 


For this series of molecules the potential barriers 
to internal rotation are relatively low, and the 
moments of inertia are always fairly large; con- 
sequently, it is convenient to calculate the 


% 


classical free internal rotational contribution ani 
then to correct it in all cases for the potentid 
barrier present. The alternate procedure of ca- 
culating directly the thermodynamic contrib: 
tions for restricted internal rotations would l 
equally correct. 

All the molecules considered in this section fal 
in the class studied by Pitzer and Gwinn {i 
(symmetrical tops attached to a rigid fram 
except ethylbenzene. However, it may be cor 
sidered as a methylene group with a symmetrial 
methyl group and a balanced, but unsymmetrical 
phenyl group attached. By the term “balanced 
it is meant that the center of gravity of th 
phenyl group is on the axis of rotation but the 
the two moments of inertia of the phenyl grow 
perpendicular to the axis of rotation are n0 
equal. An analysis of this type of problem, whic 
shows that in such cases the various moments ¢! 
inertia should be calculated for the equilibrium 
configuration of the molecule, is to be publishe 
elsewhere [4]. Although more general formulas 
are given there, they reduce to those already pul 
lished [5] for this case. 

The formulas of Pitzer and Gwinn [5], with the 
new physical constants mentioned above, com 
bined with the equations for translation and ove" 
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rotation [1], yield the following equations for 
» total contribution of translation, over-all 
ation, and free internal rotation: 


H° —H3)/T=C°,=7.9787 +0.9936m (1) 

(F°—Hs)/T=C, logio T+C, (2) 
C,=18.3026+-2.2878m (3) 
C,=6.8635 logy, M+- 


D878 logic A+2.2878 logo B—4.5757 logio o 
—2.5342m— 10.2960. (4) 


these equations, m is the number of internal 
ations in the molecule; M is the molecular 
ight; A is the product of principal moments 
inertia for over-all rotation, each in units of 
-* ¢ em’; B is the product of reduced moments 
inertia for internal rotation, each in units of 
-” ¢ em®; and o is the total symmetry number, 
e product of the over-all rotation symmetry 
mber and all internal rotation symmetry 
mbers. The reduced moments of inertia for 
ternal rotation were calculated from equations 
b and le of Pitzer and Gwinn, which give a 
ond approximation to the reduced moments 
r molecules with several tops. However, for 
he molecules considered here, except perhaps 
hylbenzene, there is no significant difference 
etween the second approximation and the first 
pproximation given by equation la of these 
thors. 


The values for the bond lengths used in the cal- 
culations were chosen after a consideration of the 
published X-ray and electron-diffraction data. 
The benzene ring was assumed to be a plane 
hexagon, with carbon-carbon bond lengths of 1.39 
angstroms. The other bond lengths, in angstroms, 
were 1.54 for the carbon-carbon bonds in the 
alkyl groups, including the bond joining the group 
to the ring [7], 1.09 for the carbon-hydrogen bonds 
in the alkyl groups, and 1.08 for the carbon- 
hydrogen bonds on the ring. Tetrahedral angles 
were assumed in the alkyl groups. These dimen- 
sions yield for the moment of inertia of the methyl 
group about its axis of rotation 5.3010-*. 
However, the moment of inertia of ethane about 
its threefold axis [8] should be double that of a 
methyl group giving for the latter 5.52x10-*. 
Although no definite decision can be reached at 
this time, the larger value was used in all calcula- 
tions for the methyl derivatives of benzene. The 
calculations for ethylbenzene had already been 
made with the smaller value and have not been 
changed.® 

The values of the constants in eq 1 to 4 are 
summarized in table 1. 


‘ The use of the higher value of the methyl moment of inertia for ethyl- 
benzene would increase the entropy for free rotation by 0.04 cal/deg mole at 
all temperatures. As a result the potential barrier for the rotation of the 
phenyl (or ethyl) group evaluated later in this report from the experimental 
entropy would be increased from 1,080 to 1,130 cal/mole, which is well within 
the limit of other errors. 


TABLE 1.—Constants for the calculation of the translational and free rotational contributions to the 
thermodynamic functions 








Quantity 


i-|1, 2, 4-Tri-j1, 3, 5-Tri- 
methyl- | methyl- 
benzene | benzene 


Ethyl- 


bensene © o-X ylene |m-Xylene 





umber of rotors 


olecular weight... . 











2 2 3 

6 18 27 

106. 160 106. 160 
18. 50 23. 47 
56.71 47.61 
68. 51 70. 02 
71, 920 78, 240 


{ aim. | 2.23 


28. 41 
65. 34 
92. 16 
171, 000 


159. 6 


48.13 
84. 62 
155, 100 


161.7 


39. 07 
64. 32 
66, 140 


| 29. 31 

















—15. 214 | —13. 752 








—4. 763 |] a 
Looe i ay 8. 670 





> Por the configuration in which the carbon atoms of the ethyl group lie 
h a plane perpendicular to the plane of the benzene ring. 

* For the moment of inertia 5.52 for the methyl group; unless otherwise 
Ddicated. 

¢ For the moment of inertia 5.30 for the methy! group; this value was used 
bh the calculation of the thermodynamic functions for ethylbenzene. 
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¢ The principal z-axis is perpendicular to the plane of the benzene ring, 
except for ethylbenzene, for which it makes an angle of 11° with the perpen- 
dicular. The orientation of the principal z- and y-axes may be determined 
by inspection, except for 1,2,4-trimethylbenzene, for which the projection of 
the z-axis on the plane of the benzene ring makes an angle of 13° with a line 
passing through the 1 and 4 positions on the ring. 





2. Vibrational Contributions 


(a) Benzene, Toluene, the Xylenes, and 
Mesitylene 


The vibrational assignments of Pitzer and Scott 
[9] for benzene,* toluene, the three xylenes, and 
mesitylene, were accepted without change. Such 
new data as have become available would not 
appear to make a complete reconsideration of the 
spectra of these compounds at this time worth- 
while. However, the harmonic oscillator contri- 
butions to the thermodynamic functions were 
recomputed at 250° K and at 100-degree intervals 
from 300° to 1,500° K, using the tables of Sherman 
and Ewell [3] and the new values of the physical 
constants [1]. 


(b) Ethylbenzene 


The vibrational frequencies of ethylbenzene 
were assigned on a semiempirical basis that in- 
volved a detailed consideration of only the lowest 
frequencies of vibration of the benzene ring and 
also of the frequencies associated with vibration 
of the ethyl group. The remaining frequencies 
were taken from the assignment for toluene [9]. 
A complete frequency assignment for ethylbenzene 
at this time is impracticable. The details of this 
assignment, together with the method of calcula- 
tion, will be presented elsewhere [6]. The result- 
ing frequencies as used in the thermodynamic val- 
culations are listed in table 2. 


(c) 1,2,3- and 1,2,4-trimethylbenzenes 


The total vibrational contributions to the 
thermodynamic functions of 1,2,3- and 1,2,4- 
trimethylbenzene were estimated empirically from 
the values calculated for 1,3,5-trimethylbenzene 
(mesitylene) and the three xylenes. The differ- 
ences in the vibrational contributions for the three 
xylenes are small, as may be seen from the values 
given in table 3 for the heat-content function, 
free-energy function, and heat capacity, at sev- 
eral temperatures. In view of the small differ- 
erences in the functions obtained with the methyl] 
groups ortho-, meta-, or para- to one another, it 
seems reasonable to conclude that these positional 
effects are approximately additive. This view 
does not take account of the effect of molecular 

* The recent vibrational assignment of Herzfeld, Ingold, and Poole [27] 
yields thermodynamic properties for benzene which do not differ signifi- 


cantly from those presented here. The largest differences are about 0.1 
cal/deg mole in the entropy and heat capacity. 


symmetry on the vibrations. However, 
metry affects the over-all distribution of freqy 


TaBLe 2.—Vibration frequencies of the ethylbenzene 
cule 








Type of vibration « | Symmetry > Freques 
Wie 


emi 





Vibrations principally characteristic of the 
benzene ring or pheny! radical: 
11 (CsHs rocking) - ...._- aginpepltabdedeniiats ; 
18b (C2Hs rocking) 





>—H stretching. _. 
Da 
20b 
Vibrations principally characteristic of the 
ethyl group: 
CsHs—C Hi—CH,; angle bending 
CH; rocking a 
CH,—CH,; stretching 
CH; rocking 
CH, rocking... ....._. 


CH; twisting__- 

CH; rocking 

CH; symmetrical bending 

CH, symmetrical bending __. 

CHs unsymmetrical bending An 

CH; symmetrical stretching. ............. , 

CH) unsymmetrical stretching 

CHa: symmetrical stretching A’ 

CHs unsymmetrical stretching S. 

Internal rotations: 
CsHs—C2Hs rotation A” 
CHs rotation ; Aa 











* The numbers refer to the normal modes of vibration for the benz 
molecule, as given in figure 6 of reference [9]. 

> The ethylbenzene molecule, in the equilibrium configuration assume. 
has the symmetry Cis, with the symmetry plane perpendicular to the plas 
of the benzene ring. Vibrations of the type A’ and A” are symmetric so! 
antisymmetric, respectively, to reflection in the symmetry plane. 
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»s, and therefore the vibrational function, to a 
uch smaller extent than it does the individual 
quencies. It is easily verified that, on the 
is of the assumption that positional effects are 
jditive, the proper formulas are 


2,3-trimethyl benzene) = (1,3,5-trimethylben- 
zene) + 2{(o-xylene) — (m-xylene)] 


2 4-trimethyl benzene) = (1,3,5-trimethylben- 
ene) +[(o-xylene) + (p-xylene) — 2(m-xylene)} (6) 


(5) 


The vibrational functions for 1,2,3- and 1,2,4- 
imethylbenzene were calculated by means of 
ese formulas at each temperature. The result- 
ng functions are compared with those for 1,3,5- 
imethylbenzene at several temperatures in 
able 4. 


. Restricted Rotational Contributions 


The calculation of the translational, total free 
over-all and internal) rotational, and vibrational 
ntributions to the thermodynamic functions 
ave been described in the preceding sections for 
nzene, toluene, ethylbenzene, the three xylenes, 
nd the three trimethylbenzenes. This com- 
letes the calculation for benzene and leaves only 
he restricted internal rotational contributions for 
he other molecules. The potential barriers 


restricting internal rotation of the methyl groups 
in toluene and the three xylenes have been de- 
termined by comparison with the experimental 
entropies and heat capacities. The potential 
barrier for the methyl rotation in ethylbenzene 
was assigned a reasonable value empirically, and 
the barrier for the ethyl rotation was then obtained 
from the experimental entropy. The barriers for 
the methyl rotations in the trimethylbenzenes 
were assigned on the basis of the barriers in the 
xylenes. The details of these calculations will be 
presented in this section. 


(a) Symmetry of Potential Barriers 


The usual assumption has been made that the 
potential barrier for each internal rotation may be 
approximated by an n-fold cosine function and that 
interaction terms in the potential energy between 
different internal rotations may be neglected. 
The barrier for the methyl rotation in toluene is 
sixfold (n=6), and the barriers for the methyl 
groups in p- and m-xylenes must have essentially 
the same magnitude and symmetry as the barrier 
in toluene. The situation in o-xylene has been 
discussed by Pitzer and Scott [9], and their 
recommendation that the methyl groups be treated 
as independent rotors, with n=3, has been fol- 


TaBLe 3.—Vibrational contributions for the zylenes 











Vibrational contribution in cal/deg mole 





(H°— H/T 


(F°— H/T Cc; 





o Xylene m-Xylene o-Xylene 


m-Xylene | p-Xylene 


m-Xylene 


o-Xylene 








7.92 7. 88 


31.84 








—4.13 
—10.17 
—27. 69 
—44. 25 


—4.33 
—10. 39 
—27.80 
—44, 29 


20. 30 
38. 30 


75. 93 

















TABLE 4.—Vibrational contributions for the trimethylbenzenes 











Vibrational contribution in cal/deg mole 





(H°—Hi)/T | 


(F°—Hi)/T Cc; 





1,3,5- 
| Trime thyl- | 
benzene 


1,2 
Trime thyt- | 
benzene 


1,2,3- 
Trimethyl- | 
benzene 


1,2,3- 
ney Trimethyl- 


1,24 | 1,2,4- 
Trimethyl- 


benzene 


1,3,5- 
Trimethyl- 
benzene 


1,2,3- 
Trimethyl- 
benzene 


1,3,5- 
Trimethyl- 


benzene benzene 








Yr 
benzene 
| 








| 
| 
| 
| 
| 


| —5.79 
—13.06 
—33. 80 
—53. 25 
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lowed. The barrier for the ethyl (or phenyl) 
group rotation in ethylbenzene is twofold (n=2). 
Although the barrier in this case may deviate 


considerably from the simple cosine form, suffi- 
cient calorimetric data are not available to justify 
a more detailed treatment. 


(b) Reduced Moments of Inertia 


Although only the product of reduced moments 
of inertia was needed for the free-rotation cal- 
culation, individual values are required for re- 
stricted rotation calculations [5], since the parti- 
tion function for each individual rotation is 
needed. In terms of our constants it is 


Q,=2.7935 (10°J7')*/n, (7) 


where J is the reduced moment of inertia, 7 the 
absolute temperature, and n the symmetry num- 
ber. The values of J were calculated from eq 
1b and Ic of reference [5], and are listed in table 5. 











Taste 5.—Reduced moments of inertia for internal 
rotation * * 


[Values in units of 10- g cm‘) 





edehcs | e518 | | 

as Ca dee ae 5.41 . be aes 
BD TE eC 

RM 

i danclnbibteis tia’ 5.44 | 








* The values refer to the groups in the order mentioned in the names for the 
trimethylbenzenes. The methyl group is given first for ethylbenzene. 
These values are rounded to three significant figures; hence their product 
will not agree exactly with that in table 1. 

> For the moment of inertia 5.52 for the methyl group, unless otherwise 
indicated. 

¢ For the moment of inertia 5.30 for the methyl group. 


(c) Calculation of Potential Barriers from 
Experimental Entropies 


The experimental entropies for these com- 
pounds, for the liquids under their saturation 


pressures, at 298.16° K, are given in the first row 


of table 6. These values have been taken as the 
standard entropy of the liquid, at 1-atmosphere 
pressure and 298.16° K, as the difference is only 
0.002 or 0.003 cal/deg mole. The references to 
the sources of these values are given in the table. 
The second row gives the standard entropy of 
vaporization, ASv°, at 298.16° K, as taken from 
the work of Osborne and Ginnings [10], and 
reduced to the standard state [14]. The third 
row (sum of first and second) is the experimental 


100 


value for the standard entropy of the gas, yj 
atmosphere pressure and 298.16° K. The {oy 
row gives the sum of the calculated contribyi 
to the standard entropy of the gas for translajj 
total free (over-all and internal) rotation, 
vibration, at 298.16° K. The fifth row @ S° (lia 
ference of third and fourth) is the decrease jy ¢ 
entropy caused by the restriction of the inte, 
rotations (except for benzene). Finally, the 
three rows give respectively the maximum ,e""" 
tential barrier consistent with the experime 
entropy and its uncertainty, the probable bam aie 
indicated by the experimental entropy value iti! “o™ 
and the minimum barrier consistent with 4 
experimental entropy and its uncertainty. 

All contributions of restricted rotation thro 
out this paper are taken from tables of Pitzer » 
Gwinn [5]. 

For ethylbenzene there are two separate pot 
tial barriers to be evaluated. The barrier { 
methyl rotation has been assumed to have the ¢ 
value of 3,400 cal/mole found in propane [1MMB). offe 
On this basis, 1.50 cal/deg mole of the entropy Scott 
ference, (S;—S°), is to be assigned to the metiyimms 
rotation, leaving 0.36 cal/deg mole as the cont 
bution of the ethyl rotation, which leads to 1, 
cal/mole as the indicated value of the potent 
barrier for the ethyl (phenyl) rotation. In th 
case, the uncertainty in the calculated entropy: 
vibration and of restricted methyl rotation ¢ 
ceeds that of the experimental value. Hence, 
culated maximum and minimum potential bw 
riers would have no meaning. 

The experimental value for the entropy : 
benzene gas at 298.16° K differs from the calculat 
value by 0.19 +0.30 cal/deg mole. Even ti 
insignificant difference is in the direction to} 
explained by the neglected anharmonicity : 
vibrations. 


ptential b 


* Referer 
> Referen 
¢ For the 
400 cal/m: 


(d) Calculation of Potential Barriers from 
Experimental Heat Capacities 


Xylene ' 


The comparison with the gaseous heat-capaci 
data of Gwinn, and of Pitzer and Scott [9], is sho 
in table 7. The experimental values are given! 
the third column, and the values calculated for th 
sum of translation, total free rotation (over-all a 
internal), and vibration are given in the four 
column. The difference of these values, given! 
the fifth column, is the increase in heat capac! 
owing to restriction of rotation (except for be 
zene), The values of the potential barriers int 
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TaBLe 6.—Polential barriers from entropy values at 298.16° K 
(Entropy, cal/degree mole; potential barrier, cal/mole]] 































| | 
Benzene Toluene Ethylbenzene oXylene | m-Xylene p-Xylene 
| | 
————————— | | 
perimental entropy: | 
a Ghee 1 CPs ~ 5.555 54 - 1h cc ctay cb cqee jn wteessssaits | 41. 49 52. 39 60. 99 58. 80 60. 42 59. 20 
Vapertantlem, AMP uc... .-----e-nereecenvenecncesreconereay | 23. 04 23. 94 25. 16 25. 40 25. 22 25. 11 
pray Ee a 8 | 64. 53+ 0.30 | 76.33 40.25 | 86.15 40.12 | $4.20 40.25 | 85.64 40.35 | 84.31 40.25 
culated entropy of translation, vibration, and over-all and free in- 
ernal rotation, Sy (gas, 1 atm)_..........-.....-....----- 64. 34 76. 56 88.01 86.15 85. 80 #4. 51 
erence (S; —S°) 0. 23 1. 86 1.95 0. 16 0.20 
ptential barrier: 
pheoteees ©. ... .cncccenceccdpoonencessbseneecngnes Jo----nncnnee 1, 600 (*) 2, 5300 1, 000 900 
wtidid aE cs coieeckdaaks+<0snseeakieleeuns SAE At SE ee 1, 000 © 1, 080 2, 200 550 600 
Detniemams 6... 22sec ce nnddn nc eennn-ovccmccccecoceseses eecccceesenes pesorccece 0 (*) 2, 000 0 0 




















+ References: Benzene [17], toluene [18], ethylbenzene [19, 20], xylenes [9]. 
> Reference [14]. 

¢ For the ethy! (phenyl) rotation in ethylbenzene assuming a barrier of 
400 cal/mole for the methy] rotation. 









ated by the data are shown in the last column. 
hese values are obviously subject to considerable 
neertainty, in view of the estimated uncertainties 
the experimental heat capacities, and the possi- 
ble effects of error in vibrational assignment. 

Scott and Mellors [11] have reported the value 
8.06 +0.20 cal/deg mole for the gaseous heat 
apacity (C,) of ethylbenzene at 373.16° K. 
This value was obtained by applying a calibration 

















TasLe 7.—Potential barriers from gaseous heat-capacity 
data 











| | Heat capacity C°, | Difference 
































- ; | Caleulated 
—— = — (expere | pote atial 
pound 7 | | | mental | barrier 
Experi- Calculated minus (Vv) 
mental * | (V=0) | calculated) | 
i, i. | eid nee ; 
°K | cal/deg mole\ cal/deg mole! cai/mole 
388 | 26.0-40.3 | 25.96 | O08 LL. 
Benzene : 417 | 28. 10.3 27.85 | x 4 
481 | 31.4+0.3 31.73 —. 33 
| | 
393 | 26.440.3 26. 29 | Hu 
enzent 428 | 28. 440.3 28. 55 | —.15 |. 
463 30. 340. 3 | 30. 70 | | re a 
| {393 32. 840.3 | 32. 57 | 23 1, 000 
Toluene ¢ 428 | 35.740.4 35. 27 .43 1, 500 
463 | 38.040. 4 | 37. 87 | .13 750 
303 | 40.240. 4 38. 62 1.58 1, 900 
Xylene 428 | 43.5-40.4 | 41.72 1.78 2, 200 
463 46. 040.5 | 44.72 | 1. 28 1, 800 
| [393 | 39. 140.4 | 39. 12 —0.02 0 
m-Xylene > 428 | 42.440.4 | 42.25 15 600 
| (463 45.440. 4 | 45. 24 .16 600 
| | 
| 393 | 39.0+0.4 | 38. 78 22 650 
p-Xylene >... ___ | (428 | 42.640.4 41.90 | .70 1, 300 
| [463 | 45.2+0.4 44. 89 31 850 








* The first series for benzene is the data of Gwinn [9]; the remaining data are 
those of Pitzer and Scott [9]. 

‘Toluene, m-xylene, and p-xylene calculated for a potential, barrier of 6 
maxima; o-xylene with a barrier of 3 maxima. 
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4 Allowed by the experimental uncertainties. 
* Uncertain—see text. 





correction to the observed value of 37.71. The 
calculated value at this temperature for transla- 
tion, total free rotation, and vibration is 36.56, 
leaving 1.50 cal/deg mole as the contribution of 
restricted rotation. Of this, 1.20 is to be assigned 
to the methyl rotation, assuming, as before, a 
3,400 cal/mole barrier. The contribution of the 
ethyl rotation is therefore 0.30 cal/deg mole, which 
yields a value of 850 cal/mole for the barrier. 


(e) Selected Values of Potential Barriers 


The potential barriers selected on the basis of 
the entropy and heat-capacity data, and used in 
the calculation of the final values of the thermo- 
dynamic functions, may now be stated. The 
barriers for the methyl groups in toluene, p- and 
m-xylenes were assigned a value of 750 cal/mole, 
as the best average value indicated by the data 
on these compounds,’ and a value of 2,100 
cal/mole was assigned to the barriers in o-xylene. 
The differences between these values and the 
values reported previously [9] are attributable 
principally to changes in the calculated vibra- 
tional contributions resulting from the use of the 
new values of the fundamental physical constants 
[1]. The barrier for the ethyl group rotation in 
ethylbenzene was assigned the value 1,080 cal/- 
mole indicated by the entropy. The entropies 
calculated with these barriers are compared with 
the experimental entropies at 298.16° K in table 8. 





’ This value of 750 cal/mole is to be regarded as an empirical parameter for 
best reproduction of thermodynamic functions. The low-temperature heat- 
capacity curves of toluene and m- and p-xylene make it improbable that 
the real value is this large—see Pitzer and Scott [9], page 826. The thermo- 
dynamic data are in reasonable accord with the 500 cal/mole value adopted 
previously [9}. 
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TABLE 8.—Comparizon of experimental and calculated 
entropies at 298.16 ° K 





Entropy, S° (gas, 1 atm, 
298.16° K) 
Compound 





Experimental | Calculated 





cal/deg mole 
64. 533+0. 30 
76. 3320. 25 76. 42 
86. 154-0. 12 86.15 
84. 20-0. 25 84. 31 
85. 6440. 35 85. 49 
84. 3140. 25 4.23 


cal/deg mole 
4.44 











The barriers for the trimethylbenzenes were 
assigned empirically on the basis of the barriers 
in the xylenes, as follows: (1) 1,2,3-, and 1,2,4- 
trimethylbenzenes, one 750, and two 2,100 
cal/mole barriers, each; (2) 1,3,5-trimethylben- 
zene, three 750 cal/mole barriers. The deriva- 
tion of these values is obvious except in the case 
of the 750 cal/mole barrier assigned to the central! 
methyl group in 1,2,3-trimethylbenzene. It is 
likely that the interactions of the central methyl 
group with the adjacent groups on either side 
increase about equally the maximum and mini- 
mum of the potential energy for the rotation, 
and do not, therefore, alter greatly the potential 
barrier (variation of potential energy) for the 
rotation. 


4. Calculations by the Method of 


Increments 


The thermodynamic functions for the C,H, 
alkyl benzenes, aside from the three trimethyl- 
benzenes, were estimated from the functions for 
appropriate lower alkylbenzene and paraffin hy- 
drocarbons. The formula used for the three 
methylethylbenzenes was 
1-methyl-z-ethylbenzene= 1 ,z-dimethylbenzene 

+ethylbenzene—toluene+R In 2, (8) 
where z=2,3, and 4 For propylbenzenes, 
n-propyl benzene = ethylbenzene + n-butane 

—propane, (9) 
isopropyl benzene = ethylbenzene + isobutane 
—propane+ R In (3/2). (10) 


The higher normal alkylbenzenes, from n-butyl- 
benzene (C\pH,) to n-hexadecylbenzene (C.,H;.), 
were calculated from the formula 


normal alkylbenzene (C,,H;,_.) ethylbenzene 
(C,H,)+normal paraffin (C,_;H,,_s) 
—propane (C;H,), (n=9). (11) 


When n =9 this reduces to the formula for n- 
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propylbenzene (CjxHy») eq 9. The thermody 
mic functions for isobutane were taken {y 
reference [23], and for the normal paraffin hy 
carbons from reference [16]. The logarith 
terms correct for the discrepancies in total gy 
metry numbers for over-all and internal rota; 
in the formulas as written. They are t | 
included in the calculation of the entropy and 
negative of the free-energy function, —(P 
H,)/T, but omitted in the calculation of the ly 
content function and the heat capacity. 

Some correction may be desirable for ¢ 
changes in the barriers for rotation of the 
groups attached to the benzene ring, especial 
in 1-methyl-2-ethylbenzene and isopropylbenze 
but it was not considered worthwhile to estim 
the corrections in view of the uncertainty , 
values calculated by the present method. 

An indication of the accuracy of the values 
n-propylbenzene may be obtained by compariy 
the directly calculated values for ethylbenz 
with values calculated from the equation 


ethylbenzene=toluene+ propane— ethane. (|! 


The differences, expressed as the value calculate 
from eq 12, less the directly calculated value,» 
cal/deg mole, at 300°, 900°, and 1,500° K, » 
1.17, 0.10, and —0.17, for the free-energy fun 
tion; —1.24, —0.67, and —0.42 for the hew 
content function; and —0.92, —0.13, and —( 
for the heat capacity. However, in this case tle 
difference in the contribution of the methyl rw 
tion in toluene and the ethyl] rotation in ethylbe- 
zene is considerable. Correction for the chang 
in this contribution leads to differences of (4 
—0.32, and —0.42, for the free-energy function; 
—0.29, —0.28, and —0.18, for the heat-conte! 
function; and —0.59, —0.09, and —0.01, for th 
heat capacity. A correction of this type shoul 
be unimportant in eq 9 for n-propylbenzene. 


5. Tables of values of the thermod;- 
namic properties 

The resulting values of thermodynamic pm 
erties for the 14 alkylbenzenes through C,H, a 
for the higher normal monoalkylbenzenes are pi 
sented in tables 9 to 18, which give values of te 
heat-content function, (H°—H)/7, the {re 
energy function, (F°—Hg)/T, the entropy, © 
the heat content, H°—H%, and the heat capaci) 
C3, for the gas in the standard state at a pressut 
of 1 atmosphere. 
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TaBLE 9.—Values * of the heat-content function, (H°—H),)/ T, for the 14 alkylbenzenes, C, to Cy, for the ideal 
gaseous state, to 1,500° K 





Temperature * in ° K 








cane f=[= [= [= [=| [=| = [vem sae]iae [ve 


Heat-content function,» (H°—H})/T, in cal/deg mole 








CoH 17. 50 | 20.48 | 23.24 | 25.76 | 28.07 | 30.16 | 32.07 | 33.82 
Cris 21.94 | 25.56 | 28.92 | 32.03 | 34.86 | 37.45 | 39.82 | 41.99 


CsHw , 26.99 | 31.33 | 35.35 | 39.08 | 42.38 | 45.45 | 48.25 | 50.82 
Dimethylbenzene (o-xylene).....-| CsHe ‘ 27.61 | 31.78 | 35.66 | 39.23 | 42.52 | 45.53 | 48.30 | 50.84 
Dimethylbenzene (m-xylene)......| CsHe 26.72 | 30.94 | 34.88 | 38.52 | 41.86 | 44.93 | 47.74 | 50.32 
Dimethylbenzene (p-xylene) CsHw 26.66 | 30.83 | 34.74 | 38.36 | 41.69 | 44.76 | 47.56 


CoHu . ° 32.2 | 37.2 | 41.8 | 46.0 | 49.8 | 53.4 | 56.6 | 50.5 
CosHn . . 31.4 | 36.6 | 41.3 | 45.6 | 49.5 | 53.1 | 564 | 50.4 
CHa . 32.7 | 37.6 | 421 | 46.2 | 530.0 | 53.5 | 56.7 | 50.7 
C;.Ha . . 31.8 | 36.7 | 41.3 | 45.5 | 49.4 | 529 | 562 | 50.2 
CsHu . 31.7 | 36.6 | 41.2 | 45.4 | 49.2 | 528 | 560 | 50.0 


Csi 32.0 | 36.6 | 41.1 | 45.2 | 48.9 | 524 | 55.6 | 58.6 





CoH , 32.1 | 36.8 45.3 | 40.1 | 52.6 58.7 
CoH 36.04 | 40.55 | 44.74 | 48.50 | 52.13 . 58. 39 























8. 5Trimethylbenzene (mesitylene) _ 





























* Interpolation to other temperatures in the interval 298.16° to 1,500° K may be made by appropriate graphical or analytical methods. 

ween 200° and 298.16° K, values may be estimated by extrapolating the values for 300°, 400°, 500°, and 600° K. The values in this table are given to more 
ificant figures than are warranted by the absolute accuracy of the individual values in order to retain the internal consistency of the several thermodynamic 
ctions of a single substance, and also to retain the significance of the increment with temperature of a given thermodynamic function. 

> The heat-content function (H°—Ho)/T, is the heat content at the given temperature less the heat content at 0°K, divided by the absolute temperature 
K) of the given hydrocarbon in the thermodynamic standard gaseous state of unit fugacity (1 atmosphere). 


ABLE 10.—Values* of the heat-content function, (H°— Hj)/T, for the normal alkylbenzenes, Cs to Cx, for the ideal gaseous 
state, to 1,500° K 








Temperature * in ° K 





Compound (gas) 0 | 208.1 300 | 400 500 | 600 | 700 | 800 900 | 00 | 1,300 |, 20 | 2,00 | 1,400 | 





Heat-content function,» (#7°—Hj)/T, in cal/deg mole 





17. so] 20.48) 23.24) 25.76; 28.07) 30.16, 32.07) 33. 
21.94) 25.56) 28. 32.03) 34.86) 37.45) 39.82) 41.99 
26.99) 31.33) 35.35) 39.03) 42.38) 45.45) 48.25) 50.82 
32.2/ 37.2) 41.8) 460] 49.8) 53.4) 566/ 59.5 
37.1) 428) 48.0/ 528] 57.1] 61.1) 64.8) 68.2 
42.0) 483; 42) 3.5) 64.4) 689) 73.0) 76.8 
46.9; 53.9) 60.4) 66.3) 71.7) 767) 81.2) 85.4 
51.8) 50.5| 66.6) 73.1/) 79.0) 841.4) 89.4) 940 
56.8) 65.1) 728) 79.9/ 863); 922) 97.6) 1026 
61.7 | 70.7) 79.0) 86.6; 93.6 | 100.0 | 105.8 | 111.2 


Cols 
Cris 
CsHw 
CsHw 
-| Cro 
.| CuHwe 
Cyulin 
CuH» 
CuHa 
Cun 


escocooooosososo 


66.6) 76.3) 85.3 | 93.4 | 100.9 | 107.8 | 114.1 | 119.9 
81.9 91.5 | 100.2 | 108.2 | 115.5 | 1223 | 128.5 
76.4 | 87.4) 97.7 | 107.0 | 115.4 | 123.3 | 130.5 | 137.1 
$1.3 | 93.0 | 103.9 | 113.8 | 122.8 | 13].1 | 138.7 
86.2 | 98.6 | 110.1 | 120.5 | 130.0 | 138.8 | 146.9 
91.2 | 104.2 | 116.3 | 127.3 | 137.3 | 146.6 | 155.1 | 163.0 
96.1 122.5 | 134.1 | 144.6 | 154.4 | 163.3 


Cues 
Cus 
CisHe 








eoocoooceoso 












































Ss 


6.21) 6.78) 7.29, 7.77) 821 
| 





* See footnote“ a” of table 9. 
_* The heat-content function, (H°—Hg)/T, is the heat content at the given temperature less the heat content at 0°K, divided by the absolute temperature 


K) of the given hydrocarbon in the thermodynamic standard gaseous state of unit fugacity (1 atmosphere). 
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Tasie 11.—Values * of the free energy function, (F°— H)/T, for the 14 alkylbenzenes, Cy to Cy, for the ideal gaseous sic, 
1,500° K 








Temperature * in ° K 





Peels [= [=e] = |= | [|e fifa 





Free-energy function, (F°—H)/T. in cal/deg mole 








CoHs —53. 00 —60. 24) —63. 70) —67.06| —70.34| —73. 50) —76. 57 leat —82. 40) —85. 18 
Methylbenzene (toluene).....| CrHs —62. 07) —T71. 20|—75. 52) —79.72| —83.79| —87. 72) —91. 53) —95. 21) —98. 77)—102. 21 


CsHw — 68. 37 —79. 64|—84. 94) —90. 08) —95.05) —99. 84) —104. 47/— 108. 04) —113. 25) —117. 42) —121. 44-1 
1,2-Dimethylbenzene (o0-xyl- 

ene). CiHw —65. 73 —77. 40| —82. 81] —88. 01) —93. 01) —97. 82! —102. 46} —106. 93) —111. 24) —115. 39) —119. 42-1» 
1,3-Dimethylbenzene (m-xy]- 

ene)... CsHw —67.74 —78. 95} —84. 20) —89. 28) —94. 18) —98. 91/—103. 48) —107. 89) —112. 15] —116. 27| — 120. 25 -1y 
1,4-Dimethylbenzene (p-xyl- 
CsHe —66. 37 —77. 58) —82. 83) —87. 89, —92. 76) —97. 48|— 102. 02) —106. 42) — 110. 66| —114. 77|—118, 74/-19, 


n-Propylbenzene.............| Celis —74. 19 —87.8 |—94. 1 |—100.1 |—106.0 |—111.6 |—117.1 |—122.3 |—127.4 |—132.2 |—136.9 |; 
Isopropylbenzene (cumene) Cow —72. 4 . 2 |—85.6 |—91.8 | —97.8 |—103.6 |—109.2 |—114.6 |—119.8 |—124.8 |—129.6 
1-Methy!-2-ethylbenzene CHa —73. 41 —87 2 |—93.6 | —99.8 |—105.6 |—111.3 |—116.8 |—122. 0 |—127.1 |—132.0 
1-Methyl-3-ethylbenzene..___.| Cellu —75. 42 —88.8 |—95.0 |—101.0 |—106.8 |—112.4 |—117.8 |—123.0 |—-128.0 —132.9 | 
1-Methyl-4-ethy! benzene. - CsHe —74. 05 —87.4 |—93.6 | —99.6 |—105.4 |—111.0 |—116.3 |—121.5 |—126.5 |—131.4 
1,2,3-Trimethylbenzene 
(hemimellitene) Css —T1 53 —&.9 |—91.2 | —97.2 |—102. 9 |—108.5 |—113.8 |—119.0 |—123.9 |—128.7 
1,2,4-Trimethylbenzene 
(pseudocumene)............| CsHn ~—72. 70 —86.2 |—92.4 | —98.4 |—104.2 |—100.8 |—115.1 |—120.3 |—125.3 —130.1 
1,3,5-Trimethylbenzene | } 
(mesitylene)................| CoH —71. 06 Rene Paes —96. 08| — 101. 77, — 107. 26) —112. 56) —117. 68) —122. 62 —127. 40 —132. 03 -1m 


| | | | 


















































* See footnote “a” of table 9. . 
> The free-energy function, (F°—H,)/T, is the free energy (exclusive of nuclear spin) at the given temperature less the heat content at 0° K, divided by 
absolute temperature (° K), of the given hydrocarbon in the thermodynamic standard gaseous state of unit fugacity (1 atmosphere). 


Tasie 12.—Values * of the free-energy function, (F°— H§)/T, for the normal alkylbenzenes, Cy to C22, for the ideal gas 
state, to 1,500° K 








Temperature * in ° K 





Compound (gas) < o| 208.16 | 300 | wo | 00 | oo | 700 | 800 900 | 1,000 | 1,100 | 1,200 | 1,300 | 1.400 | 





Free-energy function >, (F°—H?)/T, in cal/deg mole 


| | | | 
Calls 0; —52.93 | —53.00 | —56. 69) —60. 24) —63. 70) —67. 06, —70. a —73. —76. 57| —79. 54) —82. 40) —85. 18) —87. 85 -a6 





Methylbenzene (tolu- | | 
ene)... sewee C:Hs 1 —61.98 | —62.07 | —66.74) —71. 20) —75, 52} —79. 72, —83. 79) —87. 72) —91. 53] —95, 21) —98, 77; —102, 21) —105, 53/-1K) 


J 
Ethylbenzene.__. CsHw | —68. 26 | —68.37 | —74. 14) —79. 64) —84. 94) —90. 08) —95. 05) —99, 84)— 104. 47|— 108, 94) — 113. 25) —117. 42, —121. 44-18 
n-Propylbenzene ..--| CoH | —74.05 | —74.19 | —81.2 | —87.8 | —94.1 |—100.1 |—106.0 |—111.6 |—117.1 |—122.3 |—127.4 |—132. 2 |—136.9 -i 
n-Butylbenzene. CwHu —79. 79 | —79.94 | —88.0 | —95.6 |—102.8 |—109.8 |—116.6 |—123.0 |—129.3 |—135.3 |—141.1 |—146.6 |—152.0 |-1%! 


n-Amylbenzene. _ . -| CuHiw —85. 54 | —85.72 | —04.9 |—103.5 |—111.7 |—119.6 |—127. 2 |—134.5 |—141.6 |—148.3 |—154.8 |—161.1 |—167.1 -19 
n-Hexylbenzene_ CuHu | —91.30 | —91.49 |—101.7 |—111.4 |—120.5 |—129.4 |—137.8 |—146.0 |—153.8 |—161.3 |—168.6 |—175.6 |—182.3 -I8! 
n-Heptylbenzene Cue 0, —97.05 | —97. 27 |—108.6 |—119.2 |—129.4 |—139.1 |—148. 4 |—157.4 |—166.1 |—174.4 |—182.3 |—190.0 |—197.4 -2K/ 
n-Octylbenzene -| CHa | 0)—102. 80 |—103.05 |—116.5 |—127. 1 |—138.2 |—148.9 |—159. 1 |—168.9 |—178.3 |—187.4 |—196. 1 |—204.5 |—212.5 -m! 
n-Nonylbenzene ..| CuHe 0|—108. 56 |—108. 82 |—122. 3 |—135.0 |—147.1 |—158. 6 |—169.7 |—180. 4 |—190.6 |—200. 4 |—209.9 |—218.9 |—227.7 -2 


n-Decylbenzene. _ . Cul duis —114. 60 |—129, 2 |—142.9 |—155. 9 |—168. 4 |—190.3 |—191. 8 |—202.8 |—213. 4 |--223.6 |—233.4 |—242.8 -21! 
n-Undecylbenzene CuHis 0| —120. 06 —120. 38 |—136. 1 |—150.8 |—164.8 |—178. 1 |—191.0 |—203.3 |—215.1 |—226. 4 |—237. 4 |—247.9 |—257.9 -2! 
n-Dodecylbenzene.......| CisHse | 0) —125. 81 |—126. 16 |—142. 9 |—158. 7 |—173. 6 |—187.9 |—201.6 |—214. 7 |—227.4 |—239. 5 |—251. 1 |—262.3 |—273.1 -2i! 
n-Tridecylbenzene CyHlas | 0|—131. 57 |—131. 93 |—149.8 |—166.6 |—182. 4 |—197.6 |—212. 2 |—226. 2 |—239. 6 |—252. 5 |—264.9 |—276.8 |—288.2 -™ 
n-Tetradecylbenzene Coola | 0)—137. 32 |—137. 71 |—156.7 |—174. 5 |—191.3 |—207.4 |—222. 8 |—237.6 |—251.9 |—265. 5 |—278.6 |—201.3 |—303.3 -vi4! 

















n-Pentadecylbenzene CuHgs | 0|—143. 07 |—143. 49 |—163.5 |—182.4 |—200. 1 |—217. 1 |—233. 5 |—249. 1 | -264.1 |—278.5 |—292.4 |—305.7 j—318 5 |-i! 
nHexadecylbenzene Calis | 0)/—148. 83 |—149. 26 |—170.4 —190. 2 |—209. 0 —226.9 |—244.1 taoe hanes —291.6 —306. 2 |—320.2 |— 333.6 —s6! 


| 
—7.89, —8.84) —9. 75) aq —11. 46, —12. 26, —13.02 —13. 76, —14. 46 —15.13 =i 




















Increment 








| 


| 





0) —5. 753 —5.777| —6.87 
| 


' 





® See footnote “a” of table 9. 
» The free-energy function (#°—Hg)/T, is the free energy (exclusive of nuclear spin) at the given temperature less the heat content at 0° K, divide: " 
the absolute temperature (°K), of the given hydrocarbon in the thermodynamic standard gaseous state of unit fugacity (1 atmosphere). 
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Tape 13.—Values * of the entropy, S°, for the 14 alkylbenzenes, C, to Cy, for the ideal gaseous state, to 1,500° K 





Temperature * in ° K 





Compound (gas) o| 8.16] 300 | 400 | 0 | oo | roo | sao | wo | 3,00 | 





Entropy,® S°, in cal/deg mole 





84.17) 90.30] 96.10) 101.57) 106.73 
101. 08} 108. 64) 115.81) 122. 58) 128. 98 


116. 28) 125. 43) 134.08) 142.23) 149. 92 
114. 60} 123. 67) 132. 24) 140. 34) 147. 99 
D 115. 14} 124. 16} 132.70) 140. 78) 148. 41 
Dimethylbenzene (p-xylene) - - - -- - . 113. 66) 122. 63) 131.13) 139.17) 146. 78 


ypylbenzene : 131.2 | 141.9 | 152.0 | 161.5 | 170.5 
propylbenzene (cumene) - ' . 128.3 | 139.0 | 149.2 | 158.7 | 167.6 
ethy]-2-ethylbenzene - - ro 119.9 | 131.2 | 141.8 | 151.9 | 161.4 | 170.3 
ethyl-3-ethyl benzene - - - Sisce L hb . 120.5 | 131.7 | 142.3 | 152.3 | 161.8 | 170.7 
ethyl4-ethylbenzene 119.1 | 130.2 | 140.8 | 150.8 | 160.2 | 169.1 
s-Trimethylbenzene (hemimelli- 





116.9 | 127.8 | 138.3 | 148.1 166. 2 





Oo) ee . ‘ 118.2 | 129.2 | 139.6 | 149.5 167.7 
hylbenzene (mesitylene) - 1 126. 22) 136. 64) 146. 51 164. 69 























5-Trimet 


























* See footnote “a” of table 9. 
» S* is the entropy (exclusive of nuclear spin) of the given hydrocarbon in the thermodynamic standard gaseous state of unit fugacity (1 atmosphere) at the 


pperature indicated. 


Tapie 14.—Values * of the entropy, S°, for the normal alkylbenzenes, Cz to Cr, for the ideal gaseous 











Temperature * in ° K 


Compound (gas) | 1,000 | 1,100 | 1,200 | 1,300 1, 400 | 1, 500 





Entropy,® S°, in cal/deg mole 


j | | | 
wa--ee-] CoHle 0 04. 34 | 64.46 77.74 84. i 90. 30) 106.73) 111.61] 116.22! 120.59] 124.74) 128. 68 
thylbenzene (toluene) - | CrHs | 0! 76.42 | 76.57 | 84.91) 93. 13) 101.08) 108. 64) 22. 58) 128.98) 135.03] 140.76) 146. 19] 151. 34) 156. 25 
bylbenzene : CsHw| 0) 86.15 | 86.34 . 59] 106.63) 116.28) 125. 43) 1: 2 23} 149.92) 157. 19} 164. 07) 170. 50] 176. 79) 182. 67 
CoH | 0) 95.74 | 95.97 | 108.1 | 119.9 | 131.2 | 141.9 | 152. .5 | 170.5 | 178.9 | 186.9 | 194.5 | 201.7 | 208.6 


CusHu | 0104.91 |105. 16 | 119.2 | 132.7 145.6 | 157.8 | 169. .2 | 190.4 | 200.1 | 200.2 217.9 | 226.1 | 234.0 


CuHe | 0/114. 09 |114.38 | 130.2 | 145.5 | 160.0 | 173.8 | 186.7 | 198.9 | 210.4 | 221.3 | 231.6 | 241.4 | 250 6 | 280.4 

CuHis | 0/123. 28 123.61 | 141.3 | 158.3 | 174.5 | 180.8 | 204.1 | 230.5 | 242.5 | 254.0 | 264.8 | 275.1 | 284.8 

CisHae | 0/132. 46 (132.83 | 152.4 | 171.1 | 188.9 | 205.7 | 221.5 | 236.4 | 250.5 | 263.8 | 276.4 | 288.3 | 299.6 | 310.3 

CuHm | 0/141. 64 |142.06 | 163.5 | 183.9 | 208.3 | 221.7 | 239.0 | 255.2 | 270.5 | 285.0 | 208.7 | 311.7 | 324.0 | 335.7 
2 


CyuHe | 0150.82 151.28 | 174.5 | 196.7 | 217.8 | 237.7 . 290.6 | 306 321. 335.2 | 348.5 | 361.1 
| | | 


' | | 

| CiuHine | 0,160.01 160. 50 | 185.6 | 209.5 232.2 | 253.6 | 273.8 | 310.6 | 327.5 | 343.5 | 358.6 | 373.0 | 386.5 

--| CH | 0/169. 19 |169. 73 | 196.7 | 222.3 | 246.6 | 269.6 | | 330.6 | 348.7 | 365.9 | 382.1 | 397.5 | 412.0 

| CisHao | 0/178. 37 |178.95 | 207.8 | 235.1 | 261.0 | 285.6 30.2 | 350.6 | 370.0 | 388.2 | 405.6 | 422.0 | 437.4 

CwHa | 0/187. 56 |188.18 | 218.8 | 247.9 | 275.5 | 301.5 | 326. .9 | 370.7 | 391.2 | 410.6 | 429.0 | 446.4 | 462.8 

Petradecylbenzene. _ -.-.---| CooHas| 0/196. 74 ‘ie | 220.9 | 260.7 | 289.9 | 317.5 | .7 | 390.7 412.4 | 433.0 | 452.5 | 470.9 | 488.2 

entadecylbenzene......__..__......| CarEfse | 0/205. 92 |206.62 | 241.0 | 273.5 | 304.3 | 333.5 | 360.8 | | 410.7 | 433.7 | 455.4 
Heradecylbenzene.. j 0/215. 11 [215.85 | 252.1 | 286.3 | 318.8 | 349.4 | 378.2 | 405.2 | 430.8 | 454.9 | 477. 


| 
11. = 12. . 14.43) 15.97, 17.40) 18.75 
| 





| 476.0 | 495.4 | 513.7 
| 499.4 | 519.9 | 539.1 
erement per CH; group 0 9.183) 9. 204 


20.03) 21.24) 22.38) 23.46) 24.48 25.43 





* See footnote “a”’ of table 9. 
> S° is the entropy (exclusive of nuclear spin) of the given hydrocarbon in the thermodynamic standard gaseous state of unit fugacity (1 atmosphere) at the 


mperature indicated. 





TaBLe 15.—Values * of the heat content, (H° — Hj), for the 14 alkylbenzenes, Cy to Co, for the ideal gaseous state, to 1,507 








Compound (gas) ze 0 | 208 10| 300 400 | 500 | 600 | 700 800 900 | 1.000 | 1,100 | 1,200 | 1,00 | 1, ts 





Heat content, (H°—H}), in cal/mole 





Benzene... a 12285 | 16267 | 20612 | 25280 - 


Methylbenzene (toluene)... Z 4306 15334 | 20247 | 25621 hylbe 


yiben: 

imet 
Dimet 
Dimet 


BUIERIEED. cncccnnesecvecocencsesed 5335 18799 
1,2-Dimethylbenzene (o-xylene) . ...... 19070 | 24962 
1,3-Dimethylbenzene (m-xylene) . ..... 18563 
1,4Dimethylbenzene (p-xylene)... .... 18499 


2 


“opylb 
propy!! 
ethyl- 
ethyl- 
ethyl- 
3-Trit 
ie)... 
4Trit 
ne) 

Trim 


fe 22300 | 29250 
Isopropylbenzene (cumene) __ - 28900 
1-Methy}2ethylbenzene. - 22530 | 20460 
1-Methyl-3-ethylbenzene._ _ 22030 
1-Methy}-+ethylbenzene __ - 


28820 

1,2,3-Trimethylbenzene (hemimell 
tene) 21990 | 28750 
28830 
28386 


2 3 
geGl Gf 


1,2,4Trimethylbenzene 





22060 
1,3,5-Trimethylbenzene (mesitylene) - . 21623 















































* See footnote “a” of table 9. : - > 
* (7°—H,) is the heat content at the given temperature less the heat content at 0° K of the given hydrocarbon in the thermodynamic standerd p perata 


state of unit fugacity (1 atmosphere). 


TABLE 16.—Values* of the heat content, (H°—H}), for the normal alkylbenzenes, Cy to Cx, for the ideal gaseous state, 
1,500° K 





Temperature * in ° K 





Compound (gas) mula o|28.16| 200 | wo | soo | coo | 70 | sao | ooo | 1,000 1,100 | 





prene 

hylbe: 
bylbenz 
| 64130 ropylb 
77480 utylbe 


Heat content,» (H°—H}), in cal/mole 





Benzene... ses 
Methylbenzene (toluene) 
Ethylbenzene 
n-Propylbenzene 
a-Butylbenzene 


NT 


103800 1) mylbe 
lexyl be 

57960) 68900! 116800 leptylt 

64530; 76670) 129000 Detylbe 

84440 143000 | 

77650} 92210 156100) {Tia 

84220) 99980) 


3 


n-Amylbenzene - 
n-Hexylbenzene 
n-Heptyl benzene 
n-Octylbenzene 
n-Nonylbenzene 


90780} 107760) 
97340) 115530) 
103910) 123300 
110470; 131080) 
117040) 138850 


n-Decyl benzene 
n-Undecylbenzene. - - 
n-Dodecylbenzene..__ - 
n-Tridecylbenzene . _ . 
n-Tetradecylbenzene _- 








n-Pentadecylbenzene._. 
130160) 154400¢ 


n-Hexadecylbenzene._.. 








g EEE: 
eis d3088 GREG? GREEE 
































Increment per C H; group 6564| 7773 





* See footnote “a” of table 9. 
> (H° Hy) is the heat content at the given temperature less the heat content at 0°K of the given hydrocarbon in the thermodynamic standard gaseous state @ 


unit fugacity (1 atmosphere). 
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o 1,500 pie 17.—Values * of the heat capacity, Cy, for the 14 alkylbenzenes, Cy to Co, for the ideal gaseous state, to 1,500° K 
Y | 








Temperature * in ° K 





Compound (gas) 


=== [== lo] ™ |= | [ime] nm ome] oe 





Heat capacity,> C>, in cal/deg mole 





| Cole A a 19. 65 ae 41.75 | 45.06 | 47.83 | 50.16 | 52.16 






































51640 
64129 | amrnyibensene (toluene)... ---.--..-- CrHs 24.80 | 24.95 \ 40. 54 | 46.58 | 51.57 | 55.72 | 59.22 | 62.19 | 64.73 
, : , CsHw 30. 69 | 30.88 L 49.35 | 56.44 | 62.28 | 67.15 | 71.27 | 74.77 | 77.77 
7 Dimethy .-....| CeBe 81.85 | 32.02 | 41.03 | 49.11 | 55.98 | 61.76 | 66.64 | 70.80 | 74.35 | 77.40 
sad Dimethylbenzene (m-xylene) _. CsHw 30. 49 | 30.66 . 48.43 | 55.51 | 61.43 | 66.41 | 70.63 | 74.23 | 77.31 
+4 Dimethylbenzene (p-xylene) - CsHiwe 30. 32 | 30.49 \ 48.06 | 55.16 | 61.12 70. 39 | 74.02 | 77.13 
opylbenzene........---- . CHa 36.73 | 36.99 : 66.0 | 72.7 8. 83.1 | 87.1 | 90.6 
washed sropylbenzene (cumene) . CHa 36.26 | 36.47 | 48. 9 | 662 | 72.9 83.3 | 87.3 | 90.8 
nee ethy!-2-ethylbenzene - - CsHu 37.74 | 37.04 . 65.8 | 72.5 82.8 | 86.9 | 004 
eneey ethyl-3-ethylbenzene.. | CoE | 0 | 36.38 | 36.59 | 47. 2 | 65.4 | 721 | 77.8 | 827 | 96.8 | 904 
an ethyl4-ethyl benzene. . -. CHa 36.22 | 36.42 A 65.0 ‘ 4 82.4 | 86.6 | 90.2 
on 3-Trimethylbenzene (hemimelli- 
- .| CHa 36.85 | 37.04 64.0 . . 81.6 | 85.9 | 80.5 
89400) 
adeeanehoattesunntinting CHa 37.10 | 37.28 . . 64.2 > 3 b 86.0 | 89.6 
oxen 5-Trimethylbenzene (mesitylene)..| CsHi 35.91 | 36.10 . . 64. 08 . . 84 | 81. 86.07 | 89.72 














39180 





* See footnote “‘a’’ of table 9. 
> Cis the heat an at constant pressure of the given hydrocarbon in the thermodynamic standard gaseous state of unit fugacity (1 atmosphere) at the 

































































dard g perature indicated. 
sed BLE 18.—Values* of the heat capacity, Cl for the normal alkylbenzenes, Cs to Cx, for the ideal gaseous state, to 1,500° K 
—— Temperature * in ° K 
Compound (gas) 0 | 2 1| 300 | 400 | 500 | eon | 700 | 800 | 900 | 000 | 1,100 | 1,200 | 1,300 | 1,400 | 1, 0 
2 
400 ||, Heat capacity,» C$, in cal/deg mole 
a... 
Pe oe et ee _ 
a i. Cos ’ . ’ 41.75] 45.06) 47.83) 50.16) 52.16) 53.86) 55.32) 56.58) 57. 67 
hylbenzene (toluene)..............| CrHs i . . ’ 58) 51.57) 55.72) 59.22) 62.19) 64.73) 66.90) 68.77) 70.38) 71.78 
gen fhe bylbenzene . ; | CsHie ; . ' 44) 62.28, 67.15) 71.27) 74.77) 77.77| 80.35) 82.57, 84.49) 86.16 
om ropylbenzene. g-=---+--| Collis , 0} 72.7 83.1 | 87.1 | 90.6] 93.6| 96.2| 98.5 | 100.4 
—e en utylbenzene pense | \ | 75. 83.1 L 94.9 | 99.5 | 103.5 | 106.9 | 109.9 | 112.5 | 114.7 
3800 1! mylbenzene.. . Lannaee] Onsllles 54 | 62. .7 9 | 93.4 | 106.7 | 111.9 | 116.3 | 120.2 | 123.6 | 128.5 | 128.9 
erylbenzene = CuHis ' ’ ; 4 | 103.8 | 118.4 | 124.2 | 129. 2 | 133. 5 | 137.2 | 140.4 | 143.2 
eon | ylbenze seveeeeeeeee=| Cras . 114.2 | 190.2 | 198.6 142.1 | 146.8 | 150.9 | 154.4 | 157.4 
coon b Petylbenzene. . -| Cisse | . 124.5 | 134.0 | 142.0 | 149.0 | 154.9 | 160.1 | 164.6 | 168.4 | 171.7 
~ . Nonylbenzene | Cis | 71. ; . 134.9 1 | 153.8 | 161.4 | 167.8 | 173.4 | 178.2 | 182.4 | 185.9 
9200 IMDecylbenzene.... lea ; ; ; .2 | 145.3 165.6 | 173.7 | 180.6 | 186.7 | 191.9 | 196.4 200. 2 
---------| CrrHie . ; ’ .7 | 155.6 4 |r 186.1 | 193.5 | 200.0 | 205.6 | 210.4 | 214.4 
-_ = z Ses | CisElse . .2 | 166.0 199.2 | 198.5 | 206.4 | 213.3 | 219.3 | 224.4 | 228.7 
ce y Be swacceace-.| Opis . . .6 | 176.4 201.0 | 210.8 | 219.2 | 226.6 | 232.9 | 238.4 | 242.9 
3500) 21 es ....| CroHs ’ 1 | 186.7 212.8 | 223.2 | 232.1 | 239.8 246.6 | 252.3 | 257.2 
1600 24% | wn. 
1700, SEP entadecylbenzene aimed 105. .6 | 197.1 | 212.0 | 224. 4 235.6 | 244.9 | 253.1 260.3 | 266.3 | 271.4 
| 207.5 | ea oem | 248.0 | 257.8 | 266.4 | 273.9 saael Fa 
60) SMB rement per CH group............. Risdieal 5. } : . ws 11. . 11.79) 12.37! 12. 86| 13. 20 13. 67 — 14. 25 
| | 
ial aS | 
an! * See footnote “a” of table 9. 
> C; is the heat capacity at constant pressure of the given hydrocarbon in the thermodynamic standard gaseous state of unit fugacity (1 atmosphere) at 
onl te mpe rature indicated. 
us 


search Igeee@*s and Free Energies of Alkylbenzenes 


IV. Heat of Formation, Free Energy of Formation, and 
Equilibrium Constant of Formation 


1. Method and Data Used in the elements, solid carbon (graphite) and gas 
Calculations hydrogen, all at the given temperature, are p 
sented in tables 19 to 24, which give values of; 
heat of formation, AHf°, the free energy of fom 
tion, AFf°, and the logarithm of the equilibria 
constant of formation, logy Af, to 1,500° K. 
Figure 1 shows the thermodynamic stability 4 
the normal monoalkylbenzenes in the gas 
state as a function of temperature, in the form# 
a plot of the standard free energy of formatiy 
> sae per carbon atom, divided by the absolute tempe 
The resulting values of the thermodynamic prop- ture. This plot may be compared with com 
erties for the formation of the 14 alkylbenzenes _ sponding plots for the normal paraffins {13}, i 
through C,H, and of the higher normal mono- _ normal 1-alkynes [15], and thenormal 1-alkenes |) 
alkylbenzenes, in the gaseous state, from the 


The same method of calculation was used as 
described in Section IV—1 of reference [1]. 

For the heats of formation at 25° C of the alkyl- 
benzenes in the gaseous state from the elements, 
solid carbon (graphite) and gaseous hydrogen, the 
values given in reference [2] were used. 
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Figure 1.—Thermodynamic stability of the normal monoalkylbenzenes in the gaseous state as a function of temperawt 


The scale of ordinates gives the value of (1/n) (AFf°/T) in calories per degree mole, where n is the number of carbon atoms per molecule, 7 is the absolute & 
perature in degrees Kelvin, and AFf° is the standard free energy of formation of the hydrocarbon from the elements, solid varbon (graphite) snd gaa 
hydrogen, all at the given temperature, as given in tables 15 and 16. The scale of abscissas gives the temperature in degrees Kelvin. 
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V. Free Energies and Equilibria of Isomerization 


rom the values in table 21 and 23, calculations 
» made of the values of the free energies and 
librium constants for the isomerization of the 
» and the C,H» alkylbenzenes, according to 
reactions 
(ethylbenzene, gas) C,H,» (isomeric alkyl- 
benzene, gas) (13) 
(n-propylbenzene, gas)=C,H,, (isomeric 
alkylbenzene, gas). (14) 
e resulting values are given in tables 25 and 


under the following headings, for the isomeri- 
ion reaction as written: A F°/T, the standard 


10 


12 


BLE 25. 


free-energy change divided by the absolute temp- 
erature; and K, the equilibrium constant. In 
tables 27 and 28 are given values of N, the mole 
fraction of the given isomer present at equilibrium 
with its other isomers. For any two isomers, the 
ratio of the corresponding values of K (or of N) 
in tables 25 and 27 gives the ratio of the amounts 
of those two isomers present at equilibrium with 
one another in the gas phase at the given tempera- 
ture. For the purpose of retaining the significance 
of their change with temperature, the values in 
tables 25 to 28 are written with more figures than 
are warranted by the absolute uncertainty. 


Values of the free energies and equilibrium constants for the isomerization of ethylbenzene to the isomeric Cg 


alkyl benzenes in the ideal gaseous state to 1,500° K 








Ethylbenzene=o-xylene 


Ethylbenzene= m-xylene 


Ethylbenzene= p-xylene 





A F°/T 


4SF°/T xK 


AF°/T 





® 


—6. 813 
—6. 760 
—4. 645 
—3. 403 
—2. 572 


: 





—1. 962 
—1. 489 
1.113 
—0. 809 
—. 555 
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—. 145 
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- 145 
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+. 208 
- 453 
- 662 
. 836 

















BLE 26.—Values of the free energies and equilibrium constants for the isomerization of propylbenzene to the isomeric Cy alkyl 
benzenes in the ideal gaseous stale to 1,500° K 
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Tasie 27.—Values of the equilibrium concentrations for the isomerization of the Cy alkylbenzenes in the ideal gaseow 
to 1,500° K 








Composition, in mole fraction, of equilibrium mixture of isomers 





Ethylbenzene 


oXylene 


m-Xylene p-Xylene 





% > 
> 


seeez2 e888 


-_-~ 





gees 




















TaBLe 28.—Values of the equilibrium concentration for the isomerization of the Cy alkylbenzenes in the ideal gaseous us 
to 1,500° K 
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In figures 2 and 3 are plotted, as a function of 
the temperature, the values of AF°/7' given in 
tables 25 and 26 for the isomerization of the 
C,H and the C,H,» alkybenzenes. These plots 
may be compared with corresponding plots for 
the paraffins [21, 22], the acetylenes [15], and the 
monoolefins [24]. From these charts, one may 
see at a glance, for any temperature in the given 
range, and within the limits of uncertainty of 
the calculations, which of the isomers is thermo- 
dynamically most stable (lowest value of AF°/T) 
and which is the least stable (highest value 
of AF°/T). 

In figures 4 and 5 are plotted, as a function of 
temperature, for the C,H, and the C,H, alkyl- 
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benzenes, respectively, the amounts, in m0 
fraction of each of the isomers present at equ 
librium with its other isomers in the gas pha 
as given in tables 27 and 28. The vertical witt 
of each band gives the mole fraction for tb 
isomer at the selected temperature. The mé 
fractions of the several isomers are plotted adé 
tively, so that their sum is unity at all 
peratures. 

In table 29 the directly measured equilibrium 
compositions for several reactions [9, 25] are com 
pared with those calculated in this paper. 
values agree within the combined limits of ¥ 
certainty of the experimental and calculate! 
values. 
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Figure 4.—Equilibrium concentrations of the CsHy 
alkylbenzenes. 


The scale of ordinates measures the amount in mole fraction, and the scaie 
of abscissas gives the temperature in degrees Kelvin and degrees centigrade. 
The vertical width of a pand at a given temperature meesures the mole 
fraction of the given isomer present when at equilibrium with all of its 
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IGURE 3.—Free energy of isomerization of the CH» FIGURE 5.— Equilibrium concentralions of the C,Hy 
alkylbenzenes. 


alkylbenzenes. 

The scale of ordinates measures the amount in mole fraction, and the scale of 
ae scale of ordinates gives the value of AF?/T, in calories per degree mole, abscissas gives the temperature in degrees Kelvin and degrees centigrade. 
for the isomerization of n-propylbenzene into the other isomers, in the The vertical width of a band at a given temperature measures the mole 
sous state, as indicated. The scale of abscissas gives the temperature in fraction of the given isomer present when at equilibrium with all of its 
degrees Kelvin. other isomers, in the gas phase. 
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TaBite 29.—Comparison of calculated and experimental 
equilibria for the alkylbenzenes in the liquid state at 
823° K 








Equilibrium composition in mole percent 





Component 
Experimental ® Calculated > 





Three xylenes: 
o-Xylene (liquid) _.......... 
m-Xylene (liquid) 
p-Xylene (liquid)-.......... 


Benzene, toluene, and xylenes 
(initial com position): 
Benzene (liquid) =50....... Aa) . 57.5 
Toluene (liquid) =50___.._- 4i 35.0 
Total xylenes (liquid) =0_.. 4.7 7.5 


Benzene (liquid) ~55_. __. 58.5 54.5 
Toluene (liquid) =36__ .__- 365 37.0 
‘Total xylenes (liquid) =9 _. 5.0 . 8.5 














* The data in the first column are those of Pitzer and Scott (9); the data in 
the second column are those of Norris and Vaala [25], with the uncertainties 
estimated by the present authors. 

> These values are cm iculsted from the thermodynamic functions given 
in this report and the vapor pressure data given in the tables of the American 
Petroleum Institute Research Project 44 [26]. 
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VI. Equilibrium Constants for Some Reactions Involving 
Alkylation, Cyclization, and Trimerization 


Table 30 presents numerical values of the plotted, as a function of temperature, values o 
equilibrium constant, and of its logarithm, for the Jogarithm of the equilibrium constant for thes 
some reactions of alkylation (addition of an olefin _reactions of alkylation, cyclization, and trimer 
to benzene to form an alkylbenzene), cyclization zation, respectively. The effect of changig 
(conversion of normal paraffin to alkylbenzene plus isomers, either reactants or products, can & 
hydrogen), and trimerization (of acetylene to determined readily from the free energies ani 
benzene and of methylacetylene to 1,3,5-tri- equilibrium constants of isomerization previous 
methylbenzene). In figures 6, 7, and 8 are given. 
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Study of the Damaging Effect of Fatigue Stressing 
on SAE X4130 Steel 


By John A. Bennett 


Rossini 


Ml Bureay 


The damaging effect of fatigue stressing above the endurance limit was investigated with 
notched specimens of SAE X4130 steel. The damage was measured by the decrease in 
endurance at another stress. A deflection method for detecting the formation of a fatigue 
crack permitted the damage measurement to be limited chiefly to the precrack stage. The 
results showed that the apparent rate of damage depends on the stress history. If the 
prestress is higher than the test stress, the damage occurs rapidly at first, then more slowly. 
The reverse is true if the damaging stress is lower than that used to measure the damage. 

Tests also were made with smooth specimens in an effort to determine the cumulative 
damage caused by fatigue at more than one stress. Two methods were developed for 
extrapolating to determine the point at which fatigue cracking starts, so the damaging cycle 
ratios could be based on the life of the specimens prior to cracking. Complete S-N curves 
were determined for specimens after each of eight different damaging treatments. With 
these curves a method of expressing damage was developed that permitted the direct addi- 
tion of damage occurring at different stresses. The reliability of this method was checked 
by testing specimens after fatigue loading at two or more different stresses and comparing 
the results with predictions based on the addition of the indicated damage. The agreement 


was within the experimental error. 
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maging Effect of Fatigue Stressing 





I. Introduction 


Most ferrous metals exhibit a fatigue limit. If 
a sample of the metal is subjected to a pulsating 
stress greater than the fatigue limit, a crack will 
form that eventually will lead to fracture. Prior 
to formation of the first crack a change is taking 
place in the metal that is detectable only by 
means of a fatigue test. For example, if a speci- 
men is subjected to an overstress for one-half the 
number of cycles required to cause fracture, the 
mechanical properties of the material, other than 
fatigue, appear to be unchanged. However, the 
fatigue properties, both the endurance at stresses 
above the fatigue limit and the fatigue limit itself, 
will be lowered. 

The purpose of the work described in this 
report was to study some of the phenomena 
associated with this damaging process, particu- 
larly with a view to making the information more 
readily applicable to steel structures against 
failure under alternating loads. 

Previous investigations of fatigue damage have 
dealt largely with smooth specimens, although the 
majority of fatigue failures in service propagate 
from points of stress concentration. Therefore, 
it was thought worthwhile to make a thorough 
study with notched specimens for comparison 
with the results with smooth specimens. This 
study comprised the first part of the present inves- 


tigation, and the results are discussed in go 
IV of this report. The work was done unde; 
auspices of and with financial assistance from; 
National Advisory Committee for Aeronay 
and has been described in a publication of 
committee [1].* 

A large proportion of the highly stressed mg: 
elements in service probably are stressed 4) 
their fatigue limit, occasionally during their J 
It is important, therefore, to know how » 
damage is done by this overstressing. As the 
of damage in a simple test is known [2] to 
complicated function of the stress history, jt | 
not been possible to predict systematically 
performance of material in service where the st 
variation follows a random schedule. The seq 
part of this investigation dealt with the study 
the damage caused by combinations of diffe 
ranges of stress. It was hoped that the labors 
tests could be made to conform more closely 
the random fluctuation of stress usually encog 
tered in service by making the loading sched 
more and more complicated. The results of i 
part of the investigation are described in section 
of this report. 


! Figures in brackets indicate the literature references at the en! 
paper. 


Il. Materials and Specimens 


All the fatigue tests included in this report were 
made on R. R. Moore machines. The specimens 
were of normalized SAE X 4130 steel. A different 
heat was used for each part of the investigation. 
The mechanical properties of the steels are listed 
in table 1. 


TaBLeE 1.— Mechanical properties of the steels 








| M 390 used 
for notched 
| specimens 


M 357 used 
for smooth 
specimens 





Yield strength, 0.2-percent offset, Ib/in.?_. 63, 000 
Ultimate tensile strength, Ib/in.*. _ . 98, 000 
Elongation in 1 in., percent 29 
Reduction of area, percent wasnbd 60 
Rockwell hardness, B scale ide 89 


62, 500 








The notched specimens used in the tests 
scribed in section 1V had a cylindrical test sect 
0.35 in. in diameter, into which was cut a 
circular circumferential notch 0.05 in. deep 
0.05 in. in radius. The stresses listed for t 
specimens were calculated for the smallest sect 
0.25 in. in diameter, without regard to st 
concentration. The notches were finished wit 
copper wire, slightly smaller in diameter than! 
notch, charged with a slurry of No. 302 emey 
water. The specimens were rotated slowly = 
lathe while the wire, held at right angles ' 
axis of the specimen, was rotated rapidly. 1 
the fine polishing scratches at the bottom of! 
groove were substantially parallel to the axis 0! 
specimen. 





tically 4 


» the str 
"he see 


es 


closely 
ly encog 
y sched 


Its of t 


) Section 


om off 
axis off 


he smooth specimens used in the tests de- 
‘hed in section V had a minimum section 0.25 
in diameter, with a contour radius of 9% in. 
ey were polished with Aloxite polishing paper 
ts mounted on a rubber drum 5% in. in diam- 


+. The specimens were mechanically rotated 


Ficure 1.—Specimen comparator. 


out 15 rpm while being held lightly against the 
tating drum, the axis of the specimen being at 
zht angles to the axis of the drum. The drum 
eed was 1,750 rpm, and the specimens usually 
buld be polished satisfactorily in less than 1 
inute. The only precaution necessary was to 


avoid heavy pressure of the specimen against the 
abrasive. This polishing method has proved 
economical in preparing standard R. R. Moore 
specimens. Although the resulting surface is not 
highly polished, the results obtained with these 
specimens are as consistent as those with a higher 
polish. All the specimens were examined with a 
low-power microscope to make sure that none 
with circumferential scratches were tested. 

The minimum diameter of the specimens was 
measured, with a probable error of about 0.0001 
in. For the notched specimens this was done 
with the aid of calibrated gage wires slightly less 
than 0.1 in. in diameter placed in the groove, the 
distance between the outsides of the wires being 
measured with micrometer calipers. The smooth 
specimens were measured with a vertical com- 
parator using @ Mo,ooo-in. dial gage. A sketch of 
this instrument is shown in figure 1. In place 
of the usual anvils, contact with the specimen was 
made with two ball bearings about 1-in. outside 
diameter. The specimen was placed between the 
bearings with the axis of the specimen at right 
angles to those of the bearings. In this way the 
specimen could be moved back and forth parallel 
to its length to locate the minimum section with- 
out danger of scratching the surface. 

For some of the work described in section V it 
was necessary to have a reference mark at the 
minimum section of the specimens. To locate this 
point, a light source was set up so that an edge 
of the projected beam coincided with the line 
between the centers of the comparator bearings. 
Thus when a minimum reading was obtained on 
the dial, the edge of the beam passed across the 
minimum section of the specimen. A reference 
mark was made at this position with a very small 
drop of rubber cement placed on the specimen 
with a fine wire. 


Ill. Testing Machines and Special Apparatus 


The R. R. Moore machines used in this investi- 
tion were equipped with ball-bearing spindles, 
hd all tests were run at 3,600 rpm. A variable 
totransformer was used to reduce the speed of 
he 10,000-rpm motors with which some of the 
achines were equipped. 


1. Deflection Measuring Micrometer 


This has been described in detail in a previous 
ublication [3]. Briefly, the upper one of a pair 
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of contacts was fastened to the specimen end of a 
bearing housing on the R. R. Moore machine, 
the lower contact being carried on a micrometer 
screw mounted on the bed of the machine. The 
contacts operated a signalling device through a 
tube circuit, so the position of the upper contact 
could be determined by raising the lower contact 
and noting the reading of the micrometer drum 
when the circuit closed. The contacts were set 
slightly apart so that any increase in the deflec- 
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tion of the specimen during test would close the 
circuit, operating the signal. The best results 
with this apparatus were obtained by mounting 
the fatigue-testing machine on springs to minimize 
vibrations from extraneous sources. No specimen 
was used that caused excessive vibration when 
the machine was running. With these precautions, 
it was possible to detect a deflection of 0.001 to 
0.002 mm. 


2. Stroboscopic Viewing Device 


One of a pair of contacts was mounted in a 
Bakelite disk and fitted to the shaft of the R. R. 
Moore machine, the other (stationary) contact 
being mounted so as to be rotatable about the 
center line of the machine. These contacts were 
used to flash a stroboscope, and by adjusting the 
stationary contact circumferentially the flash 
could be made to occur at any position in the 
revolution of the specimen. The flashing light 
was directed, by means of a condenser and mirror, 
onto the bottom surface of the specimen, and a 
low-power microscope fitted with a prism in front of 
of the objective was used to view the illuminated 
area. Because this area was stressed in tension, 
the fatigue cracks were opened and made readily 
visible. 


3. Apparatus for Counting and Measur- 
ing Fatigue Cracks 


A traveling microscope was used for this work, 
figure 2. The specimen was removed from the 
fatigue machine and mounted in spindles on the 
stage so that it was parallel to the traversing screw 
of the microscope. Weights hung on the ends of 
the spindles produced a bending moment which 
stressed the upper surface in tension. The area 
being examined was illuminated by an intense 
beam incident at a grazing angle in the plane of 


IV. Test Results with Notched Specimens 


Previous work, both at this Bureau and else- 
where, revealed two significant factors for satis- 
factorily evaluating fatigue damage. First, the 
dispersion of the experimental results necessitates 
an average of about eight specimens to give a 
sufficiently precise result for any one condition. 
Second, the evaluation of damage by running tests 
to complete fracture is limited to cycle ratios ? 


2“Cycle ratio” is defined as the ratio of the number of cycles run at a given 
overstress to the number necessary to cause failure at that stress. 
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the axes of the specimen and microscope. | 
scratches seen under the microscope were aly 
this plane, little light was reflected into the nj 
scope. A fatigue crack, being at right angls 
the incident beam, showed up brightly againg 
generally dark background. 
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Ficure 2.—Equipment for measuring and counting 
fatigue cracks. 













The eyepiece of the microscope was fitted wi 
both a cross hair and an adjustable mask. § 
measuring the length of cracks, a subdivided dm 
was fitted on one of the spindles. The cent 
angle subtended by the crack was obtained | 
rotating the specimen under the cross hair 
noting the drum reading at each end of the en¢ 
For counting cracks the mask in the eyepiece 
set so that the field observed on the specimen 
just one millimeter wide (the edges of the 
were at right angles to the length of the specim 
The surface could then be scanned in 1 mm wi 
bands, starting from the reference mark at 
minimum section, by traversing the micros 
1 mm after each revolution of the specimen. 
ratchet was set to give a click to indicate 
completion of each revolution. 


















sufficiently small so that the chance of crack fom 
tion during prestressing is negligible; otherm 
the damage in some cases would be dependé 
only on the size of an initial fatigue crack mi 
than on the effect of prestressing. 

The methods followed in the study of fati 
damage with notched specimens were chosel 
the basis of these considerations. In order 
guard against crack formation it was necessaly 
have a means of detecting fatigue cracks 
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stage without stopping the test. Such a 
hod was developed based on the fact that the 
tion of a rotating beam specimen under 
stant load increases when a crack forms. The 
sratus and technic used to detect and measure 
deflection are described in section III of this 
rt. The deflection is a function of the size 
rack, and the criterion of failure chosen was a 
ection of 0.005 mm, corresponding to a crack 
, about 12-percent of the original cross sec- 
(3]. The specimens were actually run until 
deflection had increased to 0.01 mm in order 
nake certain that the deflection was due to a 
k. but the data in this section are based on 
number of cycles run at the time the specimen 
ected 0.005 mm. This number is referred to 


Ar 


ests were conducted in the conventional man- 
(to fracture) on a group of notched specimens 
etermine the S-N curve and the fatigue limit. 
r stresses in the range between the fatigue 
it (+39,000 lb/in.? for notched specimens) and 
yield strength of the steel (62,500 Ib/in.*) 
e used in the investigation. The four stresses 
sen were +42,000, +48,000, +54,000, and 
0,000 Ib/in.? 
ests were made to determine N, at each of 
se four stresses. At least 10 specimens were 
ed at each stress, and the median of each group 
sused as the value of endurance at that stress. 
»median was considered to be more representa- 
p of a group of endurance values than the 
rage, as it is less affected by an occasional exces- 
ly divergent value. Also the median can be 
mined in cases where some elements of the 
up cannot be expressed numerically but are 
bwn to lie in a certain range. For example, in 
emining the endurance at a given stress, if 
y one specimen of the group did not fail, the 
rage endurance of the group would be infinite, 
reas the median value would be finite. 
is the scatter of the results was not uniform 
he four groups, it was necessary to test more 
‘imens in some cases than in others in order to 
in approximately the same precision for each 
The number of tests in each group 
ged from 11 to 17. 
he values of N, (the number of cycles to form 
ack of definite size) for each selected stress are 
ed in table 2. 
i the next stage of the investigation each speci- 
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men was first stressed for a predetermined number 
of cycles at one of the four chosen stresses, then 
the stress was changed as given in table 3 and the 
number of cycles to failure determined. The dif- 
ference between this value and the median (see 
footnote 3) value of endurance for the original 
material at the same stress (table 2) gave a 
measure of the damage caused by the prestressing. 
This difference was expressed as a percentage of 
the undamaged endurance. 

For each of the combinations in table 3, tests 
were made with the prestressing carried to 10, 25, 
50, 75, and 90 percent of the median value of N.. 
The number of tests made for each prestress con- 
dition was between 6 and 10, depending on the 
scatter of the individual values. The median 
values for each group are listed in table 4.5. The 
variation of damage with percentage of prestress 
is shown graphically in figures 3 and 4 for a test 
stress of +48,000 lb/in.?, and in figures 5 and 6 
for a test stress of + 54,000 Ib/in.? In the graphs 
the broken lines join the upper and lower quartile 
points of each group of values and thus give an 
indication of the scatter of the data. The curves 
from figures 3 to 6 are combined in figure 7 for 
comparison. 

In the test groups in which prestressing was 
carried to 75 and 90 percent of the endurance, 
some specimens failed before the prestressing was 
completed, that is, the damage was greater than 
100-percent. The inclusion of these results caused 
no difficulty in determining the median, although 
it did make the comparison of precision between 
the groups more uncertain. 

While the results of the damage tests as shown 
in figure 7 are not directly comparable with those 
obtained with smooth specimens [1, 2], there are 
certain similarities. The tendency with both 
types of specimen was for the curve representing a 
prestress above the test stress to lie above the 45° 
line. The opposite tendency was noted when the 
magnitude of the prestress was below that of the 
test stress. (The 45° line may be considered as 
the damage curve for a prestress equal to the test 
stress.) 


? The middle number of a series arranged in order of magnitude is the 
median. The middle one of the numbers that lie below the median is the 
lower quartile, the middle one of those above is the upper quartile. Half 
the difference between the upper and lower quartiles is the semi-interquartile 
(SIQ) range. This is a measure of the dispersion of the individual values 
which expresses the same limits relative to the median as the probable error 
expresses relative to the average. 
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ne 7.—Damage as measured by loss of endurance (N,). 


from figures 3 to 6, notched specimens. P= Prestress in pounds per 
square inch. 7T'=Test stress in pounds per square inch. 


TaBLe 2.—Endurance of notched specimens 





Number of 
cycles to crack 








I 
i 
Stress, Ib/in.? 
| 
! 


963, 000 
264, 000 
93, 000 
44, 000 


+42, 000 


+54, 000 
60, 000 





i 
| 





TaBLE 3.— Damage test conditions 


——= 


| Prestress, Ib/in.? 








Test stress 
Ib/in.? 





+42, 000 
+54, 000 
+48, 000 
+60, 000 


yr 


— 











order to provide a basis of comparison with 
esults of fatigue damage reported by others 
lly based on the fatigue limit), some tests of 
type were conducted with two values of pre- 
s, +42,000 and +60,000 lb/jin2 The effect 
ious amounts of prestressing on the fatigue 
subsequently determined by fracture tests is 

in table 5 and figure 8. The number of 
s of prestress used for the 75-percent and 90- 
at cycle ratios was the same as the number 
in the previous phase of the investigation, 
correction was made in determining the cycle 
stepresented by these amounts of prestress- 
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ing. The method of making this correction will be 
shown for the case of the highest cycle ratio at 
+42,000 Ib/in.? 


TABLE 4.—Effect of prestress on endurance of notched speci- 
mens of SAE X 4130 steel stressed as rotating beams 








| 
Cycle | 
- Endur- 
ratio, per- 
ent ot | ance, 10% 


| N. value cycles 


Num- | 


ber of 
speci- Test stress | Prestress 





+48, 000 
+48, 000 


48,000 | 
+48, 000 


+48, 000 
48,000 | 
+48, 000 
+48, 000 


+54, 000 
+54, 000 
+54, 000 
+54,000 | 
+54, 000 
+54, 000 


+54. 000 
+54, 000 
+54,000 | 
+54, 000 
+54, 000 
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Ficure 8.—Effect of prestress on fatigue limit. 
Notched specimens. 





Taste 5.—Effect of prestress on fatigue limit of notched 
specimens of SAE X 4130 steel stressed as rotating beams 








| 
| Prestress +42,000 Ib/in.* Prestress +60,000 Ib/in.* 


| Cycle | 
| ratio 





Damage, 
| Percent- 
| age of | 
| maxi- | 
| mum 


| Damage, | | 
| Pezeent” || cyete | 
| maxi- || ‘atlo 


mum 


Fatigue 
limit 


Fatigue 
limit 








blina | 
+38,000 | 
+35,000 | 
| +31, 000 
+26, 500 
| +28. 500 


lbjia. 
| +39, 500 
+39,000 | 
| +37,000 | 
| 37,500 | 
| +38, 500 
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The median N, for this stress was 963 10°. 
When it was desired to prestress a group of speci- 
mens to 90 percent of this value, or 867X10* 
cycles, some of the specimens failed before pre- 
stressing was completed. Only those which did 
not fail in prestressing could be used for the 
determination of fatigue limit, so this determina- 
tion was not truly representative. Therefore, in 
the original series of values of N,, the median of 
all) values greater than 867X10* was taken as 
being representative of the specimens which were 
used to determine the fatigue limit after pre- 
stressing. The cycle ratio as plotted in figure 8 
is based on this second median (1,057 X10 cycles), 
and is consequently less than 90 percent. A 
similar correction was made for the 75 percent 
prestress at +42,000 lb/in.*, and for the 90 percent 
value at +60,000 lb/in.2 The corrected values of 
cycle ratio at +42,000 lb/in.? were 68 percent and 
80 percent; at +60,000 lb/in.? the highest value 
was 77 percent. The slight repair of damage 
above 50 percent cycle ratio during stressing at 
+42,000 lb/in.’, if real, is surprising, and should 
be verified by additional work. 

In order to make the data on the effect of 
prestressing on fatigue limit (fig. 8) more directly 
comparable with those on the effect of prestress- 
ing on endurance (N,, fig. 7) an attempt was made 
to determine the fatigue limit of cracked speci- 
mens. The accurate calculation of stress on a 
specimen containing a fatigue crack is virtually 
impossible because of the irregular and unpre- 
dictable shape of the crack. A rough estimate of 
the increase in stress with spreading of the crack 
was made by assuming the stress inversely pro- 
portional to the ratio of uncracked area to original 
area. 

The specimens used for this determination had 
deflected 0.010 mm in the fatigue test. This 
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corresponded to a crack area equal to 17 pene, 
of the original minimum section and the fatig, 
limit of the cracked specimen was found to } 
+23,500 lb/in.? In figure 9 the data of figur 
have been replotted, with the damage express 
as percentage of the damage caused by a crack, 
this area ratio. The straight line drawn throw 
the data points for + 60,000 Ib/in.? was determi 
by the method of least squares. 
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Ficure 9.—Damage as measured by loss of fatigue |i 
notched specimens. 


L,~=tatigue limit of original specimen. L,—Fatigue limit of prestrese 
specimen. L,=fatigue limit of cracked specimen. 


Figure 10 shows the S-N curves for specime 
with three types of stress concentration: smo 
(zero concentration), notched, and cracked (ma 
mum concentration). The points marked Y 
+ in figure 10 are the results from cracked spe 
mens having crack areas other than 17 percel 
It will be noted that the value for a specim 
having an area ratio of 50 percent lies very 4 
above the 17 percent ratio curves, while the val 
for specimens having a 12 percent ratio are 
the same range or lower. This suggests that! 
stress concentration at high stresses is less 
large than for small cracks. More experime! 
work would be required before any general st 
ment could be made with certainty. 

In table 6 are shown the values of stress 
centration factors given by the ratios of the {att 
limits of the specimens under three conditiom 
stress concentration. Almen [4] has shown 
the slope of the falling part of the S-N 
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reases With increasing stress concentration. If 
» slope of the falling part of the S-N curve is a 
sure of effective stress concentration, then the 
tios of the slopes of the three lines of figure 10 
ould give the same values for these factors as 
tained from the fatigue limits. There was 
pstantial agreement of these factors as deter- 
ined by the two methods (table 6). The value 
theoretical stress concentration factor given in 
ble 6 for the notched specimen was calculated 
om a table given in reference [5]. 


BLE 6.—Stress concentration factors in specimens of SAE 
X4130 steel stressed as rotating beams 
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Ficure 10.—S-N curves for SAE X4130 steel under three 
conditions of stress concentration. 


Large open points represent the median values from a group of specimens 
tested at the same stress. Small closed points are the results of individual 
tests. 


V. Test Results With Smooth Specimens 


. Methods of Estimating the Start 
of Cracking 


he preceding experiments showed that fatigue 
mage results are more reproducible when based 
the number of cycles to form a small crack 
n when based on the number of cycles to frac- 
¢. It was considered that even greater repro- 
ibility might be obtained if it were possible to 
emine the number of cycles at which a crack 
tstarts. Then the damaging cycle ratios could 
based on the precrack stage oniy. This was 
le in the experiments described below, all tests 
ig made on smooth (9%-in. contour radius) 
imens. 
i the discussion of methods for estimating the 
of cracking, the relationship shown in figure 
should be kept in mind. This shows a semi- 
hithmic plot of test stress against the number 
cracks formed prior to fracture. The points 
ted are based only on cracks lying in a band 
wide at the center of the specimen. The 
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number of cracks increased rapidly as the stress 
was raised above the fatigue limit. Because of 
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Ficure 11.—Influence of applied stress on the number 
of cracks formed prior to fracture. 


the large variation in the number of cracks formed, 
two different methods were used for estimating the 
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start of cracking, depending on the number of logarithmic base. This is shown graphically 
cracks observed. figure 13 where (L—C) is plotted logarithmi. 
By means of the stroboscopic apparatus de- against N for one of the curves shown in figury 
scribed in section III, the surface of a specimen When curves similar to figure 13 were plotte 
was observed during the tests. The field of the various nominal stress values, it was found 
microscope was large enough so that the chance of — the slope was a simple function of the stres, 
a crack forming outside the field was small. By shown in figure 14. 
adjusting the stationary contact, any point on the 
circumference could be brought into view, thus 
permitting observation of the initial crack irrespec- 
tive of circumferential position. A small crack 
could be seen readily when in the bottom posi- 
tion because it was opened up by the tensile stress. 
As soon as a crack was observed, the test was 
stopped and the circumferential length of the 
crack was measured. This procedure was re- 
peated at frequent intervals until fracture oc- 
curred. Figure 12 shows the growth of two cracks 
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Figure 13.—Semilogarithmic plot of crack growth cu 
The ordinate is the are length (L) of the crack minus a constant (( 
figure | 
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10,000 CYCLES 
Figure 12.—Curves showing change of crack length at each 
of two different stresses. 
The midpoint of the crack when first measured is taken as zero. 
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each under a different stress. It will be noted that 
the curves appear to start from a finite value of 
crack length, rather than from zero. This is a 
surprising observation, and one which should be i | 
checked further, but in the present investigation Bary “s 
there appeared to be no question of its validity, owmnes, wanf, taveeens 
so all extrapolation was performed on this basis. Fiaure 14.—IJ nfluence of stress on the rate of crack gr 
The fact that small, though numerous cracks 
frequently have been observed at stresses consider- 
ably in excess of the fatigue limit suggests that the | to determine by extrapolation the point at w! 
large value of the initial crack length shown in cracking starts, assuming that the cracks ¢ 
figure 12 is associated with stresses nearer the exponentially from an original arc length of 
fatigue limit. 6°, as indicated in figure 12. The defle 
It was found that the crack length, ZL, could be = measuring apparatus described in section Ii} off 
expressed by the equation report was used to indicate the presence 
crack, but this did not operate until the 
log (L—C)=aN, (1) reached 50° or 60° in the circumferential lex 
where C and a@ are constants, the latter dependent the growth curves (fig. 13) being used to & 
on stress; N is the number of cycles, and A the _ polate back to the start of the crack. 
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Curves of the type shown in figure 13 were 
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The above method was used only in cases where 
e total number of cracks was not more than five, 
nce more numerous cracks might change the 

ace stress sufficiently to invalidate the extra- 
lation. For specimens which had numerous 
acks the extrapolation was performed as follows: 

means of the masked eyepiece described in 
tion III, the cracks were counted in circum- 
rential bands 1 mm wide. Practically all the 
acks occurred within 6 mm of the center of the 
cimen and for the purpose of totaling the 
mber of cracks this area was divided into three 
mes. The first extended 2 mm each side of 
e center, the second from 2 to 4 mm on each 
je and the third 4 to6 mm. The average stress 
each zone was calculated from the geometry 
the specimen and found to be 99.8, 98.8, and 
0 percent of the stress at the center for the 
st, second, and third zones, respectively. The 
uare root of the number of cracks in a band 
lotted against the average stress in that band, 
ve a nearly linear relationship which could be 
trapolated back to the stress at which no cracks 
uld be formed. Curves of this type are shown 
figure 15. 


100 


© © © 
a 4 oa rr) 


STRESS, PERCENT OF MAXIMUM 


ro 
vy 


4 16 36 64 100 144 
NUMBER OF CRACKS IN 4MM BAND 


URE 15.—Results of crack counts on one smooth specimen 
at various stages between the beginning of cracking and 


'shown on figure 16, both the stress and number 
cycles being plotted logarithmically. The 
ught line through the points was determined 
the method of least squares. A similar graph 
sed on the points representing fracture is also 
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included in this figure (the points for this curve are 
not shown). 
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Ficure 16.—S-N curves for start of cracking and for fracture 
in normalized SAE X 4130 steel. 
The experimental points are for cracking only. 


2. Damage Tests, Including Strength- 
Loss Line Hypothesis 


In the next phase of the work, two stresses were 
chosen between the fatigue limit and the yield 
strength of the material, and specimens were 
stressed for various cycle ratios at these stresses. 
The stresses chosen were +48,000 and +54,000 
Ib/in.2, and the cycle ratios were based on the 
number of cycles to the start of cracking. Eight 
prestress conditions were used, as shown in table 7. 


TaBLe 7.—Prestress conditions 








Cycle ratio, 
percentage of 

















The prestressed specimens were used to deter- 
mine a complete S-N curve for each condition, a 
typical curve of this type being shown in figure 17. 
This method has been used recently by Kommers 


[2] to measure fatigue damage. These curves 
were determined in the usual manner, that is by 
tests to fracture. The experimental data are 


4 The term N, will be used to designate the number of cycles to the start of 
cracking. This value is determined from the cracking curve of figure 16. 
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[Cycle ratios are based on Nz. 
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TABLe 8.—Results of damage tests 


“FP” indicat 





that speci 


“U”" indicates that the specimen did not break.) 


fractured, 





PRESTRESS, +48,000 LB/IN.? 


PRESTRESS, +48,000 LB/IN.? 





Cycle ratio, 10 percent 


Cyele ratio, 6634 percent 











Stress, Ib/in.?| 10° cycles Stress, Ib/in.?| 10° cycles 
+54, 000 102 F +, 000 73 F 
+54, 000 70 F +54, 000 63 F 
+54, 000 «6 F +4, 000 131 F 
+52, 000 225 F +52, 000 129 F 
+50, 000 146 F +50, 000 6 F 
+46, 200 1, 187 F +46, 000 1,131 F 

+46, 200 514 F +46, 000 516 F 
46, 000 22, 339 U +44, 600 240 F 
+46, 000 27,721 U +44, 400 1,038 F 
+45, 800 44, 088 U +44, 400 20, 635 U 

+44, 200 36, 233 U 





Cycle ratio, 334% percent 


Cycle ratio, 90 percent 





+, 000 61 F 
+54, 000 81 F 
+, 000 142 F 
| 52,000 86 F 
| 50, 000 MIF 
+47, 000 313 F 
+46, 000 743 F 
+45, 300 1, 608 F 
+45, 100 22,711 U 
+45, 000 831 F 
+44, 900 40, 787 U 
+44, 800 21,942 U 
+44, 500 20, 027 U 





+54, 000 9 F 
+54, 000 “Fr 
+54, 000 58 F 
+52, 000 112 F 
+50, 000 1ilF 
+46, 000 622 F 
+44, 000 132 F 
42, 600 410 F 
42, 400 416 F 
+42, 200 #4, 516 U 
+42, 000 17, 709 U 











PRESTRESS, +54,000 LB/IN.? 


| PRESTRESS, +48,000 LB/IN.? 





Cycle ratio, 10 percent 








| Cycle ratio, 6634 percent 











Cycle ratio, 3344 percent 





+52, 000 7OF 
+50, 000 109 F 
+48, 000 267 F 
+46, 000 1,276 F 
+45, 600 327 F 
+44, 000 &53 F 
+43, 800 23, 882 U 
+43, 800 19,220 U 
+43, 600 21,451 U 





Stress, Ib/in.?| 10° cycles Stress, Ib/in.?| 10° cycles 
+52, 000 110 F || 82,000 60 F 
+50, 000 318 F || 50,000 186 F 
+48, 000 720 F ! +48, 000 132 F 
+45, 800 2,003 F || +48,000 86 F 
+45, 800 84 F || +46,000 454 F 
+45, 600 33, 016 U | +44, 000 19, 227 U 
+45, 600 4 F || +44,000 19, 133 U 
+45, 400 21,033 U || 43,000 18, 986 U 

|| 42,000 459 F 
| +40, 400 16, 101 U 





Cycle ratio, 90 percent 











| 
| 59, 000 12F 
+52, 000 “MF 
+50, 000 oF 
+48, 000 128 F 
+45, 000 641 F 
| +44, 000 162 F 
|| 43,000 25, 304 U 
| +42, 000 144 F 
+42, 000 16, 009 U 
| +41, 000 575 F 
+40, 400 176 F 
| 40, 200 249 F 
+40, 000 21, 338 U 
| +40, 000 11, 643 U 














shown in table 8. The curves for the eight 
stress conditions are shown in figures 18 and 
The cycle ratio for the largest number of cy 
should be less than 90 percent because of 
specimens which cracked before reaching 
predetermined number of cycles. However. 
like the tests described in section IV, these rs 
could not be accurately corrected for the pre 
ture failures, so no attempt was made to do w, 
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Figure 17.—S-N curve for specimens that had bees: 
jected to 3345 percent cycle ratio at + 54,000 lb/in: 
The square point is at the stress used for prestressing and the num) 
cycles obtained by subtracting the prestress cycles from expected lile « 
stress. 





FIGURE 


An additional curve is shown on figure 
representing the S-N relationship for specim 
which had been improved by understressing. 
five specimens were used to determine this ¢ 
and the prestress treatments in the five cases w 
not identical, but as the scatter appeared t 
small, the curve is probably reasonably relia 
The complete history of each of the five s 
mens is given in table 9. 

TaBLe 9.—Tests of specimen improved by understress 


[“F” indicates that specimen fractured, ‘‘U” indicates that specimen & 
fracture] 








10 * cycles 
Test stress at 
test stress 


10 # cycles 
at 
prestress 





lbjin.? 
+56, 000 
+545, 000 
+54,000 
+52, 000 
+50, 000 


Dijin? 
+45, 000 
+46, 000 
+45, 900 
+45, 900 
+45, 800 


so F 
70 F 
14 F 
48 F 
3,933 U 


21, 700 
24, B9 
18, 567 
26, 700 
31, 578 














In figure 19, a dashed line represents the § 
relationship for specimens which had been da 
100 percent, that is, specimens which wer: | 
ready to crack. This curve was obtained 
follows: Points were obtained at + 48,00) 
+54,000 lb/in.? by subtracting the numb 
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Fieure 18.—S-N curves for specimens overstressed at + 48,000 lb/in.2 and for specimens improved by understressing. 


The figures on the curves represent the cycle ratio (based on N,) to which the specimens were run in prestress. 
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Fiaure 19.—S-N curves for specimens overstressed at + 54,000 lb/in.? 


The figures on the curves represent the cycle ratios (based on N,) to which the specimens were run in prestress, 
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cycles to form a crack from the total number to 
fracture. Two additional points were obtained at 
+42,000 and +38,000 lb/in.? by testing specimens 
having very small cracks (2° to 3° long). The four 
points obtained in this way fell very close to the 
straight line shown in figure 19. There is no way 
to obtain directly the fatigue limit for 100-percent- 
damaged specimens. 
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Figure 20.—Strength loss relationship derived from curves 
of figures 18 and 19. 


It will be noted in figures 18 and 19 that the 
curves appear to radiate from approximately the 
same point. This suggests that the damage could 
be measured by a simple geometrical property of 
the S-N curves. In figure 20 the fatigue limit is 
plotted against the slope of the falling part of the 
S-N curve. Although there is considerable scat- 
ter, the points lie approximately on a straight line, 
indicating that the damage is of the same nature 
regardless of the stress at which it was inflicted. 
The horizontal line at 3.27 represents the slopes of 
the 100-percent-damage curve and the intersection 
of this line with that through the points should 
serve as a means of estimating the fatigue limit of 
the 100-percent-damaged specimens. There ap- 
pears to be no method for checking the value 
obtained. 

It will be noted that the point representing 
specimens that had been improved by understress- 
ing lies near the straight line through the damage 
points. This indicates that the changes due to 
understressing are of a similar type to those oc- 
curring on overstressing, but the results produced 
are in the opposite direction. 

If the type of damage does not depend on the 
stress at which it is inflicted, then the amount of 
damage should be represented by distance along 
the curve of figure 20, being expressed as percent 
of the distance from the value for the original 
material to that for 100-percent damage. This 
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line therefore will be referred to as the strony: 
loss line. The distance along this line is ploty 
in figure 21 against the cycle ratio for the ¢ 
values of prestress used, the distance being ny 
ured to the foot of the perpendicular from ej 
point to the line. The curves of figure 2) , 
significant as they permit the prediction of 
complete S-N curves for specimens that have beg 
damaged by any amount of prestress at +48) 
Ib/in.?, +54,000 lb/in.? or by any combintion ¢ 
these stresses. 





100 





























a 
°o 








ALONG STRENGTH LOSS LINE 
n 
° 
































Pe et 
oO 20 40 60 
PERCENT CYCLE RATIO AT PRESTRESS 





Fiaure 21.—Cycle ratio versus damage relationship for 
stresses. 


Dashed lines show method of predicting results for one compiler * 
schedule. 


3. Multiple Prestress Tests 


In order to investigate the reliability 0: | 
predictions based on the curves of figure 2! 
at the same time to study the damage cause 
more complex loading programs, the six pres! 
ing schedules shown in table 10 were used. 
each schedule, three to five specimens were ‘ 
to fracture at the test stress, and for schedul 
an attempt was also made to determine the {i! 
limit after prestress. All results of these tests 
shown in table 10, together with the predic’ 
based on figure 21. The method of making 
predictions is illustrated by the dotted line 
the figure, showing the steps involved in sche! 
V. A cycle ratio of 85 percent at +* 
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Je ratic 
percent 
caused 


veles t 
000, a 
ts the 1 


stress. 


»2 results in 44 percent damage. A 44 per- 
t damage corresponds to 69 percent cycle 
io at + 54,000 Ib/in.?, and 69-+-20=89 percent 
le ratio at 54,000 lb/in.*, which corresponds to 
percent damage. <A 69 percent damage would 
caused by 96 percent cycle ratio at +48,000 
n2, so a 4 percent cycle ratio should remain 
r the prestressing is complete. The number 
veles to the start of cracking at this stress is 
000, and 0.04 333,000= 13,000, which repre- 
ts the number of cycles to form a crack at the 
stress. In addition, it requires 78,000 cycle 
his stress for the crack to grow to fracture, so 
total predicted life under these conditions is 
00 cycles. 


TaBLe 10.—Multiple prestress tests 


Median’ Predic- 
Test results, 10° cycles to fracture 102 | tion 103 
cycles cycles 

— | | 
16] : 207| 163| 177 
65 | 57 37 | 21 57 | 56 
7 2 5 : pened 7) 34 
101 sabkibwoal 101 | 148 
<o 5 36 | 91 


490 





percent cycle ratio at +-48,000 Ib/in.? + 3344 percent cycle ratio at 
in.*, test +48,000 Ib/in.? 
percent cycle ratio at +54,000 Ib/in.? + 3344 percent cycle ratio at 
Ib/in.2, test 454,000 Ib/in 2. 
75 pereent cycle ratio at +54,000 Ib/in.* + 15 percent cycle ratio at 
*, test + 54,000 Ib/in ?. 
percent cycle ratio at +48,000 Ib/in.? + 3344 percent cycle ratio at 
» lb/in.*, test +48,000 Ib/in 2. 
5 percent cycle ratio at 448,000 Ib/in.? + 20 percent cycle ratio at 
b/in.?, test +48,000 Ib/in 2. 
15 pereent cycle ratio at +48,000 Ib/in.*, + 10 pereent cycle ratio at 
n.*, then repeat twice more. ‘Test +48,000 Ib/in.?, 
ration of fatigue limit for schedule IV. 
odieates that specimen fractured. ‘‘U’’ indicates that specimen did 


tion to those listed, three specimens cracked before the prestressing 
eted, 








Fatigue-limit | 
10 3 eveles prediction, 
Ib./in.? 


Stress, Ib./in.? 


42,00 | 4,808 U +43, 200 
442,200 | 22,715 U 
$42,800 | 56,459 U 
+43, 000 390 F 
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The results listed in table 10 show reasonable 
agreement with the predictions, considering the 
large probable error of the tests involved. A 
smal] change in the rate of damage from one speci- 
men to the next, particularly in the region where 
the curves of figure 21 are steep, may cause a 
large variation in the damage due to a certain 
prestress. Other methods tried for predicting 
these results, including the simple cumulative 
addition of cycle ratios suggested by Miner [6], 
failed to give reasonably close values for more 
than one or two of the scheduies. 


4. Repair of Damage 


Two exploratory tests were made to verify the 
work of others on the repair of fatigue damage. 
In each case the prestress was 90 percent cycle 
ratio (based on N,) at +48,000 Ib/in *. The first 
specimen was then stressed at +42,200 Ib/in? 
(just under the new fatigue limit) for 3410 ° 
cycles. When tested at +48,000 Ib/in.? this 
specimen did not fail in 4010° cycles. The 
understressing had not only repaired the prior 
damage but raised the fatigue limit above that of 
the material in its original condition. As a 
matter of interest, this same specimen was sub- 
sequently “coaxed’’ by stressing 20 to 30 million 
cycles at each of seven successively higher stresses. 
At a stress of +62,000 lb/in.? the specimen en- 
dured 2110 ° cycles, then failed after 2,500,000 
cycles at +64,000 Ib/in*. Thus a _ specimen 
damaged to the extent where the fatigue limit 
was less than +42,500 lb/in.2 was improved by 
suitable stress treatment, bringing its fatigue 
limit above +62,000 Ib/in ?. 

The second damaged specimen was normalized 
in a vacuum, reproducing as closely as possible 
the temperature cycle of the original heat treat- 
ment. When tested at +48,000 lb/in.’? the speci- 
men fractured in 126,000 cycles. The normal 
expectancy for specimens after this prestress would 
be 143,000 cycles, indicating that the heat treat- 
ment resulted in no repair of the fatigue damage. 





VI. Discussion 


The life of a specimen under fatigue stressing 
consists of two distinct phases. During the 
first, part of the material is progressively damaged 
until a small crack is formed, and during the second 
this crack grows until fracture occurs. In this 
report the term “fatigue damage” is used to 
designate only the change occurring during the 
first stage, material being considered 100 percent 
damaged at the time the first crack forms. If 
one wishes to study the first-stage damage, some 
method of determining the start of cracking is 
required. Lacking this, the effects of the second 
stage will confuse the results. Within the stress 
range used in this study, the proportion of the 
total life represented by the second stage varies 
from as much as 50 percent to as little as 10 
percent, so it can be seen that a large error would 
be involved if the two stages were not separated. 

With the SAE X4130 steel used in this investi- 
gation, the ratio of the second stage to the total 
life of the specimen was larger at stresses con- 
siderably above the fatigue limit than at stresses 
close to the fatigue limit. This is contrary to what 
certain other investigators have reported, but 
such a discrepancy may be expected, since the 
ratio depends on many factors such as the slope 


of the S-N curve for the original material, the 
of crack propagation in the metal, and the siz 
the specimen. 

The damage resulting from a particular num) 
of cycles of overstress depends on the past hist 
of the specimen. Thus, the performance can 
be predicted accurately for metals in service whig 
are subjected to various ranges of stress, unless 
complete loading schedule is known. Since ¢} 
generally is not possible, the best that can be dy 
for increasing the safety of a given steel structy 
is to apply a large number of cycles of stress bel 
the fatigue limit after each period of overs 
This usually is taken care of by the nature of ¢ 
service, but in some applications involving ¢ 
ceptionally high stresses it might be worthwhil 
to apply low stresses purposely to counteract t 
effect of unusually long periods of high stress. 

The results of the limited number of t« 
made with specimens improved by understress 
suggest that this improvement can be measur 
by distance along a strength loss line in a man 
similar to that used for measuring damage (fig. 2 
If this can be done, then the results of any load 
schedule can be predicted, regardless of wheth 
the stresses are above or below the fatigue limi 


VII. Conclusions 


The results given in this paper lead to several 
important conclusions regarding fatigue damage 
occurring prior to the start of cracking. 

L. If the damage occurring at one stress was 
determined by measurement at another, the 
apparent rate of damage depended on the value 
of both stresses. When the first was larger than 
the second, the damage occurred rapidly at first, 
then more slowly: the reverse was true when the 
first stress was the smaller. 

2. The S-N curves for specimens after various 
damaging treatments showed that with increasing 
damage, the absolute value of the slope of the 
falling part of the curves increased and the fatigue 
limit decreased. 

3. These differences in the S-N curves provided 
a measure of damage by which the damage occur- 
ring at different stresses could be directly added. 


The performance of specimens tested under ca 
plex loading schedules could be ascertained | 
the use of cycle-ratio versus damage curves 
the stresses involved. 

4. The improvement due to undersiress 
apparently can be measured in the same way 
damage due to overstress. 

5. Fatigue damage prior to the start of crac 
was repaired by understressing but not by le 
treatment. 

Certain other conclusions regarding fatigue p! 
nomena in general may be drawn from the res 
obtained. 

6. The slopes of the falling part of the 5 
curves for the same material in smooth, note 
and cracked specimens showed approximately! 
same (stress concentration) ratios as the fal 
limits. 





from | 
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The number of cracks formed in a specimen 
ore fracture increased rapidly as the test stress 
s raised above the fatigue limit. 

s The length, Z, of fatigue cracks may be 
presented by the equation 


log (L—C)=aN 
ere C and @ are constants, and N is number of 


les. 
_In specimens stressed considerably above 


» fatigue limit, the number of cracks, M, in a 


narrow circumferential band of the specimen can 
be expressed by the equation 


/M=KS, 


where K is a constant, and S is the average stress 
in the band. 

10. The development of a small crack in a 
specimen resulted in a large decrease in the value 
of the fatigue limit as compared with specimens 
which had been stressed just short of cracking. 
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urification, Purity, and Freezing Points of 31 


Hydrocarbons of the API-NBS Series’ 


By Augustus R. Glasgow, Jr., Evelyn T. Murphy, Charles B. Willingham, and 
Frederick D. Rossini 


This report describes the purification and determination of freezing points and purity 
of 31 hydrocarbons of the API-NBS series, including 2 pentanes, 5 hexanes, 2 heptanes, 
1 octane, 3 alkyleyclopentanes, 3 alkyleyclohexanes, and 15 alkylbenzenes. 


I. Introduction 


In May 1943, the Advisory Committee for the 
verican Petroleum Institute Hydrocarbon Re- 
arch Committee, operating the API Hydro- 
ubon Research Project* at the Ohio State 
uiversity, requested the American Petroleum 

Research Project 6 at the National 
ureau of Standards to begin purification of a 
res of hydrocarbon compounds of the highest 
ucticable purity, under the label API-NBS 
drocarbons. This work was begun July 1, 
43, with a large number of compounds of 
atively high purity being supplied by the API 


This investigation was performed at the National Bureau of Standards 
part of the work of the American Petroleum Institute Research Project 6 
tbe Analysis, Purification, and Properties of Hydrocarbons. 
Research Associate on the American Petroleum Institute Research Project 
the National Bureau of Standards. 
On September 1, 1944, this project became the American Petroleum Re- 
h Project 45 on the “Synthesis and Properties of Hydrocarbons of Low- 
lecular Weight”, with Ceeil E. Boord continuing as Supervisor and Rob- 
F. Marschner becoming chairman of the Advisory Committee. D. P. 
ard was chairman of the former API Hydrocarbon Research Committee, 
hich Was merged on September 1, 1944, with the API Advisory Committee 
Fundamental Research on the Composition and Properties of Petroleum, 
¥hich J. Bennett Hill was chairman. 


rification of Hydrocarbons 


Hydrocarbon Research Project for further puri- 
fication. During the first year of operation, the 
purification of 31 compounds of the API-NBS 
series was completed. This report describes the 
purification and determination of freezing points 
and purity of these hydrocarbons, which are 
being made available on loan to qualified inves- 
tigators for the measurement of needed properties.* 


4 The allocation of these samples, by loan to qualified investigators for the 
measurement of needed properties, is handled by the Advisory Committee 
for the API Research Project 44 at the National Bureau of Standards on the 
“Collection, Analysis, and Calculation of Data on the Properties of Hydro- 
carbons” (W. E. Kuhn, chairman, Otto Beeck, Gustav Egloff, and 8. 
Kurtz, Jr., with F. D. Rossini as Supervisor of the Project). 
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The starting materials were supplied as follows: 

By the API Hydrocarbon Research Project, 
now the API Research Project 45 on the “Syn- 
thesis and Properties of Hydrocarbons of Low- 
Molecular Weight’, at the Ohio State University, 
under the supervision of Cecil E. Boord: Cyclo- 
pentane, ethylcyclopentane (A), methylcyclohex- 


ane (approximately half), ethyleyclohexane, 
m-xylene (A), p-xylene (A), 1,3,5-trimethylbenzene, 
n-butylbenzene (A), isobutylbenzene (A), sec- 
butylbenzene (A), and tert-butylbenzene (A). 

By the Barrett Division of the Allied Chemical 
and Dye Corporation: Cyclohexane, methylcyclo- 
hexane (approximately half). 

By the Standard Oil Co. (Indiana) and the 
M. W. Kellogg Co.: 2,3-Dimethylbutane. 

By the General Motors Corporation: 2,2,3- 
Trimethylbutane. 

By the Houdry Process Corporation: Methyl- 
cyclopentane. 

By the Humble Oil & Refining Co.: Toluene. 

By the Monsanto Chemical Co.: Ethylbenzene, 
isopropylbenzene. 

By the Standard Oil 
o-Xylene (A). 

By the Dow Chemical Co.: n-Propylbenzene 
(A). 

By the API Research Project 6: n-Pentane, 
isopentane, n-hexane, 2-methylpentane, 3-methy]l- 
pentane, 2,2-dimethylbutane, n-heptane, 2,2,4- 
trimethylpentane, benzene, 1,2,3-trimethylben- 
zene (A), 1,2,4-trimethylbenzene (A). In this 
group, the stocks of the following compounds were 
obtained by purchase of commercially available 
materials: n-Pentane, isopentane, and 2,2-di- 


Development Co.: 


4 A letter A following the name of a compound indicates that subsequently 
a new and usually slightly purer sample (B) has been or is being prepared, 
the description of which will appear in a Jater report. 





II. Materials and Purification 








methylbutane from the Phillips Petroleum (, 
n-heptane from the Westvaco Chlorine Prody 
Co.; 2,2,4-trimethylpentane from the Rohm 4 
Hass Co.; benzene, from the Koppers Co. 

Information regarding the volume of the stay, 
ing material, details of the purification by disti 
tion, and volume of the selected “best’’ lot, ; 
given in table 1. Additional details of the disti 
ing columns are given in preceding reports [1, 2) 
In addition to the purification by distillation, ,@ °"”"~ 
benzene, m-xylene, and p-xylene were subject" 
to a purification by crystallization with centrifug 
ing by B. J. Mair and A. J. Streiff of this ls) 
oratory [7]. 

For this group of 31 compounds, time was no 
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Pentane .. 
Methylbut 


Methylpen 


-Dimethy| 
Dimethy! 
Heptane.. 
2 3-Trimetl 


- a 4 24-Trimetl 
available for the detailed examination of the p 
rity, by measurement of freezing points, of tha“ 
distillate as a function of its volume,® and the besf{iifyciopentam 


sample was taken as an appropriately selecte 
“heart” cut (not necessarily the middle part ¢ 
the distillate) having the volume given in the las 
column of table 1. This volume may be com 
pared for each compound with the volume of th 
starting material given in the third column. 1) 
heart cut was selected largely on the basis o 
refractive index and boiling point, although, 
will be shown in the next report on the purifi 


ethyleycloy 
thyleyclope 
ALKYLCY 
yclohexane 
ethyleyclol 
thyleyclohe 


ALKYL 


pnzene .. 
pluene 





tion of API-NBS hydrocarbons [3], this is q@ivivenzene 

procedure which frequently results in the discard rae 

ing of the best material. Each final lot wagiptyiene (4) 

filtered through silica gel (4] to remove water an ~_ 

nonhydrocarbon impurities. -Trimeth 

> B+Trimeth: 

*Figures in brackets indicate the literature references at the end of O33}.5-Trimeth, 

paper. Butylbenze) 

*In the next report on the purification of API-NBS and API-Stané Butylbenz 

hydrocarbons, in connection with the cooperative program of the Natio -Butylbenz 

Bureau of Standards and the American Petroleum Institute on Stand -Butylbenz 
Samples of hydrocarbons, a procedure will be described wherein enough [re 

ing points are measured tc *~termine the purity of the distillate as a functs i 

of its volume (3). *Aletter A 

mew and ust 

deseriptior 

* The abbre 

APIRP 45; 

American 

yothesis an 

Ohio State 

APIRP 6; A 

, purificat i 

Standards, 

Parrett; Ba 

™ York, N. 


Dow; Dow ( 
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TaBLeE 1.—Information on the starting materials and purification 





Distillation ¢ 





Starting material pro- 


Compound * vided by laboratory > 


Number 
of theo- 
retical 


Azeotrope-forming plates at 
total 


substance if used ! (approx.) 


reflux 
(approx.) 





APIRP 6 
APIRP 6 


APIRP 6 
APIRP 6 
APIRP 6 


APIRP 6 
Std. (Ind.); Kellogg-.-. 


General Motors 
oy, 


APIRP 45 
Houdry 
APIRP 45. ._.... 


Barrett, APIRP 45__.. 


ethylcyclohexane 
: APIRP 45 


thyleyclohexane . . . - 


ALKYLBENZENES 


Mp eee Gk co cessed 
Dow....- 

Monsanto 

APIRP 6 
B.+-Trimethylbenzene (A)....| APIRP 6...--.- 
3+-Trimethylbenzene_.......| APIRP 45 
Butylbenzene (A) APERP @......<..... 
Butylbenzene (A)__. APIRP 45 
-Butylbenzene (A) of 
-Butylbenzene (A) APIRP 45 











Ethanol 
Cellosolve 


100/1 
100/1 
100/1 
100/1 
100/1 
100/1 
100/1 
100/1 
165/1 
165/1 
145/1 
145/1 
145/1 
145/1 
145/1 


Cellosolve___.... as 
.| Methyl Carbitol_..... 





eons eo wto ~~ “es 
—) & 

















*A letter A following the name of a compound indicates that subsequently 

new and usually slightly purer sample (B) has been or is being prepared, 
description of which will appear in a later report. 

* The abbreviations represent the following laboratories: 

APIRP 45; American Petroleum Institute Research Project 45 (formerly 
American Petroleum Institute Hydrocarbon Research Project) on the 
yothesis and Properties of Hydrocarbons of Low Molecular Weight,” at 
Ohio State University, Columbus, Ohio; C. E. Boord, Supervisor. 

APIRP 6; American Petroleum Institute Research Project 6 on the “Anal- 
®, purification, and properties of hydrocarbons”, at the National Bureau 
Standards, Washington, D. C.; F. D. Rossini, Supervisor. 

wrett; Barrett Division of the Allied Chemical & Dye Corporation, 


General Motors; General Motors Corporation, Detroit, Mich. 
Houdry; Houdry Process Corporation, Marcus Hook, Pa. 
Humble; Humble Oil & Refining Co., Houston, Texas. 
Kellogg; M. W. Kellogg Co., New York, N. Y. 

Monsanto; Monsanto Chemical Co., Dayton, Ohio. 

Std. Oil Dev.; Standard Oi] Development Co., Elizabeth, N. J. 
Std. (Ind.); Standard Oil Co. (Indiana), Whiting, Ind. 

© See references {1} and [2] for further details. 

¢ Three separate distillations of this volume were made. 

© See text regarding additional purification by crystallization. 

t Methyl Cellosolve is ethylene glycol monomethy! ether; Cellosolve is 


York, N.Y. ethylene glycol monoethyl ether; methyl Carbitol is diethylene glycol 
ow; Dow Chemical Co., Midland, Mich. monomethy] ether. 





III. Freezing Points and Purity 


time-temperature freezing or melting curves ¢ 
tained with the apparatus and procedure recep 
described [5]. From these data, the freezing poiy 
for zero impurity was determined according to th 


For each of the compounds except 3-methyl- 
pentane, which failed to yield crystals, the freezing 
point of the actual sample was determined from 


TABLE 2.—Freezing points and purity of 31 API-NBS 


hydrocarbons 








Kind of 
time- 
tempera- 
ture ob- 








ry 
tVaUuous 


used to 





the freez- 
ing point* 


phe mee 

int o 

Ls 
air at 1 
atm.> 





2-Methylpentane 
3-Methy!pentane 
2,2-Dimethylbutane 


2, 2,3-Trimethylbutane - - 
2, 2,4-Trimethylpentane_. 


ALE YLCYCLOPENTANES 


Ethylcyclopentane (A) - - 


ALKYLCYCLOHBXANES 


1,3,5-Trimethylbenzene.. 
n-Butylbenzene (A) 
iso-Butylbenzene (A).... 
sec-Butylbenzene (A)... 
tert-Butylhenzene (A) __. 





M 
M 





+5. 533+. 010 
—%4. 9914. 012 
— 94. 950-4. 025; 
—~ 25. 175. 010) 
—47. 872+. 020 
+13. 263. 012 


—4.7301)}.............-- 


— 87.904 
— 51.527 
— 75.572 
— 57.872 








- 0387+. 020 
-05 +. 03 
-16 +.09 
- 0824. 027 
15 +. 05 
-06 +. 03 
29 +. 08 
-04 +. 03 


-20 +. 04 


-32 4.20 
-04 +. 02 
- 09 +. 08 
-4 +.09 
31 +. 06 
-05 +. 03 





* F indicates freezing and M indicates melting. See reference [5] for details. 

* A roman numeral I indicates that the value is for that crystalline form 
having the highest freezing point. 

¢ The values in this column were calculated as described in reference [5], 
using the values of the cryoscopic constants and freezing points for rero im- 
purity given in the ‘z” tables of the American Petroleum Institute Research 
Project 44, reference [6]. For the 19 compounds for which values of the freez- 
ing point for zero impurity are reported in this paper, the values given in 
reference (6) are, at the time of this writing, identical to those in this table. 

¢ The amount of impurity in the 3-methylpentane was estimated by an- 
alogy with 2methylpentane. 
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procedure described in reference [5]. 

For benzene and p-xylene, in which the solu}jjj 
of water is quite appreciable, and for which ¢ 
freezing points are above that of water, the samp 
was introduced, by filtration through a small ty) 
containing silica gel, directly into the freezing ty) 
with substantially no contact with the air of i 
room after water in the sample was removed } 
the silica gel. Dry air was kept flowing into ¢ 
top of the freezing tube during the experiment. 

The purity of the actual sample was calculats 
according to the method described in section VI] 
of reference [5], using the values of the freeziy 
point of the actual sample and the freezing pois 
for zero impurity determined in this investigation 
(except as otherwise indicated), together wi 
values of the cryoscopic constants obtained fro 
the “z” tables of the American Petroleum Ins 
tute Research Project 44 [6] or from unpublisk 
data from references [3] and [8]. For each of th 
compounds except 3-methylpentane, table 2 ind 
cates whether freezing or melting experiments we 
used, and gives the values of the freezing point 
the actual sample, the freezing point for zero in 
purity where determined in this investigation, a 
the calculated amount of impurity. The amow 
of impurity in the 3-methylpentane was estimate 
by analogy with 2-methylpentane. From 
manner of preparation and purification of the 
compounds, it is believed that the impurities! 
most cases are substantially all close-boilis 
isomers. 


Grateful acknowledgment is made to the follo 
ing persons and laboratories for supplying 
materials for purification listed in table 1: C. i 
Boord, Supervisor of the American Petrolet 
Institute Research Project 45 (formerly the A! 
Hydrocarbon Research Project) at the Ohio St 
University, Columbus, Ohio; W. J. Sweet 
Standard Oil Development Co., New York, N. ! 
T. A. Boyd and W. G. Lovell, General Mote 
Corporation, Detroit, Mich. ; E. A. Smith, Hout 
Process Corporation, Marcus Hook, Pa.; B.! 
Dreisbach, Dow Chemical Co., Midland, Mit 
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rett Division of the Allied Chemical & Dye 
rporation, New York, N. Y.; Standard Oil Co. 
diana), Whiting, Ind.; M. W. Kellogg Co., 
w York, N. Y.; Monsanto Chemical Co., 
yton, Ohio; Humble Oil & Refining Co., Hous- 


Especial acknowledgment is due George Calin- 
gaert, formerly chairman of the subcommittee on 
purification and properties under the API Hydro- 
carbon Research Project, who made the arrange- 
ments for starting the series of API-NBS hydro- 
carbons. 
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hsorption Measurements for Broad Beams of 1- 


and 2-Million-Volt X-Rays 


By George Singer, Carl B. Braestrup*, and Harold O. Wyckoff 


Experimentally determined concrete-absorption curves are given for wide-angle X-ray 


beams produced by a 1- and 2-million-volt resonance generator. 


Data on the variation 


of X-ray intensity with the distance from a concrete barrier are included. With this 
information as a basis, the thickness of a concrete barrier necessary for protection against 
a broad X-ray beam has been computed for 1- and 2-million-volt X-rays and for various 


tube currents and barrier-to-target distances. 


I. Introduction 


Dne- and two-million-volt X-radiation is now 
1 routinely for industrial radiography. In 
application the cone of radiation used is much 
ader than is common in other applications. 
is is so for two reasons. The first is that the 
of a large field effects a saving in the time 
ured to radiograph a large object or a number 
small objects in one lot. The second is that 
massive cones needed to restrict the field 
very penetrating radiation must necessarily 
it the mobility of the X-ray generator; their 
is therefore commonly avoided, although 
rable for other reasons. 
onerete, being structurally self-supporting, 
tively inexpensive, and easily fabricated, has 
ome the standard material for the construction 
harriers for protection against very penetrating 
ation. Theoretical absorption curves appli- 
ble to such barriers are not currently available 
wide-angle beams. In the calculations for 
i absorption curves a correction must be made 
scattering which retains the scattered photons 


ding to reduce the attenuation for wide-angle 
is in comparison with that for narrow beams 
de same radiation. 


uot Physicist, Department of Hospitals, City of New York. 


Absorption of Broad Beams 


This problem has been recognized in several 
papers relating to X-ray protection.'?** How- 
ever, there are considerable technical difficulties 
involved in obtaining such data experimentally. 
For this reason all discussion in the literature 
relating to the increase in barrier thickness needed 
for wide-angle beams over that for narrow beams 
is of a general and indirect nature, and is not 
sufficiently detailed for specific application in the 
development of protection specifications. 

In the design of protective barriers, information 
is also needed on the variation of X-ray-dosage 
rate with distance from the protective barrier 
to the position at which personnel will be sta- 
tioned. Because of scattering from the barrier 
it is to be expected that for very penetrating radia- 

1 R. Jaeger and A. Trost, Elektrotech. Z. 61, 1,025 (1940). 

?T. BR. Folson and E. F. Focht, Am. J. Roentgenol. Radium Therapy 
51, 76 (1944). 

3 G. Singer and George C. Laurence, Ind. Radiology 4@, Fall (1945); Radiol- 


ogy 46, 69 (1946). 
‘GQ. Failla, Am. J. Roentgenol. 54, 553 (1945). 
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tion this variation will not follow the inverse- 
square law for the distance measured from the 
tube target. It is to be expected that this 
variation will be more rapid than indicated by 
this law because both the X-ray-tube target and 
the irradiated volume of concrete act as sources 
of radiation. 

It is the purpose of this paper to present (1) 


experimental data on absorption of X-y 
concrete for wide-angle beams generated by 
ages of 1 and 2 million volts, (2) prelimi 
data on the variation of the dosage rate with 
distance from a protective barrier, and (3) 
on concrete protection requirements for |. 
2-million-volt radiation under certain condi 
which obtain in industry. 


II. Method 


The experimental arrangement used in obtaining 
the basic data given below is indicated in figure 1. 
The concrete wall, C, is the main wall of the 
building housing the generator and is 41 inches * 
thick; D, is the covered maze within which the 
ionization measurements were made. 

Dosage measurements were made by calibrated 
dosage instruments designed by Braestrup.® 
Three sizes of ionization chambers were used; 
their characteristics are listed in table 1. The 
ionization chamber potentials were measured 
before and after irradiation by means of a string 
electrometer which is an integral part of the dosage 
meter. All ionization chambers were calibrated 


§ Because the thickness of protective barriers is often of practical interest 
to persons unfamiliar with metric units, we have specified the dimensions 
of all barriers in feet and inches. 

*C. B. Braestrup, Radiology 38, 22, 207 (1942). 


SCALE 
‘apres | 
Psenea| 


Ficure 1.—Ezperimental arrangement showing two sets of 
concrete test blocks in place. 


A total of three sets were used. 


with both 175-kv X-rays and gamma-rays { 
radium. Both calibrations agreed very 

with the theoretical values based upon th: 
volume and capacitance. At the National 
reau of Standards the small thimble chan 
and the electrometer head were checked a 
a similar instrument combination with 1-nilj 
volt radiation and were found to be in agreen 
therewith to within 2 percent. 


TABLE 1.—Jonization chambers 





Inside dimensions 
Wall mate- 
rial 





Diameter} 





em 
CX-1, Cylindri- | Bakelite__.. 12.05 

eal. 
SX-11, Spherical...) Polystyrene. 10.5 
D-1, Thimble 0.7 

















1} Dose in roentgens per scale division of electrometer. 


The barrier material consisted of concrete bi 
of density 130 lb/ft *, measuring 5% by 7% by 
in. These were laid in the doorway, A, sho 
figure 1. Measurements were made through 
of double layers of these blocks; one 5% in. 
the other 7% in. thick. The layers were stag 
so that the planes of the cracks between | 
did not pass through either the tube targe 
ionization chamber. The opening remaining 
tween the top of the test barrier and the lint 
the doorway to the maze was filled with lead b 
For this arrangement both the radiation le 
through the doorway, B—produced by the 


barrier were found to be negligible in compat 
with the dosage rate at the position of the 
ionization chambers. 








The distance, Z, in figure 1 was 1.83 m in all cases; 
js the barrier thickness. This thickness was 
nverted to thickness of concrete of a density of 
7 lb/ft * before plotting. The wall-to-center-of- 
umber distance, G, varied from 2 to 4 in. for all 
otted points. The ionization chambers were 
laced as near as practicable to the radiation 
ier. It will be shown later that this position 
es the minimum apparent absorption. 
The X-ray generator used is the 2-million-volt 
-ray resonance generator described by Charlton 


and Westendorp,’ which is operative at both 1 and 
2 million volts. All measurements were made 
with X-radiation transmitted through the tube 
target along the axis of the tube. The anode cone 


consisted of a cylinder with a lead wall 5 cm thick, 
the end of which was shaped so as to confine the 
X-ray beam to a cone of radiation measuring 
approximately 70° at its apex. 


’ E, E. Chariton and W. F. Westendorp, Electronics 17, 128 (1944). We 
are indebted to the General Electric Co. and to E. E. Chariton, of their 
Research Laboratory, for facilities that were made available to us in carrying 
out these experiments. 


Ill. Results 


The attenuation curves obtained are shown in 
re 2. As a check on the reliability of the 
thod used, transmission measurements for 
nillion-volt radiation through a solid concrete 

42 in. thick of density 152 Ib/ft* were also 














is 20 25 30 35 
CONCRETE THICKNESS (UNCHES) 


vrE 2.—Concrete absorption curves for 70° beams of 1- 
and 2-million-volt X-rays. 


Measurements were made with the radiation transmitted through a 

mgsten target along the tube axis. A resonance generator was used. 
voltage waveform of such generators is sinusoidal. All indicated 
és are peak or crest values, 
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made. The single point so obtained is plotted in 
figure 2 at a thickness of about 43 in. 

The effect of chamber-to-wall distance on the 
apparent attenuation is illustrated by the curves 
in figure 3. The ordinate is in arbitrary units. It 
is obtained by dividing the dose per milliampere 
minute by the square of the sum of the distances 
E, F, and G. The abscissa is the chamber-wall 
distance, G, shown in figure 1. If there were 
neither scattering in the wall nor absorption in the 
air, each set of plotted points should give a straight 
line with zero slope. Actually, because of wall 
scattering, the value of the ordinate decreases 
with increase of abscissa. These plotted points 
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WAL(-TO-CHAMBER DISTANCE (METERS) 


Fieurse 3.—Variation of X-ray dose with distance from a 

concrete protective barrier. 

The abscissa is the distance to the ionization chamber measured from the side 
of the barrier from which the radiation emerges. The ordinate is ob- 
tained by dividing the dose per milliampere minute by the square of the 
distance from the X-ray source to ionization chamber. 





fall along a straight line for each experimental lished by Folson and Focht.* Accordingly 
setup. The slopes measured on a linear scale— may be concluded that within these limits ; 
not the logarithmic scale of figure 3—are approxi- beam angle has no important effect on the chy 


mately 0.25. This is so not only for our own’ in apparent attenuation with chamber-to. 
l-million-volt and 2-million-volt data for 70° distance. 
beams but also for the 1 data on a 4° beam pub- * Folson and Focht, Am. J. Roentgenol, Radium Therapy 51, % s 


IV. Conclusions 


From the data presented above it is possible to _million-volt X-ray generator, 120 r/min and th 
compute for broad X-ray beams the thickness of _ the person to be protected will normally be sg. m= 
concrete barriers necessary for protection against tioned near the protective barrier. gh 
1- and 2-million-volt X-radiation from targets of TABLE 2.—Thickness of concrete of density 147 Wy 


the ergum 
the transmission type. Such targets are standard required for protection against broad beams of X eel 


in this quality range. In table 2 the thickness of for the current, distance, and voltage indicated 2.5000 t 
such barriers is listed for 1- and 2-million-volt . a te 


X-radiation and for the various tube currents and Target ee ae 2million-volt X-rays TAR 5.00 t0 1 


target distances normally encountered in practice. a hs eee abies ee 0 malae mane: yim 
The tubular values given therein were computed | sun An 
directly from the absorption curves of figure 2. Oe el Geek aeee Geolamel aol eal ~“ 
A few points at short distances and high currents 6 | m5] 3.5) 335] 345] 43. 0 857M Tames 


. ° . . . 6 23.0; 30.0} 32.5 . ° i e 
required linear extrapolation for their determi- 3 (| east asl as / a 


nation. Absorption of X-1adiation by air has been 10 25.0] 27.0] 22.5) 30. am, $2, 
neglected. A permissible daily dose of 0.1 r has ne ae — 


been assumed.» It was further assumed that 21.0} 23.0] 25.0 5] 38 | mn. (19 
the dosage rate per milliampere at 1m for a 1-mil- Sat sat ie : oe 


lion-volt X-ray generator is 20 r/min and for a 2- 15.5 19.5 | 21. ale of 
d in bue! 
75 13.5 17.5 . . y .0) BS T. TABLE | 
100 11.5 15.5 
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*X-ray protection, National Bureau of Standards Handbook HB20 
(1936). 


” Safety code for the industrial use of X-rays, American Standards Asso- 
ciation, New York, N. Y. (1946). Wasaineaton, May 31, 1946 jens to 
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Tion To Inrencat Equations of Porentiat Tarorny. By Chester Snow, 
National Bureau of Standards. Reproduced from original handwritten 
manuscript. (1942) VII+319 pages; bound in heavy paper cover. Out 
of print. 
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1 MT11—$2 ($2.50) 
C) MT12—$2 ($2.50) 
C) MT13—$2 ($2.50) 
C) MT14—$2 ($2.50) 
C) MT16—$2 ($2.50) 





(1) MT17—$1.50($1.75) [) MT25—$0.25 ($0.30) 
C) MT18—$0.25 ($0.30) [] MT26—$0.25 ($0.30) 
C1) MT19—$0.25 ($0.30) [] MT27—$0.25 ($0.30) 
C) MT20—$0.25 ($0.30) [) MT28—$0.25 ($0.30) 
C) MT21—$0.25 ($0.30) [) MT29—$0.25 ($0.30) 
(C1) MT22—$0.25 ($0.30) [] MT30—$0.25 ($0.30) 
CO) MT24—$0.25 ($0.30) [] MT31—$0.25 ($0.30) 


1 MT32—$0.25 ($0.30) 
C1) MT33—$0.25 ($0.30) 
1 MT34—$0.25 ($0.30) 
C) MT35—$0.25 ($0.30) 
0) MT36—$0.25 ($0.30) 





(Name) 





STANDARDS” IN 
8. CURRENCY. 


(Address) 








(City, Zone, and State) 





MATHEMATICAL TABLES—Continued 


MT17. Miscettansous Parsicat Tasies: Prance’s Rapiation Functions 
(Originally published in the Journal of the Optical Society of America, 
February 1940); and Execrromic Functions. (1941) VII+58 pages; 
bound in buckram, $1.50 ($1.75). 

MT18. Taste or tHe Zenos oy THe Lecenpae Potrnomuats oF Onper 
1-16 ano ree Weicer Cozrricients ror Gauss’ Mecuamcat Quap- 
nature Foamuta. (Reprinted from Bul. Amer. Mathematical Society, 
October 1942.) 5 pages, with cover, 25 cents (30#). 

MT19, On tae Forcrnon H (m, a, x) <axp (—ix)F (m+1—ia, 2m+2; 2ix); 
with table of the confluent hypergometric function and its first derivative. 
(Reprinted from J. Math. Phys, December 1942.) 20 pages, with cover. 
25 cents (30#). 


MT20. Tanus ov Inrmcnars [* Jo(s)dt and [,"Yolehde. (Reprinted from J. 


Math. Phys., May 1943.) 12 pages, with cover, 25 cents (30¢). 

MT21. Tamux ov Jisie) = f~ 2 4 4x Rexatep Funcnions. (Reprinted 
from J. Math. Phys., June 1943.) 7 pages, with cover, 25 cents (30¢). 

MT22. Taste oF Corrricients mm Numeeicat Inrecration Formutan. 
(Reprinted from J. Math. Phys., June 1943.) 2 pages, with cover, 25 
cents (30¢). 

MT23. Tape or Founrer Coerricixnts. (Reprinted from J. Math. Phys. 
Sept. 1943.) 11 pages, with cover. (Out of print.) 

MT24. Cogrriciznts ror Nuwenica, Dirvexentunion Wrre Centaat 
Durrerences. (Reprinted from J. Math. Phys., Sept. 1943.) 21 pages, 
with cover, 25 cents (30¢). 

MT25. Seven-Powrr Lacranciuan Inrecration Foamvias. 
from J. Math. Phys., Dec. 1943.) 4 pages, with cover, 25 cents (30¢). 

MT26. A Snont Tasce or tue Finer Five Zenos oF rue Transcenventar 
Equation. Jo(x) Yo(ke)—Jo(kx) Yo(x)=0. (Reprinted from J. Mathy 
Phys., Dec. 1943.) 2 pages, with cover, 25 cents (30¢). 


pages, with cover, 25 cents (30¢). 

MT28. Tamu oF falx)=(ZgyaJe(x). (Reprinted from J. Math. Phy, 
1944.) 16 pages, with cover, 25 cents (30¢). 

MT29. Taste of Coerricients ror Inverse Inteerotation Wm 
vancinc Dirrenences. (Reprinted from J. Math. Phys. May 
28 pages, with cover, 25 cents (30¢). 

MT30. A New Foruuts ror Inverse Internrotation. (Reprinted 
Bul. Amer. Mathematical Society, Aug. 1944.) 4 pages, with cove, 
cents (30¢). 

MTS31. Cogrricierts ror Inrerrotation Wrrain a Square Caw » 
Comptex Pians. (Reprinted from J. Math. Phys., Aug. 1944.) 1) 


(Reprinted from J. Math. Phys, 
1945.) 3 pages, with cover, 25 cents (30¢). 


vernences, vor Iwrecratinc Ssconp Ornper Dirrenentiat Fou. 
(Reprinted from J. Math. Phys, Nov. 1945.) 6 pages, with cove, 
cents (30¢). 

MT36. Foamvu.as ror Dmecr anv Inverse InreRro.ation oF 4 Cc 





GOVER 


